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Abstract 

Prolonging  the  lifetime  of  energy  constrained  wireless  sensor  networks  is  a  crucial  challenge  in 
sensor  network  research.  In  this  paper,  we  present  a  novel  approach  based  on  fuzzy  logic  theory 
to  analyze  the  lifetime  performance  of  a  sensor  network.  We  demonstrate  that  a  type-2  fuzzy 
membership  function(MF),  i.e.,  a  Gaussian  MF  with  uncertain  variance  is  most  appropriate 
to  model  a  single  node  lifetime  in  wireless  sensor  networks.  In  our  research,  we  concern  two 
basic  sensor  placement  schemes:  square-grid  and  hex-grid.  Two  fuzzy  logic  systems (FLSs):  a 
singleton  type-1  FLS  and  an  interval  type-2  FLS  are  designed  to  perform  lifetime  estimation  of 
the  entire  sensor  network.  Simulation  results  show  that  the  interval  type-2  FLS  in  which  the 
antecedent  membership  functions  are  modeled  as  type-2  fuzzy  sets  outperforms  the  singleton 
type-1  FLS  when  the  single  node  lifetime  behaves  the  nature  of  Gaussian  MFs  with  uncertain 
standard  deviation. 


1  Introduction 

Wireless  sensor  network  represents  significant  advantages  over  conventional  personnel-rich  meth¬ 
ods  of  long  term  data  collection  and  monitoring.  Starting  in  the  1980s,  networked  microsensors 
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technology  has  been  widely  used  in  military  applications.  Examples  of  these  applications  include 
the  Cooperative  Engagement  Capability  (CEC),  Remote  Battlefield  Sensor  System  (REMBASS), 
Advanced  Deployable  System  (ADS)  and  Tactical  Remote  Sensor  System  (TRSS)  [1].  Although 
a  majority  of  military  applications  deploy  sensor  nodes  randomly,  a  number  of  other  applications 
require  the  placement  of  sensors  at  desired  locations.  Such  placement-friendly  sensor  networks  are 
developed  rapidly  for  infrastructure  security,  environment  and  habitat  monitoring,  traffic  control 
etc  [1].  A  real-world  32  node  habitat  monitoring  sensor  network  system  was  deployed  on  a  small 
island  off  the  coast  of  Maine  to  study  nesting  patterns  of  Petrels  [2].  Other  applications  involve  the 
use  of  sensors  in  buildings  for  environmental  monitoring  which  may  include  chemical  sensing  and 
detection  of  moisture  problems.  Structural  monitoring  [3]  and  inventory  control  are  some  other 
applications  of  such  networks. 

Wireless  sensor  network  consists  of  certain  amount  of  small  and  energy  constrained  nodes.  Basic 
components  of  sensor  node  include  a  single  or  multiple  sensor  modules,  a  wireless  transmitter- 
receiver  module,  a  computational  module  and  a  power  supply  module.  Such  networks  are  normally 
deployed  for  data  collection  where  human  intervention  after  deployment,  to  recharge  or  replace 
node  batteries  may  not  be  feasible,  resulting  in  limited  network  lifetime.  Most  applications  have 
pre-specified  lifetime  requirements,  for  instance  the  application  in  [2]  has  a  lifetime  requirement 
of  at  least  9  months.  Thus  estimation  of  lifetime  of  such  networks  prior  to  deployment  becomes 
a  necessity.  Prior  works  on  evaluating  lifetime  have  considered  networks  where  sensor  nodes  are 
randomly  deployed.  [4]  gives  the  upper  bound  on  lifetime  that  any  network  with  the  specified 
number  of  randomly  deployed  nodes,  source  behavior  and  energy  can  reach  while  [6]  discusses  the 
upper  bounds  on  lifetime  of  networks  with  cooperative  cell  based  strategies.  Network  liftime  of 
fixed  deployment  schemes  are  recently  studied  in  [7].  Jain  and  Liang  [7]  discovers  that  in  wireless 
sensor  networks,  a  single  node  lifetime  behaves  the  nature  of  normal  distribution  which  brings  the 
first  light  of  exploring  the  network  lifetime  behavior  given  the  knowledge  of  nodes  lifetimes. 

In  this  paper,  we  deal  with  the  issue  of  lifetime  analysis  and  estimation  for  wireless  sensor 
networks  in  which  the  sensor  nodes  are  deployed  at  desired  locations.  We  propose  to  apply  interval 
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type-2  fuzzy  logic  systems  (FLSs)  for  network  lifetime  analysis  and  estimation.  Our  approach 
is  entirely  different  from  all  prior  research.  We  demonstrate  that  a  type-2  fuzzy  membership 
function  (MF),  i.e.,  a  Gaussian  MF  with  uncertain  variance  is  most  appropriate  to  model  a  single 
node  lifetime  in  wireless  sensor  networks.  Two  fuzzy  logic  systems(FLSs):  a  singleton  type-1  FLS 
and  an  interval  type-2  FLS  are  constructed  to  perform  lifetime  estimation  of  the  entire  sensor 
network.  In  our  work,  we  concern  two  basic  sensor  placement  schemes:  square-grid  and  hex-grid. 
We  believe  that  these  two  schemes  can  serve  as  basis  for  evaluating  more  complex  schemes  for  their 
lifetime  performance  prior  to  deployment  and  help  justify  the  deployment  costs. 

The  rest  of  this  paper  is  organized  as  follows.  In  section  2,  we  detail  two  sensor  deployment 
schemes  used  for  this  study  and  the  basic  concepts  on  network  coverage,  connectivity  and  lifetime. 
Section  ??  gives  an  overview  of  interval  type-2  fuzzy  logic  systmes.  In  section  4,  we  demonstrate 
that  a  single  node  lifetime  can  be  modeled  with  Gaussian  MFs.  We  found  that  a  type-2  Gaus¬ 
sian  MF  with  uncertain  variance  is  most  appropriate  to  model  a  single  node  lifetime  in  wireless 
sensor  networks.  We  apply  this  knowledge  and  design  an  interval  type-2  FLS  in  section  5  to  ana¬ 
lyze  network  lifetime.  A  singleton  type-1  FLS  is  constructed  as  well  for  performance  comparison. 
Simulation  results  and  discussions  are  also  presented  in  section  5.  Section  6  concludes  this  paper. 

2  Preliminaries 

2.1  Basic  Model  and  Assumptions 

Consider  identical  wireless  sensor  nodes  placed  in  a  square  sensor  field  of  area  A.  All  nodes  are 
deployed  with  equal  energy.  Each  sensor  is  capable  of  sensing  events  up  to  a  radius  rs,  the  sensing 
range.  Communication  range  rc  is  defined  as  the  distance  beyond  which  the  transmitted  signal  is 
received  with  signal  to  noise  ratio  (SNR)  below  the  acceptable  threshold  level.  In  this  paper,  We 
assume  the  communication  range  rc  to  be  equal  to  the  sensing  range  rs. 

Direct  communication  between  two  sensor  nodes  is  possible  only  if  their  distance  of  separation 
r  is  such  that  r  <rc.  We  call  such  nodes  neighbors.  Communication  between  a  sensor  node  and  its 
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non-neighboring  node  is  achieved  via  peer-to-peer  communication.  Thus  the  maximum  allowable 
distance  between  two  nodes  who  wish  to  communicate  directly  is  rmax  —  rc  =  rs.  A  network  is  said 
to  be  deployed  with  minimum  density  when  the  distance  between  its  neighboring  nodes  is  r  =  rmax . 

2.2  Placement  Schemes 

The  simplest  placement  schemes  involve  regular  placement  of  nodes  such  that  each  node  in  the  net¬ 
work  has  the  same  number  of  neighbors.  We  arrive  at  two  basic  placement  schemes  by  considering 
cases  where  each  sensor  nodes  has  four  and  three  neighbors.  This  leads  us  to  the  square-grid  and 
hex-grid  placement  schemes  shown  in  Fig.  1(a)  and  (b)  respectively. 

A  sensor  node  placement  scheme  that  uses  two  neighbors  per  sensor  node  has  been  described 
in  [13].  We  believe  that  these  three  elementary  placement  schemes  can  serve  as  basis  for  other 
placement  schemes,  because  a  placement  scheme  of  any  complexity  can  be  decomposed  into  two- 
neighbor,  three-neighbor  and  four-neighbor  groups.  Both  grids  shown  have  the  same  number  of 
nodes1  and  nodes  in  both  grids  are  equidistant  from  their  respective  neighbors  (with  distance  of 
separation  r). 

2.3  Coverage  and  Connectivity 

Coverage  and  connectivity  are  two  important  performance  metrics  of  networks  and  hence  a  discus¬ 
sion  on  them  becomes  imperative  before  the  lifetime  of  the  network  can  be  defined. 

Coverage  scales  the  adequacy  with  which  the  network  covers  the  sensor  field.  A  sensor  with 
sensing  range  rs  is  said  to  cover  or  sense  a  circular  region  of  radius  rs  around  it.  If  every  point 
in  the  sensor  field  is  within  distance  rs  from  at  least  one  sensor  node,  then  the  network  is  said  to 
provide  complete  or  100%  coverage. 

Various  levels  of  coverage  are  acceptable  depending  on  the  application.  In  critical  applications, 

complete  coverage  is  required  at  all  times.  Any  loss  of  coverage  leads  to  a  sensing  gap  in  the  field. 

The  Hex-grid  has  lower  density  than  the  Square-grid.  With  36  nodes  deployed,  the  network  with  Square-grid 
covers  an  area  of  25 r2,  and  the  Hex-grid  covers  an  area  of  48 r2,  almost  double  that  of  the  square  grid. 
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Figure  1:  Placement  Schemes  for  a  36  node  sensor  network:  (a)  Square-Grid  (b)  Hex-Grid. 


Such  gaps  cause  breach  of  security  in  case  of  surveillance  applications.  Also,  in  applications  which 
require  data  with  high  precision  ,  a  sensing  gap  leads  to  inaccuracies.  For  such  networks  any  loss 
of  coverage  renders  the  network  nonfunctional.  In  some  other  applications  a  small  loss  of  coverage 
may  be  acceptable. 

Connectivity  scales  the  adequacy  with  which  nodes  are  able  to  communicate  with  their  peers. 
One  of  the  strengths  of  sensor  networks  arises  from  their  ability  to  aggregate  data  collected  from 
different  sensor  nodes.  This  requires  adequate  communication  between  sensor  nodes.  Any  node 
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should  be  able  to  communicate  with  any  other  node  for  proper  functioning  of  the  network.  If  a 
single  node  gets  isolated  due  to  failure  of  all  its  neighbors,  it  will  be  unable  to  communicate  with 
the  rest  of  the  network.  If  a  large  number  of  nodes  fail  due  to  lack  of  energy,  a  part  of  the  network 
may  get  completely  disconnected  from  the  rest.  If  a  large  number  of  nodes  fail  due  to  lack  of  energy, 
a  part  of  the  network  may  get  completely  disconnected  from  the  rest.  In  our  analysis  we  assume 
that  only  100%  connectivity  is  acceptable  and  the  network  fails  with  any  loss  of  connectivity.  An 
example  of  a  sensor  placement  scheme  that  concentrates  mainly  on  coverage  as  its  parameter  of 
interest  can  be  found  in  [14],  where  a  sensor  placement  algorithm  for  grid  coverage  has  been 
proposed. 

In  our  analysis  we  require  the  network  to  provide  complete  coverage  and  connectivity.  We  give 
equal  importance  to  both  parameters  and  declare  the  network  nonfunctional  if  either  of  them  falls 
below  their  desired  levels. 

2.4  Lifetime 

The  basic  definition  of  lifetime,  or  more  precisely  the  post-deployment  active  lifetime  of  a  network 
is  the  time  measured  from  deployment  until  network  failure.  Based  on  the  levels  of  coverage 
and  connectivity  required  to  deem  a  network  functional,  network  failure  can  be  interpreted  in 
different  ways.  Since  only  complete  coverage  and  connectivity  are  acceptable  to  us,  network  failure 
corresponds  to  the  first  loss  of  coverage  or  connectivity. 

In  this  paper,  we  concentrate  on  finding  the  minimum  lifetime  of  a  network,  the  worst  case 
scenario.  To  be  able  to  evaluate  this  minimum  lifetime,  we  need  to  know  the  lifetime  of  a  single 
sensor  node  the,  the  minimum  number  of  node  failures  that  cause  network  failure,  and  the  positional 
relationship  2  between  the  failed  nodes. 

2 Positional  relationship  between  two  nodes  can  be  that  the  two  nodes  axe  diagonal,  adjacent  or  completely  unre- 
lated. 
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2.4.1  Minimum  Density  Networks 


Consider  the  square-grid  and  the  hex-grid  deployed  with  minimum  density.  Both  schemes  survive 
the  failure  of  a  single  node  without  loss  of  either  connectivity  or  coverage  implying  that  the  minimum 
number  of  node  failures  that  can  lead  to  network  failure  is  greater  than  one.  Failure  of  any  two 
neighboring  nodes  causes  loss  of  coverage  and  hence  network  failure  as  indicated  in  Figs  2(a)  and 
(b). 

Thus  the  minimum  number  of  node  failures  that  cause  network  failure  is  two  and  these  two 
nodes  must  be  adjacent  to  each  other  (neighbors).  A  network  may  undergo  multiple  node  failures 
and  still  be  connected  and  covered  if  any  of  the  failed  nodes  are  not  neighbors.  But  the  absolute 
minimum  number  of  node  failures  that  can  cause  network  failure  is  two. 

2.4.2  Network  lifetime  for  dense  placement 

Let  Nmin  be  the  number  of  sensor  nodes  deployed  in  a  sensor  field  of  area  A  with  minimum  density. 
In  a  minimum  density  network,  the  separation  between  neighbors  is  rmax ,  the  maximum  range  of 
the  sensor  nodes.  A  network  deployed  with  higher  density  would  require  a  larger  number  of  nodes 
N,  N  >  Nmin  and  smaller  separation  between  neighboring  nodes  r,  r  <  rmax.  Since  the  range  of  a 
sensor  node  rmax  is  greater  than  the  distance  to  its  immediately  adjacent  nodes,  it  is  now  be  able  to 
communicate  not  only  with  the  nodes  immediately  adjacent  to  it,  but  also  with  other  nodes  which 
fall  in  its  range.  In  other  words,  each  node  now  has  more  number  of  neighbors  that  any  node  in  a 
minimum  density  network. 

We  discussed  earlier  that  the  minimum  lifetime  of  a  square-grid  or  hex-grid  network  with  r  — 
rmax (minimum  density)  is  the  time  to  failure  of  two  neighboring  nodes.  High  density  networks  can 
lose  ( N  -  Nmiri)  nodes  and  still  be  deployed  with  minimum  density  and  fail  only  after  the  further 
loss  of  a  minimum  of  two  neighboring  nodes.  Since  the  network  after  ( N  —  Nmin )  node  failures  is 
deployed  with  minimum  density,  the  time  taken  for  the  failure  of  any  two  neighboring  nodes  would 
be  equal  to  the  minimum  lifetime  of  a  minimum  density  network  (Tmd)  which  was  discussed  in  the 
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(b) 


Figure  2:  Loss  of  coverage  due  to  failure  of  two  neighboring  nodes:  (a)Square-grid:  Failure  of  nodes 
20  and  21  causes  loss  of  coverage.  (b)Hex-grid:  Failure  of  nodes  20  and  25  causes  loss  of  coverage. 

preceding  section.  If  t^ease,  the  time  taken  for  the  failure  of  (N  —  Nmin)  nodes  can  be  found,  we 
define  the  lifetime  of  a  high  density  network  (Thd)  as: 


Thd  =  td  ense  -^md 


(1) 


3  Introduction  to  Type-2  Fuzzy  Sets  and  Interval  Type-2  Fuzzy 
Logic  Systems 

3.1  Introduction  to  Type-2  Fuzzy  Sets 

The  concept  of  type-2  fuzzy  sets  was  introduced  by  Zadeh  [15]  as  an  extension  of  the  concept  of 
an  ordinary  fuzzy  set,  i.e.,  a  type-1  fuzzy  set.  Type-2  fuzzy  sets  have  grades  of  membership  that 
are  themselves  fuzzy  [17].  A  type-2  membership  grade  can  be  any  subset  in  [0,1]  -  the  primary 
membership ;  and,  corresponding  to  each  primary  membership,  there  is  a  secondary  membership 
(which  can  also  be  in  [0,1])  that  defines  the  possibilities  for  the  primary  membership.  A  type-1 
fuzzy  set  is  a  special  case  of  a  type-2  fuzzy  set;  its  secondary  membership  function  is  a  subset  with 
only  one  element,  unity.  Type-2  fuzzy  sets  allow  us  to  handle  linguistic  uncertainties,  as  typified 
by  the  adage  “words  can  mean  different  things  to  different  people.”  A  fuzzy  relation  of  higher  type 
(e-g.,  type-2)  has  been  regarded  as  one  way  to  increase  the  fuzziness  of  a  relation,  and,  according  to 
Hisdal,  “increased  fuzziness  in  a  description  means  increased  ability  to  handle  inexact  information 
in  a  logically  correct  manner  [18]”. 

Figure  3  shows  an  example  of  a  type-2  set.  The  domain  of  the  membership  grade  corresponding 
to  x  =  4  is  also  shown.  The  membership  grade  for  every  point  is  a  Gaussian  type-1  set  contained 
in  [0, 1],  we  call  such  a  set  a  “Gaussian  type-2  set”.  When  the  membership  grade  for  every  point  is 
a  crisp  set,  the  domain  of  which  is  an  interval  contained  in  [0, 1],  we  call  such  type-2  sets  “interval 
type-2  sets”  and  their  membership  grades  “interval  type-1  sets”.  Interval  type-2  sets  are  very  useful 
when  we  have  no  other  knowledge  about  secondary  memberships. 

3.2  Introduction  to  Type-2  Fuzzy  Logic  Systems:  An  Overview 

Figure  4  shows  the  structure  of  a  type-2  FLS.  It  is  very  similar  to  the  structure  of  a  type-1  FLS  [16]. 
For  a  type-1  FLS,  the  output  processing  block  only  contains  the  defuzzifier.  We  assume  that  the 
reader  is  familiar  with  type-1  FLSs,  so  that  here  we  focus  only  on  the  similarities  and  differences 
between  the  two  FLSs. 
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Figure  3.  (a)  Pictorial  representation  of  a  Gaussian  type-2  set.  The  secondary  memberships  in  this 
type-1  fuzzy  set  are  shown  in  (b),  and  are  Gaussian.  Note  that  this  set  is  called  a  Gaussian  type-2 
set  because  all  its  secondary  membership  functions  are  Gaussian.  The  “principal”  membership 
function  (the  bold  line),  which  is  triangular  in  this  case,  can  be  of  any  shape. 

The  fuzzifier  maps  the  crisp  input  into  a  fuzzy  set.  This  fuzzy  set  can,  in  general,  be  a  type-2 
set. 

In  the  type-1  case,  we  generally  have  “IF-THEN”  rules,  where  the  1th  rule  has  the  form  “Rl  :  IF 
xi  is  Fj  and  x2  is  Fj  and  •  •  •  and  xp  is  F*,,  THEN  y  is  Gl”,  where:  XjS  are  inputs;  F*s  are  antecedent 
sets  (i  —  1, . . .  ,p);  y  is  the  output;  and  G*s  are  consequent  sets.  The  distinction  between  type-1 
and  type-2  is  associated  with  the  nature  of  the  membership  functions,  which  is  not  important  while 
forming  rules,  hence,  the  structure  of  the  rules  remains  exactly  the  same  in  the  type-2  case,  the 
only  difference  being  that  now  some  or  all  of  the  sets  involved  are  of  type-2;  so,  the  1th  rule  in  a 
type-2  FLS  has  the  form  “Rl  :  IF  x\  is  F*  and  x2  is  fl2  and  •  •  •  and  xp  is  F*,,  THEN  y  is  G;”. 

In  the  type-2  case,  the  inference  process  is  very  similar  to  that  in  type-1.  The  inference  engine 
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TYPE-2  FUZZY  LOGIC  SYSTEM 


Figure  4:  The  structure  of  a  type-2  FLS.  In  order  to  emphasize  the  importance  of  the  type-reduced 
set,  we  have  shown  two  outputs  for  the  type-2  FLS,  the  type-reduced  set  and  the  crisp  defuzzified 
value. 

combines  rules  and  gives  a  mapping  from  input  type-2  fuzzy  sets  to  output  type-2  fuzzy  sets.  To 
do  this,  one  needs  to  find  unions  and  intersections  of  type-2  sets,  as  well  as  compositions  of  type-2 
relations. 

In  a  type-1  FLS,  the  defuzzifier  produces  a  crisp  output  from  the  fuzzy  set  that  is  the  output  of 
the  inference  engine,  i.e.,  a  type-0  (crisp)  output  is  obtained  from  a  type-1  set.  In  the  type-2  case, 
the  output  of  the  inference  engine  is  a  type-2  set;  so,  we  use  “extended  versions”  (using  Zadeh’s 
Extension  Principle  [15])  of  type-1  defuzzification  methods.  This  extended  defuzzification  gives  a 
type-1  fuzzy  set.  Since  this  operation  takes  us  from  the  type-2  output  sets  of  the  FLS  to  a  type-1 
set,  we  call  this  operation  “type-reduction”  and  the  type-reduced  set  so  obtained  a  “type-reduced 
set”. 

To  obtain  a  crisp  output  from  a  type-2  FLS,  we  can  defuzzify  the  type-reduced  set.  The  most 
natural  way  of  doing  this  seems  to  be  by  finding  the  centroid  of  the  type-reduced  set;  however, 
there  exist  other  possibilities  like  choosing  the  highest  membership  point  in  the  type-reduced  set. 

General  type-2  FLSs  are  computationally  intensive,  because  type-reduction  is  very  intensive. 
Things  simplify  a  lot  when  secondary  membership  functions  (MFs)  are  interval  sets  (in  this  case, 
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the  secondary  memberships  are  either  0  or  1).  When  the  secondary  MFs  are  interval  sets,  we  call 
the  type-2  FLSs  “interval  type-2  FLSs”.  In  [21],  we  proposed  the  theory  and  design  of  interval 
type-2  fuzzy  logic  systems  (FLSs).  We  proposed  an  efficient  and  simplified  method  to  compute  the 
input  and  antecedent  operations  for  interval  type-2  FLSs,  one  that  is  based  on  a  general  inference 
formula  for  them.  We  introduced  the  concept  of  upper  and  lower  membership  functions  (MFs)  and 
illustrate  our  efficient  inference  method  for  the  case  of  Gaussian  primary  MFs.  We  also  proposed 
a  method  for  designing  an  interval  type-2  FLS  in  which  we  tune  its  parameters. 

3.3  Applications  of  Interval  Type-2  Fuzzy  Logic  Systems 

We  have  developed  theory  and  design  methods  for  the  most  useful  kind  of  type-2  fuzzy  logic 
system  (FLSs),  interval  type-2  FLSs  [21],  and  have  applied  them  to  a  number  of  very  interesting 
applications,  such  as 

1.  Fading  channel  equalization  [22]  and  co-channel  interference  elimination  [23].  The  channel 
states  in  a  fading  channel  or  channel  with  co-channel  interferences  are  uncertain,  and  we 
validated  that  an  interval  type-2  fuzzy  set,  Gaussian  primary  membership  function  with 
uncertain  mean,  can  be  used  to  represent  such  uncertainties. 

2.  Network  video  traffic  modeling  and  classification  [19],  MPEG  variable  bit  rate  (VBR)  traffic 
are  very  bursty.  We  validated  that  the  I,  P,  and  B  frame  sizes  are  log-normal  with  fixed  mean 
and  uncertain  variance,  so  an  interval  type-2  fuzzy  sets  can  be  used  to  model  the  bursty  video 
traffic  and  an  interval  type-2  fuzzy  logic  system  with  such  type-2  fuzzy  set  are  demonstrated 
performing  much  better  than  a  Bayesian  classifier. 

3.  Connection  admission  control  for  ATM  network  [20],  Connection  admission  control  is  actually 
a  decision  making  problem.  Different  factors  such  as  incoming  real-time  video/audio  packet 
sizes,  non-real  time  packet  sizes,  the  buffer  sizes  are  uncertain.  We  applied  an  interval  type-2 
fuzzy  logic  to  handle  these  uncertainties,  and  achieved  very  good  performance. 
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4  Modeling  Node  Lifetime  with  Gaussian  Membership  Functions 


Wireless  sensor  nodes  are  severely  energy  constrained  due  to  their  compact  form.  To  increase 
the  lifetime  of  sensor  networks,  hardware  design  and  protocol  approaches  for  different  layers  must 
take  energy  efficiency  into  account.  However,  a  foundamental  question  -  “what  is  the  nature  of 
sensor  network  lifetime?”  has  not  been  answered  yet.  Since  the  lifetimes  of  individual  nodes  are 
not  constants  but  random  variables,  it  follows  that  the  network  lifetime  is  also  a  random  variable. 
Recent  research  by  Jain  and  Liang  [7]  discovered  that  in  a  wireless  sensor  network  where  the 
workloads  are  well-balanced,  a  single  node  lifetime  behaves  the  nature  of  normal  distribution  as 
demonstrated  in  Fig  5. 


Figure  5:  Single  Node  Lifetime  Distribution 

In  this  paper,  we  are  first  interested  in  setting  up  fine  membership  functions  (MFs)  for  single 
node  lifetime.  From  the  original  data  of  single  node  lifetime  shown  in  Table  1,  we  decompose  the 
whole  data  sets  into  seven  segments,  and  compute  the  mean  /x*  and  standard  deviation  c j*  of  node 
lifetime  for  each  segment,  i  =  1, 2,  ...7.  The  mean  ji  and  standard  deviation  a  for  of  the  entire  data 
set  is  also  computed.  We  are  also  interested  to  know  which  value  -  mean  fii  or  standard  deviation 
(Ji  varies  more.  We  first  normalize  the  mean  and  standard  deviation  of  each  segment  using 
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fii/fi  and  Gif  a.  Then  we  compute  the  standard  deviation  of  their  normalized  values  am  and  cr5^. 
Results  are  presented  in  Table  1. 

From  the  last  row  of  Table  1,  we  see  that  am  «  a^d  which  means  standard  deviation  G{ 
varies  much  more  than  the  mean  value  fi{.  Therefore  we  conclude  that  if  the  single  node  lifetime 
follows  normal  distribution,  it  is  most  appropriate  to  be  modeled  as  a  Gaussian  MF  with  uncertain 
standard  deviation.  One  example  of  type-2  Gaussian  MF  with  uncertain  standard  deviation  is 
shown  in  Fig  6.  This  result  also  justifies  the  use  of  the  Gaussian  MFs  to  model  network  lifetime  in 
section  5. 


Figure  6:  Type-2  Gaussian  MF  with  uncertain  standard  deviation.  The  thick  solid  lines  denote 
upper  MFs,  and  the  thick  dashed  lines  denote  lower  MFs.  The  shaded  regions  are  the  footprints  of 
uncertainty  for  interval  secondaries.  The  center  of  the  Gaussian  MFs  is  5,  and  the  variance  varies 
from  1.0  to  2.0 
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Table  1:  MEAN  AND  STD  VALUES  FOR  SEVEN  SEGMENTS  AND  THE  ENTIRE  NODE 
LIFETIME,  AND  THEIR  NORMALIZED  STD 


Node  Lifetime  Data 

Mean 

Std 

Segment  1 

1027.4 

30.182 

Segment  2 

1028.9 

29.819 

Segment  3 

1026.3 

30.798 

Segment  4 

1028.7 

30.917 

Segment  5 

1028 

29.944 

Segment  6 

1027 

29.975 

Segment  7 

1027.9 

30.306 

Entire  Data  Set 

1027.7 

30.292 

Normalized  STD 

0.00082783 

0.013105 

5  Network  Lifetime  Analysis  and  Estimation  Using  Interval  Type- 
2  Fuzzy  Logic  Systems 

We  are  now  ready  to  evaluate  the  network  lifetime  using  interval  type-2  fuzzy  logic  systems  (FLSs). 
We  apply  reliability  theory  to  design  fuzzy  rules.  In  the  following  section,  we  first  treat  the  basics 
of  reliability  theory  before  extracting  the  knowledges  for  fuzzy  rules.  An  interval  type-2  FLS  is 
then  set  up  to  perform  lifetime  analysis  and  estimation. 

5.1  Fuzzy  Rules  Design  Using  Reliability  Theory 

Rules  are  the  heart  of  a  FLS,  and  may  be  provided  by  experts  or  can  be  extracted  from  numerical 
data.  In  either  case,  the  rules  that  we  are  interested  in  can  be  expressed  as  a  collection  of  IF-THEN 
statements,  e.g., 

IF  the  noise  of  the  wireless  channel  is  high,  THEN  the  quality  of  the  received  signal  is  poor. 
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For  the  lifetime  issue  studied  in  this  paper,  reliability  theory  provides  a  feasible  method  to 
design  fuzzy  rules.  To  understand  this,  we  introduce  the  reliability  block  diagram  ( RBD ).  RBD  is 
a  graphical  representation  of  the  components  of  the  system,  and  provides  a  visual  representation 
of  the  way  components  are  reliability-wise  connected.  Thus  the  effect  of  the  success  or  failure  of  a 
component  on  the  system  performance  can  be  evaluated. 

Consider  a  system  with  two  components.  If  this  system  is  such  that  a  single  component  fail¬ 
ure  can  render  the  system  nonfunctional,  then  we  say  that  the  components  are  reliability- wise, 
connected  in  series.  If  the  system  fails  only  when  both  its  components  fail,  then  we  say  that  the 
components  are  reliability-wise  connected  in  parallel.  Note  that  the  physical  connection  between 
the  component  may  or  may  not  be  different  from  their  reliability-wise  connection.  The  RBD’s  for 
both  cases  are  given  in  Figs.7.  Any  complex  system  can  be  realized  in  the  form  of  a  combination 
of  blocks  connected  in  series  and  parallel. 


Component  1 


(a) 

-T  Component  2 


(b) 


Figure  7:  Reliability  block  diagrams  (RBD)  for  a  system  of  two  components.  (a)RBD  with  series 
connected  components.  (b)RBD  with  parallel  connected  components. 

In  our  analysis,  the  wireless  sensor  network  is  the  system  under  consideration  and  the  sensor 
nodes  are  the  components  of  the  system.  We  demonstrate  bellow  how  the  knowledge  is  extracted 
for  fuzzy  rules  design  referring  to  the  two  examples  in  Fig.  7. 


Example  Is  Set  up  Fuzzy  Rules  for  Parallel  System  in  Fig.  7(a) 
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In  the  parallel  system,  the  system  (network)  fails  only  when  both  components  (sensor  nodes) 
fail.  The  rules  can  be  set  up  as  one  example  shown  bellow: 

IF  the  remaining  battery  level  of  component  1  (sensor  node  1)  is  low  and  the  remaining  battery 
level  of  component  2  (sensor  node  2)  is  moderate ,  THEN  the  lifetime  of  the  system  (network)  is 

moderate. 

Example  2:  Set  up  Fuzzy  Rules  for  Series  system  in  Fig.  7(b) 

In  the  series  system,  the  system  (network)  fails  when  either  component  fails.  The  rules  can  be 
set  up  as  one  example  shown  bellow: 

IF  the  remaining  battery  level  of  component  1  (sensor  node  1)  is  low  and  the  remaining  battery 
level  of  component  2  (sensor  node  2)  is  moderate ,  THEN  the  lifetime  of  the  system  (network)  is 

low. 

Note  that  the  parallel  and  series  systems  are  the  two  basic  ways  to  model  two  sensor  nodes.  A 
wireless  sensor  network  consisting  of  multiple  sensor  nodes  can  be  first  represented  in  the  reliability 
block  diagram.  Fuzzy  rules  are  then  set  up  based  on  the  above  examples. 

5.2  Simulation  and  Discussion 

In  our  simulations,  interval  type-2  FLSs  are  constructed  separately  for  square-grid  netowrk  and 
hex-grid  network.  We  take  the  remaining  energy  level  of  each  individual  sensor  node  as  input  to 
the  interval  type-2  FLS  and  the  output  is  the  estimated  network  lifetime.  As  we  have  discussed 
in  section  4,  we  choose  type- 2  Gaussian  with  uncertain  standard  deviation  as  the  membership 
functions  for  both  antecedent  and  consequent.  The  linguistic  variables  to  represent  the  antecedent 
remaining  energy  level  are  divided  into  three  levels:  high,  moderate  and  low  and  the  consequent 
-  estimated  network  life  is  divided  into  five  levels:  very  high ,  high ,  moderate ,  low  and  very  low. 

Simulations  are  performed  for  both  squre-grid  and  hex-grid  sensor  networks.  In  both  cases,  36 
nodes  are  deployed  and  the  distance  between  neighboring  nodes  is  assumed  to  be  the  same.  All 
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sensor  nodes  are  initialized  with  maximum  energy  level  of  10  Joules.  We  assume  the  workload  of 
the  entire  network  is  well-balanced,  therefore  energy  dissipation  of  the  single  sensor  node  can  be 
figured  as  linear  most  of  the  alive  time.  Our  simulations  are  based  on  N  =  600  lifetime  data  from 
the  above  well-balanced  sensor  network.  The  first  300  data  are  for  training  and  the  remaining  300 
data  are  for  testing.  After  training,  the  rules  are  fixed  and  we  test  the  interval  type-2  by  evaluating 
the  root  mean  square  errors  (RMSE)  between  the  defuzzified  output  of  FLS  and  the  real  lifetime 
data.  We  compare  our  interval  type-2  FLSs  with  singleton  type-1  FLSs.  For  both  FLS  schemes, 
we  run  200  Monte-Carlo  realizations  and  for  each  realization,  each  FLS  is  tuned  using  a  simple 
steepest-descent  algorithm  for  six  epochs.  Simulation  results  are  averaged  over  all  200  Monte-Carlo 
realizations. 


5.2.1  Square  Grid 


As  defined  in  Section  2.4,  the  minimum  network  lifetime  is  the  time  to  failure  of  any  two  neighboring 
nodes.  We  know  that  the  failure  of  any  single  node  does  not  cause  network  failure.  The  failure 
of  any  node  coupled  with  the  failure  of  any  of  its  neighbors  causes  network  failure.  Using  this 
definition  we  build  the  RBD  for  the  square-grid  as  shown  in  Fig  8. 
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Figure  8:  RBD  of  a  single  node  in  a  square  grid.  Nodes  belonging  to  region-1  are  modeled  as 
block-1  and  nodes  belonging  to  region-2  axe  modeled  as  block-2.  The  network  RBD  consists  of 
W Nmin  -  l)2  block-l’s  in  series  with  2(s/Nmin  -  1)  block-2’s 
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Fig  8  shows  the  RBD  block  for  a  single  node  in  the  network.  A  node  can  be  modeled  in  two 
ways  depending  on  its  position  in  the  sensor  field.  This  distinction  based  on  its  position  is  made 
due  to  a  simple  observation  that  nodes  at  the  right  edge  of  the  sensor  field  (region-2)  do  not  have 
any  right  neighbor  (node  b)  as  opposed  to  nodes  in  region-1.  Also,  nodes  at  the  bottom  edge  of  the 
sensor  field  (region-2)  do  not  have  a  bottom  neighbor  (node  c)  as  opposed  to  the  nodes  in  region-1. 
Note  that  as  every  node  in  a  square-grid,  node  a  has  four  neighbors,  but  its  relationship  with  only 
two  neighbors  is  modeled  in  its  RBD  block.  This  is  because  the  relationship  with  the  other  two 
neighbors  will  be  modeled  when  their  RBD  blocks  are  constructed.  If  this  is  not  followed  then  the 
relationship  between  every  node-neighbor  pair  will  be  modeled  twice. 

In  this  square  grid  network  shown  in  Fig  8,  we  classify  three  antecedents  based  on  the  RBD  of 
block-1: 

•  The  remaining  battery  level  of  node  a. 

•  The  minimum  remaining  battery  level  of  node  b  and  c. 

•  The  remaining  battery  level  of  node  d. 

We  set  up  27  rules  for  this  FLS  because  every  antecedent  has  3  fuzzy  sub-sets  and  there  are  3 
antecedents. 

Two  antecedents  are  chosen  based  on  the  RBD  of  block-2  in  Fig  8  and  total  9  rules  are  con¬ 
structed  in  this  case. 

•  The  remaining  battery  level  of  node  x. 

•  The  remaining  battery  level  of  node  y. 

Let  Nmin  be  the  number  of  sensor  nodes  required  to  be  deployed  with  minimum  density.  The 
network  RBD  consists  of  —  l)2  block-l’s  and  2 (y/~Nmin  —  1)  block-2’s  in  series,  the  whole 

square  grid  network  can  actually  be  decomposed  into  multiple  blocks  serial  connected  together  and 
the  method  of  setting  up  rules  can  be  applied  as  well.  Simulation  results  of  RMSE  is  shown  in  Fig. 
9. 
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Figure  9:  Square-Grid:  The  RMSE(for  the  test  data)  for  two  FLS  approaches  averaged  over  200 
Monte-Carlo  realization 

5.2.2  Hex-Grid 


The  analysis  for  the  hex-grid  is  carried  out  on  the  same  lines  as  that  of  the  square-grid.  Fig  2  (b) 
shows  that  as  in  the  case  of  a  square  grid,  two  neighboring  node  failures  cause  network  failure.  The 
RBD  block  of  a  single  node  is  shown  in  Fig.  10. 


Ob 

Oa 

Oc  Od 


Figure  10:  RBD  block  for  a  single  node  in  the  Hex-grid:  The  network  RBD  consists  of  N/2  such 
blocks  in  series. 

Since  the  relation  between  a  node  and  all  of  its  neighbors  is  modeled  by  its  corresponding  RBD 
block,  the  RBD  block’s  for  the  neighbors  is  not  constructed  as  this  causes  the  relationship  between 
the  nodes  to  be  considered  twice.  In  this  case,  we  classify  two  antecedents  and  construct  9  rules. 
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•  The  remaining  battery  level  of  node  a. 


•  The  minimum  remaining  battery  level  of  node  b,  c  and  d. 

Since  Nmin/2  such  blocks  connected  in  series  represent  the  network,  the  whole  hex  grid  network 
can  be  decomposed  the  same  way  as  in  square  grid  networks.  Simulation  results  of  RMSE  is  shown 
in  Fig.  11. 


Figure  11:  Hex-Grid:  The  RMSE(for  the  test  data)  for  two  FLS  approaches  averaged  over  200 
Monte-Carlo  realization 

Observe  Fig.9  and  Fig.ll,  interval  type-2  FLSs  outperform  singleton  type-1  FLSs.  This  results 
shows  that  the  interval  type-2  FLSs  are  more  feasible  for  real-time  energy  estimation. 


6  Conclusion  and  future  works 

In  this  paper,  we  describe  a  new  method  based  on  fuzzy  logic  theory  to  analyze  and  estimate  the 
network  lifetimes  for  wireless  sensor  networks.  Our  approach  is  illuminated  by  the  discovery  that  a 
single  node  lifetime  behaves  the  nature  of  normal  distribution.  However, we  deem  that  if  the  single 
node  lifetime  follows  normal  distribution,  it  is  most  appropriate  to  be  modeled  as  a  Gaussian  MF 
with  uncertain  standard  deviation.  We  then  set  up  the  interval  type-2  FLSs  for  energy  estimation 
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and  evaluate  their  performance  using  real  lifetime  data.  Simulation  results  justified  the  feasibility  of 
applying  type-2  FLSs  into  wireless  sensor  network  lifetime  analysis.  Interval  type-2  FLSs  provides 
a  way  to  handel  knowledge  uncertainty.  We  believe  that  our  approaches  opens  up  a  new  vision  for 
research  on  sensor  network  lifetimes. 

Our  future  work  will  focus  on  lifetime  evaluation  under  the  circumstances  that  the  task  schedul¬ 
ing  is  variable  and  how  the  estimated  network  lifetime  could  be  used  to  accomodate  the  scheduling 
change. 


Acknowledgment 

This  work  was  supported  by  the  U.S.  Office  of  Naval  Research  (ONR)  Young  Investigator  Award 

under  Grant  N00014-03-1-0466,  “Energy  Efficient  Wireless  Sensor  Networs  for  Future  Combat 

System  Using  Fuzzy  Logic” . 

References 

[1]  C.Y.  Chong,  S.  P.  Kumar,  “Sensor  Networks:  Evolution,  Opportunities,  and  Challenges”  Proc. 
IEEE  ,  vol  91,  no.  8  ,  Aug  2003,  pp.  1247  -1256 

[2]  A.  Mainwaring,  J.  Polastre,  R.  Szewczyk,  D.  Culler,  J.  Anderson,  “Wireless  Sensor  Networks 
for  Habitat  Monitoring”  Proc.  WSNA  ’02  Atlanta,  Georgia,  Sep  28,  2002. 

[3]  V.  A.  Kottapalli,  A.  S.  Kiremidjian,  J.  P.  Lynch,  Ed  Carryer,  T.  W.  Kenny,  “Two-tired  wireless 
sensor  network  architecture  for  structural  health  monitoring”  Proc.  SPIE  San  Diego,  CA,  Mar 
2003. 

[4]  M.  Bhardwaj,  T.  Garnett,  A.  Chandrakasan,  “Upper  Bounds  on  the  Lifetime  of  Sensor  Net¬ 
works”  Proc.  IEEE  International  Conference  on  Communications ,  pp.785-790,  2001. 


22 


[5]  M.  Bhardwaj,  A.  P.  Chandrakasan,  “Bounding  the  Lifetime  of  Sensor  Networks  via  Optimal 
Role  Assignments”  Proc.  INFOCOM  2002  pp:1587  -  1596,  voL3 

[6]  D.  M.  Blough,  P.  Santi,  “Investigating  Upper  Bounds  on  Network  Lifetime  Extension  for 
Cell-Based  Energy  Conservation  Techniques  in  Stationary  Ad  Hod  Networks”  Proc.  MOBI- 
COM’2002  Atlanta,  Georgia,  Sep  2002 

[7]  E.  Jain  and  Q.  Liang,  “Sensor  placement  and  lifetime  of  wireless  sensor  networks:  theory  and 
performance  analysis,”  Sensor  Network  Operations ,  edited  by  S.  Phoha,  T.  F.  La  Porta,  and 
C.  Griffin,  IEEE  Press,  Piscataway,  NJ,  2004. 

[8]  N.  N.  Karnik  and  J.  M.  Mendel,  An  Introduction  to  Type-2  Fuzzy  Logic  Systems ,  October  1998, 
USC  Report,  http://sipi.usc.edu/~mendel/report. 

[9]  N.  N.  Karnik,  J.  M.  Mendel,  and  Q.  Liang,  “Type-2  fuzzy  logic  systems”,  IEEE  Trans.  Fuzzy 
Systems ,  vol.  7,  no.  6,  pp.  643-658,  Dec.  1999. 

[10]  N.  N.  Karnik  and  J.  M.  Mendel,  “Type-2  Fuzzy  Logic  Systems  :  Type-Reduction”,  presented 
at  the  1998  IEEE  SMC  Conference ,  San  Diego,  CA,  October. 

[11]  N.  N.  Karnik,  J.  M.  Mendel,  and  Q.  Liang,  “Centroid  of  a  type-2  fuzzy  set,”  Information 
Sciences,  vol.132,  pp. 195-220,  Feb,  2001, 

[12]  J.-S.  R.  Jang,  “ANFIS:  adaptive-network-based  fuzzy  inference  system,”  IEEE  Trans,  on 
Systems,  Man,  and  Cybernetics ,  vol.  23,  no.  3,  pp.  665-685,  May/ June  1993. 

[13]  K.  Kar,  S.  Banerjee,  “Node  Placement  for  Connected  Coverage  in  Sensor  Networks”  Extended 
Abstract.  Proc.  WiOpt  2003  Sophia-Antipolis,  France,  March  2003. 

[14]  S.  S.  Dhillon,  K.  Chakrabarty,  S.  S.  Iyengar,  “Sensor  Placement  for  Grid  Coverage  under 
Imprecise  Detections,”  FUSION,  2002. 

[15]  L.  A.  Zadeh,  “  The  concept  of  a  linguistic  variable  and  its  application  to  approximate  reasoning 
-  I,”  Information  Sciences ,  vol.  8,  pp.  199-249,  1975. 


23 


[16]  J.  M.  Mendel,  “Fuzzy  logic  systems  for  engineering:  a  tutorial,”  Proc.  of  the  IEEE,  vol.  83, 
no.  3,  pp.  345-377,  March  1995. 

[17]  D.  Dubois  and  H.  Prade,  Fuzzy  Sets  and  Systems:  Theory  and  Applications,  Academic  Press, 
New  York,  USA,  1980. 

[18]  E.  Hisdal,  “The  IF  THEN  ELSE  statement  and  interval-valued  fuzzy  sets  of  higher  type,” 
Int’l.  J.  Man-Machine  Studies,  vol.  15,  pp.  385-455,  1981. 

[19]  Q.  Liang  and  J.  M.  Mendel,  “MPEG  VBR  video  traffic  modeling  and  classification  using  fuzzy 
techniques,”  IEEE  Transactions  on  Fuzzy  Systems,  vol.  9,  no.  1,  pp.183-193,  Feb  2001. 

[20]  Q.  Liang,  N.  Karnik,  and  J.  M.  Mendel,  “  Connection  admission  control  in  ATM  network 
using  survey-based  type-2  fuzzy  logic  systems,”  IEEE  Transactions  on  Systems,  Man,  and 
Cybernetics,  Part  C,  vol.  30,  no.  3,  pp.  529-539,  August  2000. 

[21]  Q.  Liang  and  J.  M.  Mendel,  “Interval  type-2  fuzzy  logic  systems:  theory  and  design,”  IEEE 
Transactions  on  Fuzzy  Systems,  vol.  8,  no.  5,  pp.  535-550,  Oct  2000. 

[22]  Q.  Liang  and  J.  M.  Mendel,  “Equalization  of  time-varying  nonlinear  channels  using  type-2 
fuzzy  adaptive  filters,”  IEEE  Trans,  on  Fuzzy  Systems,  vol.  8,  no.  5,  pp.  551-563,  Oct  2000. 

[23]  Q.  Liang  and  J.  M.  Mendel,  “  Overcoming  time-varying  co-channel  interference  using  type-2 
fuzzy  adaptive  filters”,  IEEE  Transactions  on  Circuits  and  Systems,  II,  vol.  47,  no.  12,  pp. 
1419-1428,  Dec  2000. 

[24]  L.-X.  Wang  and  J .  M.  Mendel,  “  Fuzzy  basis  functions,  universal  approximation,  and  orthog¬ 
onal  least  squares  learning,  ”  IEEE  Trans,  on  Neural  Networks,  vol.  3,  pp.  807-814,  Sept. 
1992. 

[25]  D.  Kececioglu,  “Reliability  Engineering  Handbook”  Volume  1  and  2,  Prentice  Hall,  Ney 
Jersey  1991. 


24 


[26]  L.  M.  Leemis,  “Reliability:  Probabilistic  Models  and  Statistical  Methods,”  Prentice-  Hall, 

1995 

[27]  Life  Data  Analysis  Reference.  [Online]  Available:  http://www.weibull.com/lifedatawebcontents.htm 

[28]  A.  Papoulis,  S.  U.  Pillai  “Probability,  Random  Varibales  and  Stochastic  Processes,  ”  4th  ed, 
McGraw-Hill,  New  York  2002. 

[29]  B.  Healy,  “The  Use  of  Wireless  Sensor  Networks  for  Mapping  Environmental 

Conditions  in  Buildings”  ASHRAE  Seminar,  July  2  2003  Available  Online: 

http://www.nist.gov/tc75/ASHRAESummer2003SeminarHealy.pdf 


25 


Fuzzy  Optimization  for  Distributed 
Sensor  Deployment 

Haining  Shu,  Qilian  Liang 
Department  of  Electrical  Engineering 
University  of  Texas  at  Arlington 
Arlington,  TX  76019-0016  USA 
E-mail:  shu@wcn.uta.edu,  liang@uta.edu 


Abstract —  The  effectiveness  of  distributed  wireless  sensor  net¬ 
works  highly  depends  on  the  deployment  of  sensors.  Given 
a  finite  number  of  sensors,  optimizing  the  sensor  deployment 
will  enhance  the  field  coverage  of  a  wireless  sensor  network. 
Network  lifetime  and  quality  of  communication  in  terms  of 
outage  probability  as  a  result,  will  be  greatly  ameliorated  as  the 
topology  approaches  uniformity  fast.  In  this  paper,  we  propose 
a  fuzzy  optimization  algorithm  (FOA)  to  efficiently  adjust  the 
sensor  placement  after  an  initial  random  deployment.  We  apply 
fuzzy  logic  theory  to  handle  the  uncertainty  in  sensor  deployment 
problem.  Simulation  results  show  that  our  approach  achieves  fast 
and  stable  deployment  and  maximize  the  field  coverage.  Outage 
probabiliy,  as  a  measure  of  communication  quality  gets  effectively 
decreased. 

I.  Introduction 

Sensor  networks  consist  of  certain  amount  of  small  and 
energy  constrained  nodes.  Sensor  nodes  are  deployed  in  sup¬ 
port  of  various  missions  including  environment  monitoring, 
battlefield  surveillance,  and  emergency  search  and  rescue. 

A  number  of  applications  require  the  placement  of  sensors 
at  desired  locations.  Such  placement-friendly  sensor  networks 
are  widely  used  for  infrastructure  security  [1],  where  critical 
buildings  and  facilities  such  as  airports  and  power  plants  are 
monitored  by  a  network  of  sensors  placed  adequately. 

Meanwhile,  other  applications  employ  sensor  nodes  with 
certain  mobility  like  mobile  robots.  Mobile  sensors  are  prac¬ 
tically  desirable  for  they  have  the  capability  to  move  around 
and  re-adjust  their  positions  for  high  quality  communication 
and  better  surveillance  [2]. 

Our  primary  interest  lies  in  the  wireless  sensor  network 
comprised  of  mobile  sensors.  Our  goal  is  to  optimize  the 
sensor  deployment  such  that  the  maximum  field  coverage  and 
high  quality  communication  could  be  achieved. 

Some  prior  research  proposed  a  strategy  based  on  virtual 
forces  in  sensor  deployment  and  target  localization  [3].  A 
distributed  self-spreading  algorithm  was  introduced  in  [4]  to 
improve  the  network  coverage.  Poduri  et  al  [5]  proposed  an 
algorithm  with  the  constraint  that  each  of  the  nodes  has  at 
lease  K  neighbors.  These  algorithms  have  made  lots  efforts 
to  formulate  the  virtual  forces,  however  none  of  which  can 
well  handle  the  random  move  and  unpredictable  oscillation  in 
deployment. 


Our  algorithm  is  very  different  from  all  previous  works. 
Instead  of  attempting  to  formulate  the  virtual  forces,  we 
propose  to  use  fuzzy  logic  system  in  control  of  the  sensor 
movement.  To  save  battery  life,  we  refer  our  appoach  to 
open-loop  control  systems,  in  which  the  control  action  is 
independent  of  the  physical  system  output,  e.g.  feedback 
control.  By  applying  this  fuzzy  optimization  mechanism  to 
each  individual  mobile  sensor,  uncertain  exhaustive  move 
and  oscillation  is  efficiently  avoided  and  fast  deployment  is 
achieved.  The  entire  network  as  a  result,  survives  for  longer 
lifetime  and  the  quality  of  communication  in  terms  of  outage 
probability  is  greatly  ameliorated  as  the  topology  approaches 
uniformity.  A  concept  of  coherence  time  is  introduced  for  the 
purpose  of  synchronization  among  sensors. 

This  paper  is  organized  as  follows.  In  section  II,  we  briefly 
review  the  basic  concept  of  fuzzy  logic  systems.  Section  III 
details  the  Fuzzy  Optimization  Algorithm  (FOA)  designed  for 
sensor  deployment.  Simulation  and  the  key  results  of  this 
work  are  presented  in  Section  IV.  Section  V  concludes  with 
a  summary.  Fenton-Wilkinson  method  to  tackle  the  outage 
problem  is  expatiated  in  appendix. 

II.  Overview  of  Fuzzy  Logic  Systems 

Figure  1  shows  the  structure  of  a  fuzzy  logic  system 
(FLS)  [6].  When  an  input  is  applied  to  a  FLS,  the  inference 
engine  computes  the  output  set  corresponding  to  each  rule. 
The  defuzzifer  then  computes  a  crisp  output  from  these  rule 
output  sets.  Consider  a  p-input  1 -output  FLS,  using  singleton 
fuzzification,  center-of-sets  defuzzification  [7]  and  “IF-THEN” 
rules  of  the  form  [8] 

Rl  :  IF  xi  is  F'j  and  x2  is  ¥l2  and  •  •  •  and  xp  is  F*  THEN  y 
is  G'. 

Assuming  singleton  fuzzification,  when  an  input  x'  = 
{z'i, . . . ,  xfp}  is  applied,  the  degree  of  firing  corresponding  to 
the  It h  rule  is  computed  as 

Mf'  (xi)  * Mf‘04)  *  •  •  •  * MF‘ (£p)  =  (1) 

where  *  and  T  both  indicate  the  chosen  f-norm.  There  are 
many  kinds  of  defuzzifiers.  In  this  paper,  we  focus,  for  illustra¬ 
tive  purposes,  on  the  center-of-sets  defuzzifier  [7].  It  computes 
a  crisp  output  for  the  FLS  by  first  computing  the  centroid,  cGi, 


of  every  consequent  set  Gl,  and,  then  computing  a  weighted 
average  of  these  centroids.  The  weight  corresponding  to  the 
Ith  rule  consequent  centroid  is  the  degree  of  firing  associated 
with  the  Zth  rule,  7?  j/xfi  (»{),  so  that 


Vcos  (x')  = 


1  (^i) 


where  M  is  the  number  of  rules  in  the  FLS. 


(2) 


FUZZY  LOGIC  SYSTEM 


Fig.  1 .  The  structure  of  a  fuzzy  logic  system. 


III.  Extracting  the  Knowledge  for  Fuzzy 
Optimization  Algorithm 

A.  Assumptions  and  Notations 

•  Sensor  field  is  denoted  by  a  two-dimensional  grid. 

•  Coverage  discussed  in  this  paper  is  grid  coverage. 

•  A  grid  point  is  covered  when  at  least  one  sensor  covers 
this  point. 

•  A  sensor  can  detect  or  sense  any  event  within  its  sensing 
range,  denoted  by  Rs.  Coverage  is  determined  based  on 
Rs. 

•  Two  sensors  within  their  communication  range,  denoted 
by  Rc  can  communicate  with  each  other. 

•  Neighbors  of  a  sensor  are  nodes  within  its  communication 
range. 

•  Detection  and  communication  is  modeled  as  a  circle  on 
the  two-dimensional  grid. 

B .  Fuzzy  Optimization  Algorithm 

The  Fuzzy  Optimization  Algorithm  (FOA)  is  illuminated 
by  the  powerful  capability  of  fuzzy  logic  system  to  handle 
uncertainty  and  ambiguity.  Fuzzy  logic  system  is  well  known 
as  model  free.  Their  membership  functions  are  not  based 
on  statistical  distributions.  In  this  paper,  we  apply  fuzzy 
logic  system  to  re-position  the  sensor  nodes  and  optimize  the 
network  deployment. 

Our  algorithm  starts  with  random  deployment.  Initially,  a 
given  number  of  sensors  are  randomly  deployed  in  a  square 
sensor  field.  We  have  made  the  following  assumptions: 

•  All  sensor  nodes  are  peer  to  peer. 

•  Sensor  nodes  have  certain  mobility  and  capabilities  of 
computing,  detection  and  communication. 

•  Sensor  node  knows  its  location  information. 

•  Sensors  are  synchronized  by  coherence  time.  One-time 
move  is  made  within  each  coherence  period. 

Two  critical  procedures  are  considered  in  our  algorithm: 


•  Determine  the  next-step  move  distance  for  each  sensor. 

•  Determine  the  next-step  move  direction  for  each  sensor. 

The  next-step  move  distance  is  hard  to  determine.  Too  small 

or  big  move  distance  each  step  consumes  the  network  more 
time  and  energy  to  get  stable  deployment.  Excessive  oscillation 
is  unavoidable  in  previous  work  with  no  fuzzy  control.  In  this 
paper,  we  design  a  fuzzy  logic  system  to  determine  the  next- 
step  move  distance  for  each  sensor. 

We  collect  the  knowledge  for  deployment  problem  based 
on  the  following  two  antecedents: 

Antecedent  1 .  Number  of  neighbors  of  each  sensor. 

Antecedent  2.  Average  Euclidean  distance  between  sensor 
node  and  its  neighbors 

The  linguistic  variables  to  represent  the  number  of  neighbors 
for  each  sensor  are  divided  into  three  levels:  high ,  moderate 
and  low;  and  those  to  represent  the  average  Euclidean  distance 
between  sensor  node  and  its  neighbors  are  divided  into  three 
levels:  far ,  moderate  and  near .  The  consequent  -  the  shift 
distance  normalized  by  sensing  range  is  divided  into  three 
levels:  far ,  moderate  and  near.  Table  1  summaries  the  rules 
and  consequents. 


TABLE  I 

Fuzzy  Rules  and  Consequent 


Antecedentl 

Antecedent2 

Consequent 

Low 

Near 

Moderate 

Low 

Moderate 

Near 

Low 

Far 

Near 

Moderate 

Near 

Far 

Moderate 

Moderate 

Moderate 

Moderate 

Far 

Near 

High 

Near 

Far 

High 

Moderate 

Moderate 

High 

Far 

Moderate 

One  example  of  rules  is  as  follows: 

IF  the  number  of  neighbors  of  sensor  i  is  high  and  average 
Euclidean  distance  between  sensor  i  and  its  neighbors  is 
moderate ,  THEN  the  normalized  scalar  shift  distance  of 
sensor  i  will  be  high. 


We  set  up  9  rules  for  this  FLS  because  every  antecedent 
has  3  fuzzy  sub-sets  and  there  are  2  antecedents.  Trapezoidal 
membership  functions  (MFs)  are  used  to  represent  high,  low , 
far  and  near  and  triangle  MFs  to  represent  moderate.  We  show 
these  membership  functions  in  Figure  2. 

Applying  center-of-sets  defuzzification  [7],  for  every  input 
(zi,  £2),  the  output  is  computed  using 


V(xi,  x2) 


SLi  (>2) 

ELlMF'(zl)/*F'(Z2) 


(3) 


Repeating  these  calculations  for  Vxj  €  [0,1], we  obtain  a 
control  surface  y(x  1,  xf)  as  shown  in  Figure  3. 

The  concept  of  control  surface,  or  decison  surface,  is  central 
in  fuzzy  logic  systems.  It  describes  the  dynamics  of  the 


Fig.  2.  Antecedent  Membership  Function 


Fig.  3.  Control  Surface 


controller  and  is  generally  a  time-varying  nonlinear  surface. 
From  figure  3,  we  can  see  that  although  the  number  of 
neighbors  for  a  certain  sensor  is  high,  the  move  distance  can  be 
nearer  than  some  sensor  with  fewer  ’’crowded”  neighbors,  i.e. 
very  close  average  Euclidean  distance  between  the  sensor  and 
its  neighbors.  With  the  assist  of  control  surface,  the  next-step 
move  distance  can  be  carefully  determined. 

Comparing  to  move  distance,  the  next-step  move  direction 
is  much  easier  to  decide.  Coulomb’s  law  in  physics  becomes 
a  useful  tool  to  tackle  the  problem.  Assume  sensor  i  has  2 
neighbors  as  shown  in  Figure  4. 

The  coordinate  of  sensor  i  is  denoted  as  Ci  =  (Xi,Yi)- 
The  next-step  move  direction  of  sensor  i  could  be  repre¬ 
sented  as  follows: 


V  =  >  — =r - r; - 

U\Ci-C# 

(4) 

tan(a)  =  (  } 

X(v) 

(5) 

After  getting  distance  and  direction  (angles)  ,  sensor  i 
clearly  knows  his  next-step  move  position.  In  order  to  prolong 
the  battery  life  of  each  individual  sensor,  we  introduce  a 
coherence  time  as  the  duty  cycle  of  executing  the  FOA 


Fig.  4.  Next  Step  Move  Direction 


algorithm.  Sensors  are  put  into  idle  or  sleep  mode  if  within 
the  coherence  time,  the  information  of  neighbors  remains 
unchanged. 

IV.  Simulation  and  Discussion 

We  investigate  various  number  of  sensors  deployed  in  a 
field  of  10  x  10  square  kilometers  area.  We  assume  each  sensor 
is  equipped  with  an  omni  antenna  to  carry  out  the  task  of 
detection  and  communication.  Evaluation  of  our  FOA  algo¬ 
rithm  follows  three  criteria:  field  coverage,  outage  probability 
and  convergence.  Results  are  averaged  over  200  Monte  Carlo 
simulations. 

Figure  5  shows  at  Rs= lkm  and  Rc= 2km,  the  coverage  of  the 
initial  random  deployment  and  the  one  after  FOA  algorithm. 
The  FOA  algorithm  could  improve  the  network  coverage 
by  20%  -  30%  in  average.  As  60  sensors  are  deployed, 
the  coverage  approaches  to  100%  after  FOA  algorithm  is 
implemented. 


Fig.  5.  Coverage  vs.  Number  of  Nodes  (Rc=2,Rs=l) 
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Fig.  6.  Coverage  vs.  Number  of  Nodes  (Rc=4,Rs=2) 


Figure  6  gives  the  results  when  Rs=2 km  and  Rc= 4km, 
the  coverage  comparison  between  random  deployment  and 
fuzzy  optimization  algorithm.  In  the  case  when  20  sensors 
are  deployed,  initialy  the  coverage  after  random  deployment 
is  around  85%.  After  FOA  algorithm  is  excuted,  the  coverage 
reaches  98%.  The  coverage  is  dramatically  improved  in  the 
low  density  network.  The  above  two  figures  indicate  that 
instead  of  deploying  large  amount  of  sensors,  the  desired  field 
coverage  could  also  be  achieved  with  fewer  sensors. 

In  celluar  radio  systems  the  radio  link  performance  is 
usually  limited  by  interference  rather  than  noise,  therefore, 
the  probability  of  outage  due  to  co-channel  interference  is  of 
primary  concern.  Measurements  [9]  have  shown  that  at  any 
value  of  dij  (the  Euclidean  distance  between  sensor  i  and 
sensory*),  the  path  loss  PL(dij)  is  random  and  distributed  log- 
normally  (normal  in  dB)  about  the  mean  distance  dependent 
value.  That  is: 


PL(du)[dB]  =  PL(dij)+Xa  =  PL(d0)+10nlog(^)+Xa 

do 

(6) 

and 

Pr(diij)[dBm]  =  Pt[dBm\  -  PL{di^)[dB]  (7) 

where  Xa  is  a  zero-mean  Gaussian  distribution  random 
variable  (in  dB)  with  standard  deviation  a  (also  in  dB). 

The  log-normal  distribution  describes  the  random  shadow¬ 
ing  effects  on  the  propagation  path  which  implies  that  mea¬ 
sured  signal  levels  at  certain  distance  have  a  Gaussian  (normal) 
distribution  about  the  distance- dependent  mean  and  standard 
deviation  cr.  Since  PL(dij)  follows  normal  distribution,  so  is 
Pr{di,j),  and  the  Q  function  may  be  used  to  determine  the 
probability  that  the  received  signal  level  will  exceed  (or  fall 
below)  a  particular  level. 

The  probability  that  the  received  signal  level  will  exceed  a 
certain  value  7  can  be  calculated  from  the  cumulative  density 
function  as 


Pr[Pr(dij)  >7 }  =  Q  (7  (8) 

For  sonsor  i  with  N  neighbors,  if  sensor  i  acts  as  the 
destination  node  during  one  communication,  the  signal  to 
interference  ratio  (SIR)  is  represented  as: 


SIR(i)  = 


Prjdjj) 

Etl  Pr(di,k) 


Mi 


(9) 


The  denominator  denoting  the  effect  of  co-channel  interfer¬ 
ence  is  a  sum  of  N  —  1  log-normal  signals.  Evaluating  the 
outage  probability  requires  the  probability  distribution  of  the 
interference  power.  There  is  no  known  exact  expression  for 
the  probability  distribution  for  the  sum  of  log-normal  random 
variables,  but  various  authors  have  derived  several  approaches 
which  approximate  the  sum  of  log-normal  random  variables 
by  another  log-normal  random  variable. 

In  this  paper,  we  introduce  Fenton- Wilkinson  method  [10]. 
The  co-channel  interference  can  now  be  approximated  by  one 
log-normal  random  variable.  SIR(in  dB)  as  a  result  follows 
log-normal  distribution  as  well.  We  expatiate  the  Fenton- 
Wilkinson  method  in  the  appendix.  Results  of  outage  prob¬ 
ability  are  presented  in  figure  7. 


Fig.  7.  Outage  Probability  vs.  Number  of  Sensors  (Rc:=4,Rs=2) 

From  Figure  7,  we  can  see  that  our  FOA  algorithm  suc¬ 
cessfully  reduced  the  outage  probability  by  nearly  30%  which 
implies  a  higher  probability  that  the  received  signal  level  will 
exceed  the  SIR  threshold.  The  quality  of  communications  then 
can  be  greatly  ameliorated. 

The  performance  of  FOA  algorithm  can  also  be  evaluated 
in  terms  of  convergence  speed.  We  demonstrate  the  coverage 
at  each  iteration  e.g.  coverage  at  the  ith  iteration  is  Cov(i) 
for  different  number  of  sensors  and  the  mean  square  errors 
(MSE)  between  adjacent  iterations  shown  as  below. 


MSE(i)  = 


[Cov{i)  —  Cov(i  —  l)]2  . 
Cov(i  —  l)2 


1, ...,  iteration 


(10) 


Fig.  8.  Coverage  vs.Interation  (Rc=2,Rs=l) 


Fig.  9.  Coverage  vs.Interation  (Rc=4,Rs=2) 


Figure  8  and  9  show  different  number  of  sensors  have 
approximate  same  convergence  speed.  The  coverage  increases 
vastly  within  the  first  2  to  3  iterations.  After  that  coverage 
approaches  stable  for  all  different  number  of  sensors.  Figure 
9  shows  at  Rs= lkm  and  Rc= 2km  the  coverage  mean  square 
errors  between  the  adjacent  iterations.  The  results  also  validate 
that  the  convergence  of  FOA  algorithm  is  independent  of  the 
number  of  sensors  to  be  deployed. 

Comparing  to  other  algorithm  such  as  distributed  self 
spreading  algorithm  [4]  which  takes  over  20  termination  times 
and  10  oscillations  and  status  limit  to  achieve  the  similar 
coverage  improvement,  our  FOA  algorithm  is  much  simpler 
to  implement  and  outperforms  in  a  fast  and  guaranteed  con¬ 
vergence. 

V.  Conclusions 

In  this  paper,  we  proposed  a  new  sensor  deployment  strat¬ 
egy  -  Fuzzy  Optimization  Algorithm  based  on  fuzzy  logic 
system.  Our  approach  has  a  great  advantage  to  deal  with  the 
randomness  in  sensor  deployment  which  is  particularly  useful 
when  emergency  rescue  or  redeployment  over  hostile  situation 


Fig.  10.  MSE  vs.  Coverage  (Rc=2,Rs-l) 


is  needed.  We  believe  that  in  an  energy  constraint  wireless 
sensor  network,  fast  and  efficient  deployment  strategy  is  a 
necessity  to  save  battery  power  and  extend  network  lifetime. 
Our  FOA  algorithm  is  capable  to  model  all  random  deploy¬ 
ment  with  a  fuzzy  logic  system.  The  network  coverage  as  a 
result  gets  greatly  improved  and  quality  of  communication  in 
term  of  outage  probablity  is  ameliorated.  Moreover,  the  FOA 
algorithm  brings  the  whole  network  to  a  stable  and  optimal 
deployment  very  soon  which  will  significantly  reduce  the 
energy  consumption.  Our  future  work  will  focus  on  modeling 
the  random  deployment  with  some  existing  pattern  so  that  the 
energy  consumption  can  be  further  studied  in  the  deployment 
problem. 
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Appendix 

A.  Multiple  Log-Normal  Interferers 

Consider  the  sum  of  TV/  log-normal  random  variables 


where  Z  —  eZ(dBm)-  For  example,  suppose  that  £lk>  k  — 
1, A/  have  mean  p^k,  k  =  1, Ni  and  identical  variances 
cr? .  Identical  variances  are  often  assumed  because  the  standard 
deviation  of  log-normal  shadowing  is  largely  independent  of 
the  radio  path  length.  Equating  the  means  on  both  sides  of 
(15) 


Ni  JVj 

I  =  ^  Ofe  =  ^2  10nM«Bm)/10  (H) 

k= 1  1 

where  the  flk(dBm)  are  Gaussian  random  variables  with 
mean  pnHdBrn)  and  variance  cr^fc,  and  the  Q*  = 
are  the  log-normal  random  variables.  Unfortunately,  there  is 
no  known  closed  form  expression  for  the  probability  density 
function  (pdf)  of  the  sum  of  multiple  (TV/  >  2)  log-normal 
random  variables.  However,  there  is  a  general  consensus  that 
the  sum  of  independent  log-normal  random  variables  can  be 
approximated  by  another  log-normal  random  variable  with 
appropriately  chosen  parameters.  That  is, 


nz 

I  ~  ^2 10afc^>/10  ~  ioz(dBm)/io  =  f  ^2) 

k—1 

where  Z^Bm)  is  a  Gaussian  random  variable  with  mean 
Hz  {dBm)  and  variance  o\.  The  problem  is  to  determine 
Vz{dBm)  and  variance  cr|  in  terms  of  the  pnHdBrn)  and 
variance  a ^fc,  k  —  1 ,  ...,7V/.  Several  methods  have  been 
suggested  in  the  literature  to  solove  this  problem  including 
those  by  Fenton,  Schwartz  and  Yen,  and  Farley.  Each  of  these 
methods  provides  varying  degrees  of  accuracy  over  specified 
ranges  of  the  shadow  standard  deviation  oq,  the  sum  /,  and 
the  number  of  interferes  TV/. 

B.  Fenton-Wilkinson  Method 

With  the  Fenton-Wilkinson  method,  the  mean  pz(dBm )  and 
variance  cr|  of  Z^dBm)  are  obtained  by  matching  the  first 
two  moments  of  the  sum  I  with  the  first  two  moments  of 
the  approximation  /.  To  derive  the  appropriate  moments,  it  is 
convenient  to  use  natural  logarithms.  We  write 


_  eeDfc(dBm)  _  eflk 

where  e  =  (Znl0)/10  =  0.23026  and  =  eClk{dBm) .  Note 
that  Mn k  =  €Mf2fc(dBm)  and  cr?  =  e2cr^fc .  The  nth  moment  of 
the  log-normal  random  variabfe  Q  k  can  be  obtained  from  the 
moment  generating  function  of  the  Gaussian  random  variables 
Clk  as 


=  E[entlk]  =  e1^  +(1/2)n2*k  (14) 

To  find  the  appropriate  moments  for  the  log-normal  approx¬ 
imation  we  can  use  (14)  and  equate  the  first  two  moments  on 
both  sides  of  the  equation 


Ni 


/=e< 


k=l 


eA  =  I 


Nj 

fir  =  E[I]  =  ^  E[e^]  =  E[e *]  =  E[I]  =  (16) 

k—1 

gives  the  result 
/Nx 

\fc= i 

Likewise,  we  can  equate  the  variances  on  both  sides  of 
(15)  under  the  assumption  that  the  f k  =  1,  are 

independent 

0*I=E[P]-fi  =  E[P]=0*t  (lg) 

giving  the  result 


s(i/2)^  _  m2+(i/2)4 


(17) 


{YA2^k  j  e<T"(e<T"  -  1)  =  e2M2e4(eCTl  -  1)  (19) 

By  squaring  each  side  of  (17)  and  dividing  each  side  of 
resulting  equation  by  the  respective  side  of  (19)  We  can  solve 
for  crlatz  in  terms  of  the  known  values  of  fc  =  1, ...,  iVj 
and  <t?.  Afterwards,  p%  can  be  obtained  from  (17).  This 
procedure  yields  the  following  solution: 


Fz  ~ 


(  Ni 


(20) 


—  In 


TZuYY  .A 
(E^e^)2  J 


(21) 


Finally,  Mz(dBm)  =  and  =  £_2<T|- 

The  accuracy  of  this  log-normal  approximation  can  be 
measured  in  terms  of  how  accurately  the  first  two  moments  of 
I{dB)  =  lOZo^io/  are  estimated,  and  how  well  the  cumulative 
distribution  function  (cdf)  of  I(dB)  is  described  by  a  Gaussian 
cdf.  In  problems  relating  to  the  co-channel  interference  outage 
in  celluar  radio  systems,  we  are  usually  interested  in  the  tails 
of  both  the  complementary  distribution  function  (cdfc)  Ff  = 
P{I  >  x)  and  the  cdf  FT(x)  =  1  -  Ff  -  P(I  <  x).  In  this 
case,  we  are  interested  in  the  accuracy  of  the  approximation 


Fdx)  «  P(e*  >x)  =  Q  (j—^j  (22) 

for  large  and  small  values  of  x.  It  will  be  shown  later  that 
the  Fenton-Wilkinson  method  can  approximate  the  tails  of  the 
cdf  and  cdfc  functions  with  good  accuracy. 


(15) 
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Abstract 

The  effectiveness  of  distributed  wireless  sensor  networks  highly  depends  on  the  deployment  of 
sensors.  Given  a  finite  number  of  sensors,  optimizing  the  sensor  deployment  will  enhance  the  field 
coverage  of  a  wireless  sensor  network.  Network  lifetime  and  quality  of  service  in  terms  of  outage 
probability  as  a  result,  will  be  greatly  ameliorated  as  the  topology  approaches  uniformity  fast.  In 
this  paper,  we  apply  fuzzy  logic  systems  to  efficiently  adjust  the  sensor  placement  after  an  initial 
random  deployment.  To  evaluate  the  outage  probability  due  to  co-channel  interference,  Fenton- 
Wilkinson  method  is  used  to  approximate  the  sum  of  log-normal  random  variables.  Simulation 
results  show  that  our  approach  achieves  fast  and  stable  deployment  and  maximizes  the  field 
coverage.  Outage  probabiliy,  as  a  measure  of  communication  quality  gets  effectively  decreased. 


1  Introduction 

Sensor  networks  consist  of  certain  amount  of  small  and  energy  constrained  nodes.  Sensor  nodes  are 
deployed  in  support  of  various  missions  including  environment  monitoring,  battlefield  surveillance, 
and  emergency  search  and  rescue. 

A  number  of  applications  require  the  placement  of  sensors  at  desired  locations.  Such  placement- 
friendly  sensor  networks  are  widely  used  for  infrastructure  security  [1],  where  critical  buildings  and 
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facilities  such  as  airports  and  power  plants  are  monitored  by  a  network  of  sensors  placed  adequately. 

Meanwhile,  other  applications  employ  sensor  nodes  with  certain  mobility  like  mobile  robots. 
Mobile  sensors  are  practically  desirable  because  they  have  the  capability  to  move  around  and  re¬ 
adjust  their  positions  for  high  quality  communication  and  better  surveillance  [2]. 

Our  primary  interest  lies  in  the  wireless  sensor  network  comprised  of  mobile  sensors.  Our  goal 
is  to  optimize  the  sensor  deployment  such  that  fast  deployment,  maximum  field  coverage  and  high 
quality  communication  could  be  achieved. 

Some  prior  research  proposed  a  strategy  based  on  virtual  forces  in  sensor  deployment  and  target 
localization  [3] .  A  distributed  self-spreading  algorithm  was  introduced  in  [4]  to  improve  the  network 
coverage.  Poduri  et  al  [5]  proposed  an  algorithm  with  the  constraint  that  each  of  the  nodes  has  at 
lease  K  neighbors.  These  algorithms  have  made  lots  efforts  to  formulate  the  virtual  forces,  however 
none  of  which  can  well  handle  the  uncertainties  such  as  random  move  and  unpredictable  oscillation 
in  deployment.  Various  constraints  like  oscillation  limit,  stable  status  and  number  of  neighbors  [5] 
therefore  have  to  be  imposed  to  avoid  excessive  sensor  movement. 

In  this  paper,  we  apply  fuzzy  logic  to  handle  these  uncertainties  and  design  fuzzy  logic  systems 
(FLSs)  for  distributed  sensor  deployment.  Instead  of  attempting  to  formulate  the  virtual  forces,  we 
propose  to  use  fuzzy  logic  system  to  control  the  sensor  movement.  To  save  battery  life,  we  refer  our 
appoach  to  open-loop  control  systems,  in  which  the  control  action  is  independent  of  the  physical 
system  output,  e.g.  feedback  control.  Each  individual  mobile  sensor  uses  a  FLS  to  determine  its 
moving  distance  and  direction,  and  uncertain  exhaustive  move  and  oscillation  is  efficiently  avoided 
and  fast  deployment  is  achieved.  The  entire  network  as  a  result,  survives  for  longer  lifetime  and 
the  quality  of  communication  in  terms  of  outage  probability  is  greatly  ameliorated  as  the  topology 
approaches  uniformity.  A  concept  of  coherence  time  is  introduced  for  the  purpose  of  synchronization 
among  sensors. 

The  rest  of  this  paper  is  organized  as  follows.  In  section  2,  we  briefly  review  the  basic  concept  of 
fuzzy  logic  systems.  Section  3  details  the  fuzzy  logic  systems  (FLSs)  designed  for  sensor  deployment. 
Simulation  and  discussion  are  presented  in  Section  4.  Section  5  concludes  with  a  summary.  Fenton- 
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Wilkinson  method  to  tackle  the  outage  problem  is  expatiated  in  appendix. 


2  Overview  of  Fuzzy  Logic  Systems 


FUZZY  LOGIC  SYSTEM 


Figure  1:  The  structure  of  a  fuzzy  logic  system. 

Figure  1  shows  the  structure  of  a  fuzzy  logic  system  (FLS)  [6].  When  an  input  is  applied  to 
a  FLS,  the  inference  engine  computes  the  output  set  corresponding  to  each  rule.  The  defuzzifer 
then  computes  a  crisp  output  from  these  rule  output  sets.  Consider  a  p-input  1-output  FLS,  using 
singleton  fuzzification,  center-of-sets  defuzzification  [7]  and  “IF-THEN”  rules  of  the  form  [8] 

Rl  :  IF  a?i  is  F[  and  x2  is  F^  and  •  •  •  and  xp  is  F{,,  THEN  y  is  G*. 

Assuming  singleton  fuzzification,  when  an  input  x'  =  {x[,.  ..,x'p}  is  applied,  the  degree  of  firing 
corresponding  to  the  Zth  rule  is  computed  as 

fJ'F[  ( *  ■  •  •  *  Mf lp(xp)  =  ^2=iMf 'J  (®i)  (1) 

where  *  and  T  both  indicate  the  chosen  f-norm.  There  are  many  kinds  of  defuzzifiers.  In  this  paper, 
we  focus,  for  illustrative  purposes,  on  the  center-of-sets  defuzzifier  [7].  It  computes  a  crisp  output 
for  the  FLS  by  first  computing  the  centroid,  Cqi,  of  every  consequent  set  G^,  and,  then  computing  a 
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weighted  average  of  these  centroids.  The  weight  corresponding  to  the  Zth  rule  consequent  centroid 
is  the  degree  of  firing  associated  with  the  Zth  rule,  7^/^ (a:-),  so  that 


Ucosfe  )  — 


where  M  is  the  number  of  rules  in  the  FLS. 


fi(z') 


(2) 


3  FLSs  for  Distributed  Sensor  Deployment 

3.1  Assumptions  and  Notations 

•  Sensor  field  is  denoted  by  a  two-dimensional  grid.  Detection  and  communication  is  modeled 
as  a  circle  on  this  grid. 

•  Coverage  discussed  in  this  paper  is  grid  coverage.  A  grid  point  is  covered  when  at  least  one 
sensor  covers  this  point. 

•  A  sensor  can  detect  or  sense  any  event  within  its  sensing  range,  denoted  by  Rs.  Coverage  is 
determined  based  on  Rs. 

•  Two  sensors  within  their  communication  range,  denoted  by  Rc  can  communicate  with  each 
other. 

•  Neighbors  of  a  sensor  are  defined  as  nodes  within  its  communication  range. 

3.2  Fuzzy  Logic  System  Design 

Fuzzy  logic  system  is  well  known  to  be  able  to  handle  uncertainty  and  ambiguity.  Practically 
not  all  uncertainty  is  random.  Some  forms  of  uncertainty  are  non-random  and  hence  not  suited 
to  treatment  or  modeling  by  probability  theory.  Fuzzy  set  theory  is  a  marvelous  tool  for  model¬ 
ing  the  kind  of  uncertainty  associated  with  vagueness,  or  with  a  lack  of  information  regarding  a 
pariticular  element  of  the  problem  at  hand.  Upon  concerning  the  distributed  sensor  deployment, 
the  moving  distance  and  direction  of  each  sensor  are  distributed  and  full  of  uncertainty  which  can 
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barely  be  described  by  some  random  distribution.  Fuzzy  logic  system  is  well  known  as  model  free. 
Their  membership  functions  are  not  based  on  statistical  distributions.  Therefore  we  propose  to 
apply  fuzzy  logic  system  to  the  distributed  sensor  deployment  problem.  Each  sensor  makes  fully 
distributed  decision  on  its  movement  based  on  FLS. 

Our  algorithm  starts  with  random  deployment.  In  the  initial  condition,  a  given  number  of 
sensors  are  randomly  deployed  in  a  square  sensor  field.  We  have  made  the  following  assumptions: 

•  All  sensor  nodes  are  peer  to  peer. 

•  Sensor  nodes  have  certain  mobility  and  capabilities  of  computing,  detection  and  communica¬ 
tion. 

•  Sensor  node  knows  its  location  information. 

•  Sensors  are  synchronized  by  coherence  time.  One-time  move  is  made  within  each  coherence 
period. 

Every  sensor  in  the  network  needs  to  know  which  direction  he  is  heading  to  and  how  far  he  can 
reach  in  the  next  step.  Thus  two  critical  procedures  are  considered  in  our  algorithm: 

•  Determine  the  next-step  move  distance  for  each  sensor. 

•  Determine  the  next-step  move  direction  for  each  sensor. 

The  next-step  move  distance  is  hard  to  determine.  Too  small  or  big  move  distance  each  step 
consumes  the  network  more  time  and  energy  to  get  stable  deployment.  Excessive  move  and  oscil¬ 
lation  is  unavoidable  in  previous  work  with  no  fuzzy  control.  In  this  paper,  we  design  a  fuzzy  logic 
system  to  determine  the  next-step  move  distance  for  each  sensor. 

We  collect  the  knowledge  for  deployment  problem  based  on  the  following  two  antecedents: 

Antecedent  1.  Number  of  neighbors  of  each  sensor. 

Antecedent  2.  Average  Euclidean  distance  between  sensor  node  and  its  neighbors 
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The  linguistic  variables  to  represent  the  number  of  neighbors  for  each  sensor  are  divided  into 
three  levels:  high ,  moderate  and  low ;  and  those  to  represent  the  average  Euclidean  distance  between 
sensor  node  and  its  neighbors  are  divided  into  three  levels:  far ,  moderate  and  near.  The  consequent 
-  the  shift  distance  normalized  by  sensing  range  Rs  is  divided  into  three  levels:  /ar,  moderate  and 
near.  Table  1  summaries  the  rules  and  consequents. 


Table  1:  Fuzzy  Rules  and  Consequent 


Antecedent 1 

Antecedent2 

Consequent 

Low 

Near 

Moderate 

Low 

Moderate 

Near 

Low 

Far 

Near 

Moderate 

Near 

Far 

Moderate 

Moderate 

Moderate 

Moderate 

Far 

Near 

High 

Near 

Far 

High 

Moderate 

Moderate 

High 

Far 

Moderate 

An  ideal  sensor  deployment  should  have  uniform  distribution  for  better  coverage.  But  in  random 
deployment,  coverage  uniformity  is  hardly  to  achieve  initially.  In  sensor  network  composed  of  mobile 
sensors,  each  sensor  detects  the  number  and  location  of  its  neighbors  and  decides  if  its  neighborhood 
is  sparse  or  over  crowded.  If  it  is  over  crowded,  the  sensor  makes  decision  using  FLSs  to  shift  a 
certain  distance  away  from  its  current  location.  If  it  has  very  few  and  sparse  neighbors,  the  sensor 
might  stand  still  or  shift  a  little  but  staying  more  closer  to  its  current  location. 

One  example  of  rules  is  as  follows: 

IF  the  number  of  neighbors  of  sensor  i  is  high  and  average  Euclidean  distance  between  sensor  i 
and  its  neighbors  is  moderate ,  THEN  the  normalized  shift  distance  of  sensor  i  should  be  high. 
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We  set  up  9  rules  for  this  FLS  because  every  antecedent  has  3  fuzzy  sub-sets  and  there  are 
2  antecedents.  Trapezoidal  membership  functions  (MFs)  are  used  to  represent  high ,  low,  far  and 
near  and  triangle  MFs  to  represent  moderate.  Two  antecedents  are  normalized  to  the  range  [0, 10]. 
We  show  these  membership  functions  in  Figure  2. 


Figure  2:  Antecedent  Membership  Function 

Applying  center-of-sets  defuzzification  [7],  for  every  input  (xi,  xf),  the  output  is  computed 
using 


2/(xi,X2)  —  y>9  /  ,  /  s 

Lz=l  M F*(*l)MFla(®2) 


(3) 


Repeating  these  calculations  for  V x i  G  [0, 10], we  obtain  a  control  surface  y{x\,  xf)  as  shown  in 
Figure  3. 


The  concept  of  control  surface,  or  decison  surface,  is  central  in  fuzzy  logic  systems.  It  describes 
the  dynamics  of  the  controller  and  is  generally  a  time-varying  nonlinear  surface.  From  figure  3, 
we  can  see  that  although  the  number  of  neighbors  for  a  certain  sensor  is  high,  the  move  distance 
can  be  nearer  than  some  sensor  with  fewer  ’’crowded”  neighbors,  i.e.  very  close  average  Euclidean 
distance  between  the  sensor  and  its  neighbors.  With  the  assist  of  control  surface,  the  next-step 
move  distance  can  be  determined. 
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iverage  Euclidean  Distance 


Number  of  Neighbors 


Figure  3:  Control  Surface 

Comparing  to  move  distance,  the  next-step  move  direction  is  much  easier  to  decide.  Coulomb’s 
law  in  physics  becomes  a  useful  tool  to  tackle  the  problem.  For  instance,  assume  sensor  i  has  2 
neighbors  in  its  communication  range  as  shown  in  Figure  4. 


Figure  4:  Next  Step  Move  Direction 

The  coordinate  of  sensor  i  is  denoted  as  C<  =  (XiyYi). 

The  next-step  move  direction  of  sensor  i  could  be  represented  as  follows: 
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(5) 


tan(a )  =  — r^- 
X(v) 

After  getting  distance  and  direction  (angle  a)  ,  sensor  i  clearly  knows  his  next-step  move 
information.  In  order  to  prolong  the  battery  life  of  each  individual  sensor,  we  introduce  a  coherence 
time  as  the  duty  cycle  during  which  the  changes  of  two  antecedents  can  be  ignored.  Sensors  are 
put  into  idle  or  sleep  mode  if  within  the  coherence  time,  the  information  of  neighbors  remains 
unchanged. 

4  Simulation  and  Discussion 

We  investigate  various  number  of  sensors  deployed  in  a  field  of  10x10  square  kilometers  area. 
We  assume  each  sensor  is  equipped  with  an  omni  antenna  to  carry  out  the  task  of  detection  and 
communication.  Evaluation  of  our  scheme  follows  three  criteria:  field  coverage,  outage  probability 
and  convergence.  Results  are  averaged  over  200  Monte  Carlo  simulations. 

Figure  5  shows  at  1  kilometer  sensing  range  (i?s=lkm)  and  2  kilometer  communication  range 
(i?c=2km),  the  coverage  of  the  initial  random  deployment  and  the  one  after  using  FLSs.  The  FLSs 
scheme  could  improve  the  network  coverage  by  20%  -  30%  in  average.  As  60  sensors  are  deployed, 
the  coverage  approaches  to  100%  after  FLSs  scheme  is  implemented. 

Figure  6  gives  the  results  when  Rs= 2km  and  i?c=4km,  the  coverage  comparison  between  random 
deployment  and  fuzzy  optimization  algorithm.  In  the  case  when  20  sensors  are  deployed,  the 
coverage  after  random  deployment  is  initialy  around  85%.  After  FLSs  are  used,  the  coverage 
reaches  approximate  98%.  The  coverage  is  dramatically  improved  in  the  low  density  network.  On 
the  other  hand,  the  above  two  figures  indicate  that  instead  of  deploying  large  amount  of  sensors, 
the  desired  field  coverage  could  also  be  achieved  with  fewer  sensors. 


Figure  5:  Coverage  vs.  Number  of  Nodes  (Rc~2,Rs=l) 

In  celluar  radio  systems  the  radio  link  performance  is  usually  limited  by  interference  rather 
than  noise,  therefore,  the  probability  of  outage  due  to  co-channel  interference  is  of  primary  concern. 
Measurements  [9]  have  shown  that  at  any  value  of  (the  Euclidean  distance  between  sensor  i 
and  sensor  j),  the  path  loss  PL(di}j)  is  random  and  distributed  log-normally  (normal  in  dB)  about 
the  mean  distance  dependent  value.  That  is: 


PL(ditj)[dB]  =  PL(dij)  +  Xa  -  PL(d0)  +  10 nlog(^-)  +  Xa  (6) 

do 

and 

Pr(dij)[dBm }  =  Pt[dBm ]  -  PL(dij)[dB]  (7) 

where  Xa  is  a  zero-mean  Gaussian  distribution  random  variable  (in  dB)  with  standard  deviation 
a  (also  in  dB). 

The  log-normal  distribution  describes  the  random  shadowing  effects  on  the  propagation  path 
which  implies  that  measured  signal  levels  at  certain  distance  have  a  Gaussian  (normal)  distribu¬ 
tion  about  the  distance-dependent  mean  and  standard  deviation  o.  Since  PL^dij)  follows  normal 
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Figure  6:  Coverage  vs.  Number  of  Nodes  (Rc=4,Rs=2) 

distribution,  so  is  Pr(dij),  and  the  Q  function  may  be  used  to  determine  the  probability  that  the 
received  signal  level  will  exceed  (or  fall  below)  a  particular  level. 

The  probability  that  the  received  signal  level  will  exceed  a  certain  value  7  can  be  calculated 
from  the  cumulative  density  function  as 


Pr  [Pr  (dij )  >7 }  —  Q 


7  Pr  (djj ) 


(8) 


For  sonsor  i  with  N  neighbors,  if  sensor  i  acts  as  the  destination  node  during  one  communication, 
the  signal  to  interference  ratio  (SIR)  is  represented  as: 


(9) 


E"=l  Pr(di,k) 

The  denominator  denoting  the  effect  of  co-channel  interference  is  a  sum  of  N  -  1  log-normal 
signals.  Evaluating  the  outage  probability  requires  the  probability  distribution  of  the  interference 
power.  There  is  no  known  exact  expression  for  the  probability  distribution  for  the  sum  of  log¬ 
normal  random  variables,  but  various  authors  have  derived  several  approaches  which  approximate 
the  sum  of  log-normal  random  variables  by  another  log-normal  random  variable. 
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In  this  paper,  we  use  Fenton-Wilkinson  method  [10].  The  co-channel  interference  can  now 
be  approximated  by  one  log-normal  random  variable.  SIR(in  dB)  as  a  result  follows  log-normal 
distribution  as  well.  We  expatiate  the  Fenton-Wilkinson  method  in  the  appendix.  Results  of  outage 
probability  are  presented  in  figure  7. 


Figure  7:  Outage  Probability  vs.  Number  of  Sensors  (Rc=4,Rs~2) 

From  Figure  7,  we  can  see  that  the  FLSs  scheme  successfully  reduced  the  outage  probability 
by  nearly  30%  which  implies  a  higher  probability  that  the  received  signal  level  will  exceed  the  SIR 
threshold.  The  quality  of  communications  then  can  be  greatly  ameliorated. 

The  performance  of  the  FLSs  scheme  can  also  be  evaluated  in  terms  of  convergence  speed.  We 
demonstrate  the  coverage  at  each  iteration  e.g.  coverage  at  the  ith  iteration  is  Cov[i)  for  different 
number  of  sensors  and  the  mean  square  errors  (MSE)  between  adjacent  iterations  shown  as  below. 

TkjroTPf\  [Cov(i)  -  Cov(i  -  l)]2 

MSE(i)  = - - 1  =  1,  ■■■,  iteration  (10) 

Figure  8  and  9  show  different  number  of  sensors  have  approximate  same  convergence  speed. 
The  coverage  increases  vastly  within  the  first  2  to  3  iterations.  After  that  coverage  approaches 
stable  for  all  different  number  of  sensors.  Figure  9  shows  at  Rs= lkm  and  i?c=2km  the  coverage 
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Figure  8:  Coverage  vs.Interation  (Rc=2,Rs=l) 

mean  square  errors  between  the  adjacent  iterations.  The  results  also  validate  that  the  convergence 
of  the  FLSs  scheme  is  independent  of  the  number  of  sensors  to  be  deployed. 

Comparing  to  other  algorithm  such  as  distributed  self  spreading  algorithm  [4]  which  takes  over 
20  termination  times  and  10  oscillations  and  status  limit  to  achieve  the  similar  coverage  improve¬ 
ment,  the  FLSs  scheme  is  much  simpler  to  implement  and  outperforms  in  a  fast  and  guaranteed 
convergence. 

5  Conclusions 

In  this  paper,  we  proposed  a  new  sensor  deployment  strategy  based  on  fuzzy  logic  system.  Our 
approach  has  a  great  advantage  to  deal  with  the  uncertainty  in  distributed  sensor  deployment  which 
is  particularly  useful  when  emergency  rescue  or  redeployment  over  hostile  situation  is  needed.  We 
believe  that  in  an  energy  constraint  wireless  sensor  network,  fast  and  efficient  deployment  strategy 
is  a  necessity  to  save  battery  power  and  extend  network  lifetime.  Our  FLSs  scheme  is  capable 
to  model  all  distributed  sensor  deployment  with  a  fuzzy  logic  system.  The  network  coverage  and 
quality  of  communication  in  term  of  outage  probablity  are  greatly  improved  as  a  result.  Moreover, 
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Figure  9:  Coverage  vs.Interation  (Rc=4,Rs=2) 

the  FLSs  scheme  brings  the  whole  network  to  a  stable  and  optimal  deployment  very  soon  which  will 
significantly  reduce  the  energy  consumption.  Our  future  work  will  focus  on  modeling  the  random 
deployment  with  some  existing  pattern  so  that  the  energy  consumption  can  be  further  studied  in 
the  deployment  problem. 
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Appendix 

A  Multiple  Log-Normal  Interferers 

Consider  the  sum  of  Nj  log-normal  random  variables 
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Iteration 

Figure  10:  MSE  vs.  Coverage  (Rc=2,Rs=l) 


Nr  Nr 

I  =  ^2  ^  10fiM<«W10  (11) 

fe=l  fc=l 

where  the  Qk(dBm)  are  Gaussian  random  variables  with  mean  nnk(dBm)  and  variance  a^k ,  and 

the  flk  =  10 nk(dBm)/io  are  the  log-normal  random  variables.  Unfortunately,  there  is  no  known 
closed  form  expression  for  the  probability  density  function  (pdf)  of  the  sum  of  multiple  (Nj  >2) 
log-normal  random  variables.  However,  there  is  a  general  consensus  that  the  sum  of  independent 
log-normal  random  variables  can  be  approximated  by  another  log-normal  random  variable  with 
appropriately  chosen  parameters.  That  is, 

N] 

I  =  y  10 nfc(dBm)/l°  ~  10Z(dBm)/lO  _  J  (12) 

fc= 1 

where  Z^dBm^  is  a  Gaussian  random  variable  with  mean  Hz(dBrn)  and  variance  cr|.  The  problem 
is  to  determine  nZ{dBm)  and  variance  a\  in  terms  of  the  Hnk(dBm)  and  variance  ,  k  =  1,  ...,1V/. 
Several  methods  have  been  suggested  in  the  literature  to  solove  this  problem  including  those  by 
Fenton,  Schwartz  and  Yen,  and  Farley.  Each  of  these  methods  provides  varying  degrees  of  accuracy 
over  specified  ranges  of  the  shadow  standard  deviation  <tq,  the  sum  7,  and  the  number  of  interferes 


15 


Nj. 

B  Fenton- Wilkinson  Method 

With  the  Fenton- Wilkinson  method,  the  mean  Mz(dBm)  and  variance  a\  of  Z^dBm)  are  obtained  by 
matching  the  first  two  moments  of  the  sum  I  with  the  first  two  moments  of  the  approximation  /. 
To  derive  the  appropriate  moments,  it  is  convenient  to  use  natural  logarithms.  We  write 


Clk  =  lO^B™)/10  =  ee^k(dBm)  —  g &k 

where  e  =  (Znl0)/10  =  0.23026  and  Clk  =  dlk[dBm).  Note  that  n&k  =  emk(iBm)  and  oj  = 
f2(7hk-  The  nth  moment  of  the  log-normal  random  variable  can  be  obtained  from  the  moment 
generating  function  of  the  Gaussian  random  variables  as 


E[n%]  =  E[enClk]  =  enflClk+{1/2)n2ck 


(14) 


To  find  the  appropriate  moments  for  the  log-normal  approximation  we  can  use  (14)  and  equate 
the  first  two  moments  on  both  sides  of  the  equation 


Ni 

I  =  ^Te^k  ^  ez  =  I  (15) 

fc= 1 

where  Z  =  tZ^Bm)-  For  example,  suppose  that  fi*.,  k  =  1 ,  ...,1V/  have  mean  fi ^  ,  k  =  1 ,  ...,1V/ 
and  identical  variances  a? .  Identical  variances  are  often  assumed  because  the  standard  deviation 
of  log-normal  shadowing  is  largely  independent  of  the  radio  path  length.  Equating  the  means  on 
both  sides  of  (15) 


Ni 

M/  =  E[I)  =  J2  E[eQk]  =  E[ez]  =  E[I ]  =  (16) 

k= 1 

gives  the  result 
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^e^^e(1/2)4=e^+(1/2V| 


(17) 


Likewise,  we  can  equate  the  variances  on  both  sides  of  (15)  under  the  assumption  that  the 
k  =  1 .....  iV/  are  independent 


*2I  =  E[I2}-tf  =  E{i2\  =  aj  (18) 

giving  the  result 


-  1)  =  e2^eCTl(eCTl  -  1) 


(19) 


By  squaring  each  side  of  (17)  and  dividing  each  side  of  resulting  equation  by  the  respective 
side  of  (19)  We  can  solve  for  oflatz  in  terms  of  the  known  values  of  ,  k  =  1,  ...,7V/  and  a2. 

k 

Afterwards,  can  be  obtained  from  (17).  This  procedure  yields  the  following  solution: 


V'Z 


°jrA 


+  In 


e^k 


(20) 


a 


2 

z 


=  In 


(E£i«^)a 


(21) 


Finally,  Hz{dBm)  -  e  Vz  311(1  °i  =  €_2(T|- 

The  accuracy  of  this  log-normal  approximation  can  be  measured  in  terms  of  how  accurately 
the  first  two  moments  of  I{dB)  —  10/opioI  are  estimated,  and  how  well  the  cumulative  distribution 
function  (cdf)  of  I(dB)  1S  described  by  a  Gaussian  cdf.  In  problems  relating  to  the  co-channel 
interference  outage  in  celluar  radio  systems,  we  are  usually  interested  in  the  tails  of  both  the 
complementary  distribution  function  (cdfc)  Ff  =  P(I  >  x)  and  the  cdf  Fj( x)  =  1  -Fp  =  P(I  <  x). 
In  this  case,  we  are  interested  in  the  accuracy  of  the  approximation 


for  large  and  small  values  of  x.  It  will  be  shown  later  that  the  Fenton- Wilkinson  method  can 

approximate  the  tails  of  the  cdf  and  cdfc  functions  with  good  accuracy. 
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Abstract — Prolonging  the  lifetime  of  energy  constrained  wire¬ 
less  sensor  networks  is  a  crucial  challenge  in  sensor  network 
research.  In  this  paper,  we  present  a  novel  approach  based 
on  fuzzy  logic  theory  to  analyze  the  lifetime  performance  of  a 
sensor  network.  We  demonstrate  that  a  type-2  fuzzy  membership 
funetion(MF),  i.e.,  a  Gaussian  MF  with  uncertain  variance  is 
most  appropriate  to  model  a  single  node  lifetime  in  wireless 
sensor  networks.  In  our  research,  we  concern  two  basic  sensor 
placement  schemes:  square-grid  and  hex-grid.  Two  fuzzy  logic 
systems(FLSs):  a  singleton  type-1  FLS  and  an  interval  type-2 
FLS  are  designed  to  perform  lifetime  estimation  of  the  entire 
sensor  network.  Simulation  results  show  that  the  interval  type-2 
FLS  in  which  the  antecedent  membership  functions  are  modeled 
as  type-2  fuzzy  sets  outperforms  the  singleton  type-1  FLS  when 
the  single  node  lifetime  behaves  the  nature  of  Gaussian  MFs  with 
uncertain  standard  deviation. 

I.  Introduction 

Wireless  sensor  network  represents  significant  advantages 
over  conventional  personnel-rich  methods  of  long  term  data 
collection  and  monitoring.  Starting  in  the  1980s,  networked 
microsensors  technology  has  been  widely  used  in  military 
applications.  Examples  of  these  applications  include  the  Co¬ 
operative  Engagement  Capability  (CEC),  Remote  Battlefield 
Sensor  System  (REMBASS),  Advanced  Deployable  System 
(ADS)  and  Tactical  Remote  Sensor  System  (TRSS)  [1].  Al¬ 
though  a  majority  of  military  applications  deploy  sensor  nodes 
randomly,  a  number  of  other  applications  require  the  place¬ 
ment  of  sensors  at  desired  locations.  Such  placement-friendly 
sensor  networks  are  developed  rapidly  for  infrastructure  se¬ 
curity,  environment  and  habitat  monitoring,  traffic  control  etc 
[1].  A  real-world  32  node  habitat  monitoring  sensor  network 
system  was  deployed  on  a  small  island  off  the  coast  of  Maine 
to  study  nesting  patterns  of  Petrels  [2].  Other  applications 
involve  the  use  of  sensors  in  buildings  for  environmental 
monitoring  which  may  include  chemical  sensing  and  detection 
of  moisture  problems.  Structural  monitoring  [3]  and  inventory 
control  are  some  other  applications  of  such  networks. 

Wireless  sensor  network  consists  of  certain  amount  of  small 
and  energy  constrained  nodes.  Basic  components  of  sensor 
node  include  a  single  or  multiple  sensor  modules,  a  wireless 
transmitter-receiver  module,  a  computational  module  and  a 
power  supply  module.  Such  networks  are  normally  deployed 
for  data  collection  where  human  intervention  after  deployment, 
to  recharge  or  replace  node  batteries  may  not  be  feasible,  re¬ 


sulting  in  limited  network  lifetime.  Most  applications  have  pre¬ 
specified  lifetime  requirements,  for  instance  the  application 
in  [2]  has  a  lifetime  requirement  of  at  least  9  months.  Thus 
estimation  of  lifetime  of  such  networks  prior  to  deployment 
becomes  a  necessity.  Prior  works  on  evaluating  lifetime  have 
considered  networks  where  sensor  nodes  are  randomly  de¬ 
ployed.  [4]  gives  the  upper  bound  on  lifetime  that  any  network 
with  the  specified  number  of  randomly  deployed  nodes,  source 
behavior  and  energy  can  reach  while  [6]  discusses  the  upper 
bounds  on  lifetime  of  networks  with  cooperative  cell  based 
strategies.  Network  liftime  of  fixed  deployment  schemes  are 
recently  studied  in  [7].  Jain  and  Liang  [7]  discovers  that  in 
wireless  sensor  networks,  a  single  node  lifetime  behaves  the 
nature  of  normal  distribution  which  brings  the  first  light  of 
exploring  the  network  lifetime  behavior  given  the  knowledge 
of  nodes  lifetimes. 

In  this  paper,  we  deal  with  the  issue  of  lifetime  analysis  and 
estimation  for  wireless  sensor  networks  in  which  the  sensor 
nodes  are  deployed  at  desired  locations.  We  propose  to  apply 
interval  type-2  fuzzy  logic  systems  (FLSs)  for  network  lifetime 
analysis  and  estimation.  Our  approach  is  entirely  different 
from  all  prior  research.  We  demonstrate  that  a  type-2  fuzzy 
membership  function(MF),  i.e.,  a  Gaussian  MF  with  uncertain 
variance  is  most  appropriate  to  model  a  single  node  lifetime 
in  wireless  sensor  networks.  Two  fuzzy  logic  systems(FLSs): 
a  singleton  type-1  FLS  and  an  interval  type-2  FLS  are  con¬ 
structed  to  perform  lifetime  estimation  of  the  entire  sensor 
network.  In  our  work,  we  concern  two  basic  sensor  placement 
schemes:  square-grid  and  hex-grid.  We  believe  that  these  two 
schemes  can  serve  as  basis  for  evaluating  more  complex 
schemes  for  their  lifetime  performance  prior  to  deployment 
and  help  justify  the  deployment  costs. 

The  rest  of  this  paper  is  organized  as  follows.  In  section  II, 
we  detail  two  sensor  deployment  schemes  used  for  this  study 
and  the  basic  concepts  on  network  coverage,  connectivity 
and  lifetime.  Section  III  gives  an  overview  of  interval  type- 
2  fuzzy  logic  systmes.  In  section  IV,  we  demonstrate  that  a 
single  node  lifetime  can  be  modeled  with  Gaussian  MFs.  We 
found  that  a  type-2  Gaussian  MF  with  uncertain  variance  is 
most  appropriate  to  model  a  single  node  lifetime  in  wireless 
sensor  networks.  We  apply  this  knowledge  and  design  an 
interval  type-2  FLS  in  section  V  to  analyze  network  lifetime. 
A  singleton  type-1  FLS  is  constructed  as  well  for  performance 


comparison.  Simulation  results  and  discussions  are  also  pre¬ 
sented  in  section  V.  Section  VI  concludes  this  paper. 

II.  Preliminaries 

A.  Basic  Model  and  Assumptions 

Consider  identical  wireless  sensor  nodes  placed  in  a  square 
sensor  field  of  area  A.  All  nodes  are  deployed  with  equal  en¬ 
ergy.  Each  sensor  is  capable  of  sensing  events  up  to  a  radius  rs, 
the  sensing  range.  Communication  range  rc  is  defined  as  the 
distance  beyond  which  the  transmitted  signal  is  received  with 
signal  to  noise  ratio  {SNR)  below  the  acceptable  threshold 
level.  In  this  paper,  We  assume  the  communication  range  rc 
to  be  equal  to  the  sensing  range  rs. 

Direct  communication  between  two  sensor  nodes  is  possible 
only  if  their  distance  of  separation  r  is  such  that  r  <rc.  We 
call  such  nodes  neighbors.  Communication  between  a  sensor 
node  and  its  non-neighboring  node  is  achieved  via  peer-to- 
peer  communication.  Thus  the  maximum  allowable  distance 
between  two  nodes  who  wish  to  communicate  directly  is 
rmax  =  rc  =  rs.  A  network  is  said  to  be  deployed  with 
minimum  density  when  the  distance  between  its  neighboring 
nodes  is  v  =  Ty^iax' 

B .  Placement  Schemes 

The  simplest  placement  schemes  involve  regular  placement 
of  nodes  such  that  each  node  in  the  network  has  the  same  num¬ 
ber  of  neighbors.  We  arrive  at  two  basic  placement  schemes 
by  considering  cases  where  each  sensor  nodes  has  four  and 
three  neighbors.  This  leads  us  to  the  square-grid  and  hex-grid 
placement  schemes  shown  in  Fig.  1(a)  and  (b)  respectively. 

A  sensor  node  placement  scheme  that  uses  two  neighbors 
per  sensor  node  has  been  described  in  [13].  We  believe  that 
these  three  elementary  placement  schemes  can  serve  as  basis 
for  other  placement  schemes,  because  a  placement  scheme  of 
any  complexity  can  be  decomposed  into  two-neighbor,  three- 
neighbor  and  four-neighbor  groups.  Both  grids  shown  have  the 
same  number  of  nodes 1  and  nodes  in  both  grids  are  equidistant 
from  their  respective  neighbors  (with  distance  of  separation  r). 

C.  Coverage  and  Connectivity 

Coverage  and  connectivity  are  two  important  performance 
metrics  of  networks  and  hence  a  discussion  on  them  becomes 
imperative  before  the  lifetime  of  the  network  can  be  defined. 

Coverage  scales  the  adequacy  with  which  the  network 
covers  the  sensor  field.  A  sensor  with  sensing  range  rs  is 
said  to  cover  or  sense  a  circular  region  of  radius  r$  around  it. 
If  every  point  in  the  sensor  field  is  within  distance  rs  from 
at  least  one  sensor  node,  then  the  network  is  said  to  provide 
complete  or  100%  coverage. 

Connectivity  scales  the  adequacy  with  which  nodes  are  able 
to  communicate  with  their  peers.  One  of  the  strengths  of  sensor 
networks  arises  from  their  ability  to  aggregate  data  collected 

lrThe  Hex-grid  has  lower  density  than  the  Square-grid.  With  36  nodes 
deployed,  the  network  with  Square-grid  covers  an  area  of  25 P,  and  the 
Hex-grid  covers  an  area  of  48 A,  almost  double  that  of  the  square  grid. 


Fig.  1.  Placement  Schemes  for  a  36  node  sensor  network:  (a)  Square-Grid 
(b)  Hex-Grid. 


from  different  sensor  nodes.  This  requires  adequate  commu¬ 
nication  between  sensor  nodes.  Any  node  should  be  able  to 
communicate  with  any  other  node  for  proper  functioning  of 
the  network.  If  a  single  node  gets  isolated  due  to  failure  of  all 
its  neighbors,  it  will  be  unable  to  communicate  with  the  rest 
of  the  network.  If  a  large  number  of  nodes  fail  due  to  lack  of 
energy,  a  part  of  the  network  may  get  completely  disconnected 
from  the  rest. 

In  our  analysis  we  require  the  network  to  provide  complete 
coverage  and  connectivity.  We  give  equal  importance  to  both 
parameters  and  declare  the  network  nonfunctional  if  either  of 
them  falls  below  their  desired  levels. 

D.  Lifetime 

The  basic  definition  of  lifetime,  or  more  precisely  the  post¬ 
deployment  active  lifetime  of  a  network  is  the  time  measured 
from  deployment  until  network  failure.  Based  on  the  levels 
of  coverage  and  connectivity  required  to  deem  a  network 
functional,  network  failure  can  be  interpreted  in  different  ways. 
Since  only  complete  coverage  and  connectivity  are  acceptable 
to  us,  network  failure  corresponds  to  the  first  loss  of  coverage 
or  connectivity. 

In  this  paper,  we  concentrate  on  finding  the  minimum 
lifetime  of  a  network,  the  worst  case  scenario.  Consider  the 
square-grid  and  the  hex-grid  deployed  with  minimum  density. 
Both  schemes  survive  the  failure  of  a  single  node  without  loss 
of  either  connectivity  or  coverage  however  failure  of  any  two 
neighboring  nodes  causes  loss  of  coverage  and  hence  network 
failure  as  indicated  in  Figs  2(a)  and  (b). 

Thus  the  minimum  number  of  node  failures  that  cause 


network  failure  is  two  and  these  two  nodes  must  be  adjacent  to 
each  other  (neighbors).  A  network  may  undergo  multiple  node 
failures  and  still  be  connected  and  covered  if  any  of  the  failed 
nodes  are  not  neighbors.  But  the  absolute  minimum  number 
of  node  failures  that  can  cause  network  failure  is  two. 


(a) 


□  -  □ 

(b) 


Fig.  2.  Loss  of  coverage  due  to  failure  of  two  neighboring  nodes:  (a)Square- 
grid:  Failure  of  nodes  20  and  21  causes  loss  of  coverage.  (b)Hex-grid:  Failure 
of  nodes  20  and  25  causes  loss  of  coverage. 


III.  Introduction  to  Interval  Type-2  Fuzzy  Logic 
Systems 

Figure  3  shows  the  structure  of  a  type-2  FLS.  It  is  very 
similar  to  the  structure  of  a  type-1  FLS  [16].  For  a  type-1  FLS, 
the  output  processing  block  only  contains  the  defuzzifier.  We 
assume  that  the  reader  is  familiar  with  type-1  FLSs,  so  that 
here  we  focus  only  on  the  similarities  and  differences  between 
the  two  FLSs. 

The  fuzzifier  maps  the  crisp  input  into  a  fuzzy  set.  This 
fuzzy  set  can,  in  general,  be  a  type-2  set. 

In  the  type-1  case,  we  generally  have  “IF-THEN”  rules, 
where  the  Ith  rule  has  the  form  uRl  :  IF  xi  is  ¥[  and  x2 
is  F2  and  •  •  •  and  xp  is  ¥lp,  THEN  y  is  G*”,  where:  XiS  are 
inputs;  F-s  are  antecedent  sets  (i  =  1, . . .  ,p);  y  is  the  output; 
and  G*s  are  consequent  sets.  The  distinction  between  type-1 
and  type-2  is  associated  with  the  nature  of  the  membership 
functions,  which  is  not  important  while  forming  rules;  hence, 


TYPE-2  FUZZY  LOGIC  SYSTEM 


Fig.  3.  The  structure  of  a  type-2  FLS.  In  order  to  emphasize  the  importance 
of  the  type-reduced  set,  we  have  shown  two  outputs  for  the  type-2  FLS,  the 
type-reduced  set  and  the  crisp  defiizzified  value. 


the  structure  of  the  rules  remains  exactly  the  same  in  the  type- 
2  case,  the  only  difference  being  that  now  some  or  all  of  the 
sets  involved  are  of  type-2;  so,  the  Ith  rule  in  a  type-2  FLS 
has  the  form  “Rl  :  IF  xx  is  F^  and  x2  is  P2  and  ■  ■  •  and  xp 
is  Fp,  THEN  y  is  G*” 

In  the  type-2  case,  the  inference  process  is  very  similar  to 
that  in  type-1.  The  inference  engine  combines  rules  and  gives 
a  mapping  from  input  type-2  fuzzy  sets  to  output  type-2  fuzzy 
sets.  To  do  this,  one  needs  to  find  unions  and  intersections  of 
type-2  sets,  as  well  as  compositions  of  type-2  relations. 

In  a  type-1  FLS,  the  defuzzifier  produces  a  crisp  output 
from  the  fuzzy  set  that  is  the  output  of  the  inference  engine, 
i.e.,  a  type-0  (crisp)  output  is  obtained  from  a  type-1  set. 
In  the  type-2  case,  the  output  of  the  inference  engine  is  a 
type-2  set;  so,  we  use  “extended  versions”  (using  Zadeh’s 
Extension  Principle  [15])  of  type-1  defuzzification  methods. 
This  extended  defuzzification  gives  a  type-1  fuzzy  set.  Since 
this  operation  takes  us  from  the  type-2  output  sets  of  the  FLS 
to  a  type-1  set,  we  call  this  operation  “type-reduction”  and  the 
type-reduced  set  so  obtained  a  “type-reduced  set”. 

To  obtain  a  crisp  output  from  a  type-2  FLS,  we  can 
defuzzify  the  type-reduced  set.  The  most  natural  way  of  doing 
this  seems  to  be  by  finding  the  centroid  of  the  type-reduced 
set;  however,  there  exist  other  possibilities  like  choosing  the 
highest  membership  point  in  the  type-reduced  set. 

General  type-2  FLSs  are  computationally  intensive,  because 
type-reduction  is  very  intensive.  Things  simplify  a  lot  when 
secondary  membership  functions  (MFs)  are  interval  sets  (in 
this  case,  the  secondary  memberships  are  either  0  or  1). 
When  the  secondary  MFs  are  interval  sets,  we  call  the  type-2 
FLSs  “interval  type-2  FLSs”.  In  [21],  we  proposed  the  theory 
and  design  of  interval  type-2  fuzzy  logic  systems  (FLSs). 
We  proposed  an  efficient  and  simplified  method  to  compute 
the  input  and  antecedent  operations  for  interval  type-2  FLSs, 
one  that  is  based  on  a  general  inference  formula  for  them. 
We  introduced  the  concept  of  upper  and  lower  membership 
functions  (MFs)  and  illustrate  our  efficient  inference  method 
for  the  case  of  Gaussian  primary  MFs.  We  also  proposed  a 
method  for  designing  an  interval  type-2  FLS  in  which  we  tune 
its  parameters. 

We  have  developed  theory  and  design  methods  for  the  most 
useful  kind  of  type-2  fuzzy  logic  system  (FLSs),  interval  type- 


Fig.  4.  Single  Node  Lifetime  Distribution 


2  FLSs  [21],  and  have  applied  them  to  a  number  of  very 
interesting  applications,  such  as 

1)  Fading  channel  equalization  [22]  and  co-channel  in¬ 
terference  elimination  [23].  The  channel  states  in  a 
fading  channel  or  channel  with  co-channel  interferences 
are  uncertain,  and  we  validated  that  an  interval  type- 
2  fuzzy  set,  Gaussian  primary  membership  function 
with  uncertain  mean,  can  be  used  to  represent  such 
uncertainties. 

2)  Network  video  traffic  modeling  and  classification  [19]. 
MPEG  variable  bit  rate  (VBR)  traffic  are  very  bursty. 
We  validated  that  the  I,  P,  and  B  frame  sizes  are  log¬ 
normal  with  fixed  mean  and  uncertain  variance,  so  an 
interval  type-2  fuzzy  sets  can  be  used  to  model  the  bursty 
video  traffic  and  an  interval  type-2  fuzzy  logic  system 
with  such  type-2  fuzzy  set  are  demonstrated  performing 
much  better  than  a  Bayesian  classifier. 

3)  Connection  admission  control  for  ATM  network  [20]. 
Connection  admission  control  is  actually  a  decision 
making  problem.  Different  factors  such  as  incoming 
real-time  video/audio  packet  sizes,  non-real  time  packet 
sizes,  the  buffer  sizes  are  uncertain.  We  applied  an 
interval  type-2  fuzzy  logic  to  handle  these  uncertainties, 
and  achieved  very  good  performance. 

IV.  Modeling  Node  Lifetime  with  Gaussian 
Membership  Functions 

Wireless  sensor  nodes  are  severely  energy  constrained  due 
to  their  compact  form.  To  increase  the  lifetime  of  sensor  net¬ 
works,  hardware  design  and  protocol  approaches  for  different 
layers  must  take  energy  efficiency  into  account.  However,  a 
foundamental  question  -  “what  is  the  nature  of  sensor  network 
lifetime?”  has  not  been  answered  yet.  Since  the  lifetimes  of 
individual  nodes  are  not  constants  but  random  variables,  it 
follows  that  the  network  lifetime  is  also  a  random  variable. 
Recent  research  by  Jain  and  Liang  [7]  discovered  that  in 
a  wireless  sensor  network  where  the  workloads  are  well- 
balanced,  a  single  node  lifetime  behaves  the  nature  of  normal 
distribution  as  demonstrated  in  Fig  4. 

In  this  paper,  we  are  first  interested  in  setting  up  fine 
membership  functions  (MFs)  for  single  node  lifetime.  From 


Fig.  5.  Type-2  Gaussian  MF  with  uncertain  standard  deviation 
TABLE  I 

MEAN  AND  STD  VALUES  FOR  SEVEN  SEGMENTS  AND  THE 
ENTIRE  NODE  LIFETIME,  AND  THEIR  NORMALIZED  STD 


Node  Lifetime  Data 

Mean 

Std 

Segment  1 

1027.4 

30.182 

Segment  2 

1028.9 

29.819 

Segment  3 

1026.3 

30.798 

Segment  4 

1028.7 

30.917 

Segment  5 

1028 

29.944 

Segment  6 

1027 

29.975 

Segment  7 

1027.9 

30.306 

Entire  Data  Set 

1027.7 

30.292 

Normalized  STD 

0.00082783 

0.013105 

the  original  data  of  single  node  lifetime  shown  in  Table  I, 
we  decompose  the  whole  data  sets  into  seven  segments,  and 
compute  the  mean  fa  and  standard  deviation  a*  of  node 
lifetime  for  each  segment,  i  =  1,2,  ...7.  The  mean  fi  and 
standard  deviation  a  for  of  the  entire  data  set  is  also  computed. 
We  are  also  interested  to  know  which  value  -  mean  \ii  or 
standard  deviation  <7*  varies  more.  We  first  normalize  the 
mean  fa  and  standard  deviation  <7*  of  each  segment  using 
fa/fi  and  ail  a.  Then  we  compute  the  standard  deviation  of 
their  normalized  values  am  and  astd .  Results  are  presented  in 
Table  I. 

From  the  last  row  of  Table  I,  we  see  that  am  «  astd 
which  means  standard  deviation  a i  varies  much  more  than  the 
mean  value  fa.  Therefore  we  conclude  that  if  the  single  node 
lifetime  follows  normal  distribution,  it  is  most  appropriate 
to  be  modeled  as  a  Gaussian  MF  with  uncertain  standard 
deviation.  One  example  of  type-2  Gaussian  MF  with  uncertain 
standard  deviation  is  shown  in  Fig  5.  This  result  also  justifies 
the  use  of  the  Gaussian  MFs  to  model  network  lifetime  in 
section  V. 

V.  Network  Lifetime  Analysis  and  Estimation 
Using  Interval  Type-2  Fuzzy  Logic  Systems 

We  are  now  ready  to  evaluate  the  network  lifetime  using 
interval  type-2  fuzzy  logic  systems  (FLSs).  We  apply  reliabil¬ 
ity  theory  to  design  fuzzy  rules.  In  the  following  section,  we 
first  treat  the  basics  of  reliability  theory  before  extracting  the 


(b) 


Fig.  6.  Reliability  block  diagrams  (RBD)  for  a  system  of  two  components. 
(a)RBD  with  series  connected  components.  (b)RBD  with  parallel  connected 
components. 


knowledges  for  fuzzy  rules.  An  interval  type-2  FLS  is  then  set 
up  to  perform  lifetime  analysis  and  estimation. 

A.  Fuzzy  Rules  Design  Using  Reliability  Theory 

Rules  are  the  heart  of  a  FLS,  and  may  be  provided  by 
experts  or  can  be  extracted  from  numerical  data.  In  either 
case,  the  rules  that  we  are  interested  in  can  be  expressed  as  a 
collection  of  IF-THEN  statements,  e.g., 

IF  the  noise  of  the  wireless  channel  is  high,  THEN  the 
quality  of  the  received  signal  is  poor. 

For  the  lifetime  issue  studied  in  this  paper,  reliability 
theory  provides  a  feasible  method  to  design  fuzzy  rules.  To 
understand  this,  we  introduce  the  reliability  block  diagram 
(RBD).  RBD  is  a  graphical  representation  of  the  components 
of  the  system,  and  provides  a  visual  representation  of  the 
way  components  are  reliability-wise  connected.  Thus  the  effect 
of  the  success  or  failure  of  a  component  on  the  system 
performance  can  be  evaluated. 

Consider  a  system  with  two  components.  If  this  system  is 
such  that  a  single  component  failure  can  render  the  system 
nonfunctional,  then  we  say  that  the  components  are  reliability- 
wise,  connected  in  series.  If  the  system  fails  only  when  both 
its  components  fail,  then  we  say  that  the  components  are 
reliability-wise  connected  in  parallel.  Note  that  the  physical 
connection  between  the  component  may  or  may  not  be  dif¬ 
ferent  from  their  reliability-wise  connection.  The  RBD's  for 
both  cases  are  given  in  Figs. 6.  Any  complex  system  can  be 
realized  in  the  form  of  a  combination  of  blocks  connected  in 
series  and  parallel. 

In  our  analysis,  the  wireless  sensor  network  is  the  system 
under  consideration  and  the  sensor  nodes  are  the  components 
of  the  system.  We  demonstrate  bellow  how  the  knowledge  is 
extracted  for  fuzzy  rules  design  referring  to  the  two  examples 
in  Fig.6. 

Example  1:  Set  up  Fuzzy  Rules  for  Parallel  System  in 
Fig.  6(a) 

In  the  parallel  system,  the  system  (network)  fails  only  when 
both  components  (sensor  nodes)  fail.  The  rules  can  be  set  up 
as  one  example  shown  bellow: 


IF  the  remaining  battery  level  of  component  1  (sensor  node 
1)  is  low  and  the  remaining  battery  level  of  component  2 
(sensor  node  2)  is  moderate ,  THEN  the  lifetime  of  the 
system  (network)  is  moderate . 

Example  2:  Set  up  Fuzzy  Rules  for  Series  system  in 
Fig.  6(b) 

In  the  series  system,  the  system  (network)  fails  when  either 
component  fails.  The  rules  can  be  set  up  as  one  example  shown 
bellow: 

IF  the  remaining  battery  level  of  component  1  (sensor  node 
1)  is  low  and  the  remaining  battery  level  of  component  2 
(sensor  node  2)  is  moderate ,  THEN  the  lifetime  of  the 
system  (network)  is  low. 

Note  that  the  parallel  and  series  systems  are  the  two  basic 
ways  to  model  two  sensor  nodes.  A  wireless  sensor  network 
consisting  of  multiple  sensor  nodes  can  be  first  represented 
in  the  reliability  block  diagram.  Fuzzy  rules  are  then  set  up 
based  on  the  above  examples. 

B.  Simulation  and  Discussion 

In  our  simulations,  interval  type-2  FLSs  are  constructed 
separately  for  square-grid  netowrk  and  hex-grid  network.  We 
take  the  remaining  energy  level  of  each  individual  sensor 
node  as  input  to  the  interval  type-2  FLS  and  the  output 
is  the  estimated  network  lifetime.  As  we  have  discussed  in 
section  IV,  we  choose  type-2  Gaussian  with  uncertain  standard 
deviation  as  the  membership  functions  for  both  antecedent  and 
consequent.  The  linguistic  variables  to  represent  the  antecedent 
-  remaining  energy  level  are  divided  into  three  levels:  high, 
moderate  and  low  and  the  consequent  -  estimated  network  life 
is  divided  into  five  levels:  very  high ,  high,  moderate,  low  and 
very  low. 

Simulations  are  performed  for  both  squre-grid  and  hex-grid 
sensor  networks.  In  both  cases,  36  nodes  are  deployed  and 
the  distance  between  neighboring  nodes  is  assumed  to  be  the 
same.  All  sensor  nodes  are  initialized  with  maximum  energy 
level  of  10  Joules.  We  assume  the  workload  of  the  entire 
network  is  well-balanced,  therefore  energy  dissipation  of  the 
single  sensor  node  can  be  figured  as  linear  most  of  the  alive 
time.  Our  simulations  are  based  on  N  =  600  lifetime  data 
from  the  above  well-balanced  sensor  network.  The  first  300 
data  are  for  training  and  the  remaining  300  data  are  for  testing. 
After  training,  the  rules  are  fixed  and  we  test  the  interval  type- 
2  by  evaluating  the  root  mean  square  errors  (RMSE)  between 
the  defuzzified  output  of  FLS  and  the  real  lifetime  data.  We 
compare  our  interval  type-2  FLSs  with  singleton  type-1  FLSs. 
For  both  FLS  schemes,  we  run  200  Monte-Carlo  realizations 
and  for  each  realization,  each  FLS  is  tuned  using  a  simple 
steepest-descent  algorithm  for  six  epochs.  Simulation  results 
are  averaged  over  all  200  Monte-Carlo  realizations. 

1)  Square  Grid:  As  defined  in  Section  II-D,  the  minimum 
network  lifetime  is  the  time  to  failure  of  any  two  neighboring 
nodes.  We  know  that  the  failure  of  any  single  node  does  not 
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Fig.  7.  RBD  of  a  single  node  in  a  square  grid.  Nodes  belonging  to  region- 
1  are  modeled  as  block-1  and  nodes  belonging  to  region-2  are  modeled  as 
block-2.  The  network  RBD  consists  of  tyNmin  -  l)2  block- l’s  in  series 
with  2(y/rNrnin  —  1)  block-2 ’s 


cause  network  failure.  The  failure  of  any  node  coupled  with  the 
failure  of  any  of  its  neighbors  causes  network  failure.  Using 
this  definition  we  build  the  RBD  for  the  square-grid  as  shown 
in  Fig  7. 

Fig  7  shows  the  RBD  block  for  a  single  node  in  the  network. 
A  node  can  be  modeled  in  two  ways  depending  on  its  position 
in  the  sensor  field.  This  distinction  based  on  its  position  is 
made  due  to  a  simple  observation  that  nodes  at  the  right  edge 
of  the  sensor  field  (region-2)  do  not  have  any  right  neighbor 
(node  b)  as  opposed  to  nodes  in  region- 1.  Also,  nodes  at 
the  bottom  edge  of  the  sensor  field  (region-2)  do  not  have  a 
bottom  neighbor  (node  c)  as  opposed  to  the  nodes  in  region- 
1.  Note  that  as  every  node  in  a  square-grid,  node  a  has  four 
neighbors,  but  its  relationship  with  only  two  neighbors  is 
modeled  in  its  RBD  block.  This  is  because  the  relationship 
with  the  other  two  neighbors  will  be  modeled  when  their 
RBD  blocks  are  constructed.  If  this  is  not  followed  then  the 
relationship  between  every  node-neighbor  pair  will  be  modeled 
twice. 

In  this  square  grid  network  shown  in  Fig  7,  we  classify  three 
antecedents  based  on  the  RBD  of  block- 1: 

•  The  remaining  battery  level  of  node  a. 

•  The  minimum  remaining  battery  level  of  node  b  and  c. 

•  The  remaining  battery  level  of  node  d. 

We  set  up  27  rules  for  this  FLS  because  every  antecedent 
has  3  fuzzy  sub-sets  and  there  are  3  antecedents. 

Two  antecedents  are  chosen  based  on  the  RBD  of  block-2 
in  Fig  7  and  total  9  rules  are  constructed  in  this  case. 

•  The  remaining  battery  level  of  node  x. 

•  The  remaining  battery  level  of  node  y . 

Let  Nmin  be  the  number  of  sensor  nodes  required  to  be 
deployed  with  minimum  density.  The  network  RBD  consists  of 
(V^mm  —  1)2  block-l’s  and  2(^Nmin—l)  block-2’s  in  series, 
the  whole  square  grid  network  can  actually  be  decomposed 
into  multiple  blocks  serial  connected  together  and  the  method 
of  setting  up  rules  can  be  applied  as  well.  Simulation  results 
of  RMSE  is  shown  in  Fig.  8. 

2)  Hex-Grid:  The  analysis  for  the  hex-grid  is  carried  out  on 
the  same  lines  as  that  of  the  square-grid.  Fig  2  (b)  shows  that 
as  in  the  case  of  a  square  grid,  two  neighboring  node  failures 
cause  network  failure.  The  RBD  block  of  a  single  node  is 
shown  in  Fig.  9. 


Fig.  8.  Square-Grid:  The  RMSE(for  the  test  data)  for  two  FLS  approaches 
averaged  over  200  Monte-Carlo  realization 
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Fig.  9.  RBD  block  for  a  single  node  in  the  Hex-grid:  The  network  RBD 
consists  of  N/2  such  blocks  in  series. 


Fig.  10.  Hex-Grid:  The  RMSE(for  the  test  data)  for  two  FLS  approaches 
averaged  over  200  Monte-Carlo  realization 


Since  the  relation  between  a  node  and  all  of  its  neighbors  is 
modeled  by  its  corresponding  RBD  block,  the  RBD  block’s  for 
the  neighbors  is  not  constructed  as  this  causes  the  relationship 
between  the  nodes  to  be  considered  twice.  In  this  case,  we 
classify  two  antecedents  and  construct  9  rules. 

•  The  remaining  battery  level  of  node  a. 

•  The  minimum  remaining  battery  level  of  node  6,  c  and 
d. 

Since  Nmin/ 2  such  blocks  connected  in  series  represent  the 
network,  the  whole  hex  grid  network  can  be  decomposed  the 
same  way  as  in  square  grid  networks.  Simulation  results  of 
RMSE  is  shown  in  Fig.  10. 

Observe  Fig.8  and  Fig.  10,  interval  type-2  FLSs  outperform 
singleton  type-1  FLSs.  This  results  shows  that  the  interval 
type-2  FLSs  are  more  feasible  for  real-time  energy  estimation. 


VI.  Conclusion  and  future  works 

In  this  paper,  we  describe  a  new  method  based  on  fuzzy 
logic  theory  to  analyze  and  estimate  the  network  lifetimes  for 
wireless  sensor  networks.  Our  approach  is  illuminated  by  the 
discovery  that  a  single  node  lifetime  behaves  the  nature  of 
normal  distribution.  However, we  deem  that  if  the  single  node 
lifetime  follows  normal  distribution,  it  is  most  appropriate 
to  be  modeled  as  a  Gaussian  MF  with  uncertain  standard 
deviation.  We  then  set  up  the  interval  type-2  FLSs  for  energy 
estimation  and  evaluate  their  performance  using  real  lifetime 
data.  Simulation  results  justified  the  feasibility  of  applying 
type-2  FLSs  into  wireless  sensor  network  lifetime  analysis. 
Interval  type-2  FLSs  provides  a  way  to  handel  knowledge 
uncertainty.  We  believe  that  our  approaches  opens  up  a  new 
vision  for  research  on  sensor  network  lifetimes. 

Our  future  work  will  focus  on  lifetime  evaluation  under  the 
circumstances  that  the  task  scheduling  is  variable  and  how  the 
estimated  network  lifetime  could  be  used  to  accomodate  the 
scheduling  change. 
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Abstract — In  this  paper,  we  present  a  new  approach  for 
sensed  signal  strength  forecasting  in  wireless  sensors  using 
interval  type-2  fuzzy  logic  system  (FLS).  We  show  that  a 
type-2  fuzzy  membership  function,  i.e.,  a  Gaussian  MF  with 
uncertain  mean  is  most  appropriate  to  model  the  sensed  signal 
strength  of  wireless  sensors.  We  demonstrate  that  the  sensed 
signals  of  wireless  sensors  are  self-similar,  which  means  it 
can  be  forecasted.  An  interval  type-2  FLS  is  designed  for 
sensed  signal  forecasting  and  is  compared  against  a  type-1 
FLS.  Simulation  results  show  that  the  interval  type-2  FLS 
performs  much  better  than  the  type-1  FLS  in  sensed  signal 
forecasting.  This  application  can  be  further  used  for  power 
on/off  control  in  wireless  sensors  to  save  battery  energy. 

Index  Terms — Wireless  sensors,  self-similarity,  forecasting, 
fuzzy  logic  system. 

I.  Introduction 

The  infusion  and  maturation  of  the  Micro  Mechanical 
System(MEMS),  computations,  and  wireless  communica¬ 
tion  technologies  has  advanced  the  development  of  wire¬ 
less  sensor  networks.  A  large  amount  of  low  cost  wireless 
sensor  nodes  can  be  densely  deployed  in  the  environment 
of  interest.  Usually,  such  kind  of  nodes  have  three  major 
functions:  sensing,  processing  and  communications.  The 
nodes  sense  the  environment  and  communicate  the  result 
to  the  user. 

In  this  paper,  we  use  Xbow  wireless  sensor  network 
professional  developer’s  kit  MOTE-Kit[3]  as  our  testbed 
to  get  several  data  sets,  from  different  scenarios.  First  of 
all,  we  show  that  the  sensed  signal  strength  is  self-similar 
and  long-range  dependent  using  variance-time  plotting  ,  a 
common  statistical  method  which  has  been  widely  used  to 
verify  self-similarity  of  time-series.  Since  the  sensed  signal 
strength  is  self-similar,  its  characteristics  can  be  captured. 
We  apply  a  type-1  FLS  and  an  interval  type-2  FLS  to 
do  sensed  signal  strength  forecasting.  This  study  can  be 
extended  to  power  on/off  control  in  wireless  sensors  to 
save  energy  consumption,  which  is  one  of  the  key  issues 
of  wireless  sensor  networks. 


In  the  following  sections,  Section  II  studied  the  self¬ 
similarity  of  sensed  signal  strength;  the  sensed  signal 
strength  forecasting  using  type-1  FLS  and  type-2  FLS  are 
presented  in  Section  III;  and  conclusions  and  future  works 
are  provided  in  Section  IV. 

II.  Self-Similarity  of  Sensor  Network  Data 

For  a  detailed  discussion  on  self-similarity  in  time- 
series,  see  [12]  [11].  Here  we  briefly  present  its  defini¬ 
tion  [1]. Given  a  zero-mean,  stationary  time-series  X  = 
(Xt;t  =  1,2,3,  •••),  we  define  the  ra-aggregated  series 

(m)  _  k  =  1, 2, 3,  •  •  *)  by  summing  the  original 

series  X  over  nonoverlapping  blocks  of  size  m.  Then  it’s 
said  that  X  is  H- self-similar,  if,  for  all  positive  m,  X 
has  the  same  distribution  as  X  rescaled  by  mH .  That  is, 

tm 

Xt  =  m~H  Xi  Vm&N  (1) 

i=(t— l)m+ 1 

If  X  is  H- self-similar,  it  has  the  same  autocorrelation  func¬ 
tion  r(k)  ~  E[(Xt  —  p)(Xt+k  —  m)]/^2  as  the  series  X 
for  all  m,  which  means  that  the  series  is  distributionally 
self-similar:  the  distribution  of  the  aggregated  series  is  the 
same  as  that  of  the  original. 

Self- similar  processes  can  show  long-range  dependence. 
A  process  with  long-range  dependence  has  an  autocorre¬ 
lation  function  r(k)  ~  as  k  ->  oo,  where  0  <  /?  <  1. 
The  degree  of  self-similarity  can  be  expressed  using  Hurst 
parameter  H  =  1  —  {3/2.  For  self-similar  series  with  long- 
range  dependence,  1/2  <  H  <  1.  As  H  1,  the  degree 
of  both  self-similarity  and  long-range  depence  increases. 

One  method  that  has  been  widely  used  to  verify  self¬ 
similarity  is  the  variance-time  plot ,  which  relies  on  the 
slowly  decaying  variance  of  a  self-similar  series.  The 
variance  of  X is  plotted  against  m  on  a  log-log  plot, 
and  a  straight  line  with  slope  (-/?)  greater  than  -1  is 
indicative  of  self- similarity,  and  the  parameter  H  is  given 


by  H  =  1  —  p/2.  We  use  this  method  to  verify  the  self¬ 
similarity  of  acoustic  signal. 

In  our  experiments,  8  sensors  were  deployed  in  a  lab. 
The  location  of  the  sensors  is  showed  in  Fig,  1.  We 
designed  two  scenarios,  one  is  with  a  fixed  source,  and 
the  other  is  without.  In  Fig.  2,  we  plot  the  variance  of 
X against  m  on  a  log-log  plot  for  8  sensor  nodes 
respectively  in  the  first  scenario  and  Fig.  3  is  under  the 
second  scenario.  From  the  two  figures,  it’s  very  clear  that 
the  no  matter  under  what  kind  of  condition  the  sensor 
network  data  have  self- similarity  because  their  traces  have 
slopes  much  greater  than  -1. 


Fig.  1.  The  deployment  of  the  eight  sensor  nodes  in  our  experiments. 


Fig.  2.  The  variance-time  plot  for  sensed  signal  strength  without  music 
as  background  during  3  hours.  The  sampling  rate  is  1024ms/sample. 

III.  Sensed  Signal  Strength  Forecasting  using 

INTERVAL  TYPE-2  FLS 

Since  the  sensed  signal  strength  is  self-similar,  its  char¬ 
acteristics  can  be  captured.  We  applied  a  type-1  and  an 


Fig.  3.  The  variance-time  plot  for  sensed  signal  strength  without  music 
as  background  during  3  hours.  The  sampling  rate  is  1024ms/sample. 

interval  type-2  FLS  to  sensed  signal  strength  forecast¬ 
ing.  Both  the  two  FLSs  have  16  rules  for  sensed  signal 
strength  modelling.  The  parameters  associated  with  each 
rule  are  determined  based  on  the  first  500  sensor  data, 
and  the  steepest  descent  algorithm  is  used  to  optimize  the 
parameters.  We  used  the  tuned  type-1  and  interval  type-2 
FLS  to  forecast  the  sensed  signal  strength  following  the 
training  data.  500  sensor  data  are  used  for  testing,  and  the 
forecasting  performances  are  evaluated  in  terms  of  root- 
mean  square  error. 

A.  Introduction  to  Type-2  Fuzzy  Sets 

The  concept  of  type-2  fuzzy  sets  was  introduced  by 
Zadeh  [13]  as  an  extension  of  the  concept  of  an  ordinary 
fuzzy  set,  i.e.,  a  type-1  fuzzy  set.  Type-2  fuzzy  sets  have 
grades  of  membership  that  are  themselves  fuzzy  [2].  A 
type-2  membership  grade  can  be  any  subset  in  [0, 1]  -  the 
primary  membership ;  and,  corresponding  to  each  primary 
membership,  there  is  a  secondary  membership  (which  can 
also  be  in  [0, 1])  that  defines  the  possibilities  for  the 
primary  membership.  A  type-1  fuzzy  set  is  a  special  case 
of  a  type- 2  fuzzy  set;  its  secondary  membership  function 
is  a  subset  with  only  one  element,  unity.  Type-2  fuzzy  sets 
allow  us  to  handle  linguistic  uncertainties,  as  typified  by 
the  adage  “words  can  mean  different  things  to  different 
people.”  A  fuzzy  relation  of  higher  type  (e.g.,  type-2)  has 
been  regarded  as  one  way  to  increase  the  fuzziness  of 
a  relation,  and,  according  to  Hisdal,  “increased  fuzziness 
in  a  description  means  increased  ability  to  handle  inexact 
information  in  a  logically  correct  manner  [5]”. 

Figure  4  shows  an  example  of  a  type-2  set.  The  domain 
of  the  membership  grade  corresponding  to  x  =  4  is  also 
shown.  The  membership  grade  for  every  point  is  a  Gaus¬ 
sian  type-1  set  contained  in  [0, 1],  such  a  set  was  called 


a  “Gaussian  type-2  set”.  When  the  membership  grade  for 
every  point  is  a  crisp  set,  the  domain  of  which  is  an  interval 
contained  in  [0,1],  Such  type-2  sets  were  called  “interval 
type-2  sets”  and  their  membership  grades  “interval  type-1 
sets”.  Interval  type-2  sets  are  very  useful  when  we  have 
no  other  knowledge  about  secondary  memberships. 


Fig.  4.  (a)  Pictorial  representation  of  a  Gaussian  type-2  set.  The 

secondary  memberships  in  this  type-1  fuzzy  set  are  shown  in  (b),  and 
are  Gaussian.  Note  that  this  set  is  called  a  Gaussian  type-2  set  because 
all  its  secondary  membership  functions  are  Gaussian.  The  “principal” 
membership  function  (the  bold  line),  which  is  triangular  in  this  case,  can 
be  of  any  shape. 


B.  Introduction  to  Type-2  Fuzzy  Logic  Systems:  An 
Overview 

Figure  5  shows  the  structure  of  a  type-2  FLS.  It  is 
very  similar  to  the  structure  of  a  type-1  FLS  [9].  For  a 
type-1  FLS,  the  output  processing  block  only  contains  the 
defuzzifier.  We  assume  that  the  reader  is  familiar  with  type- 
1  FLSs,  so  that  here  we  focus  only  on  the  similarities  and 
differences  between  the  two  FLSs. 

The  fuzzifier  maps  the  crisp  input  into  a  fuzzy  set.  This 
fuzzy  set  can,  in  general,  be  a  type-2  set. 

In  the  type-1  case,  we  generally  have  “IF-THEN”  rules, 
where  the  Zth  rule  has  the  form  “ Rl  :  IF  xx  is  F*x  and 
X2  is  F2  and  •  •  •  and  xp  is  Flp,  THEN  y  is  G*”,  where: 
XiS  are  inputs;  F^s  are  antecedent  sets  (i  =  1  y 

is  the  output;  and  G*s  are  consequent  sets.  The  distinction 


TYPE-2  FUZZY  LOGIC  SYSTEM 


Fig.  5.  The  structure  of  a  type-2  FLS.  In  order  to  emphasize  the 
importance  of  the  type-reduced  set,  we  have  shown  two  outputs  for  the 
type-2  FLS,  the  type-reduced  set  and  the  crisp  defuzzified  value. 

between  type- 1  and  type-2  is  associated  with  the  nature  of 
the  membership  functions,  which  is  not  important  while 
forming  rules;  hence,  the  structure  of  the  rules  remains 
exactly  the  same  in  the  type-2  case,  the  only  difference 
being  that  now  some  or  all  of  the  sets  involved  are  of 
type-2;  so,  the  Ith  rule  in  a  type-2  FLS  has  the  form  “ Rl  : 

~i  ~ 1 

IF  xi  is  F2  and  x2  is  F2  and  •  •  •  and  xp  is  Fp,  THEN  y  is 

g'”. 

In  the  type-2  case,  the  inference  process  is  very  similar 
to  that  in  type- 1 .  The  inference  engine  combines  rules  and 
gives  a  mapping  from  input  type-2  fuzzy  sets  to  output 
type-2  fuzzy  sets.  To  do  this,  one  needs  to  find  unions 
and  intersections  of  type-2  sets,  as  well  as  compositions 
of  type-2  relations. 

In  a  type-1  FLS,  the  defuzzifier  produces  a  crisp  output 
from  the  fuzzy  set  that  is  the  output  of  the  inference  engine, 
i.e.,  a  type-0  (crisp)  output  is  obtained  from  a  type-1  set. 
In  the  type-2  case,  the  output  of  the  inference  engine 
is  a  type-2  set;  so,  “extended  versions”  (using  Zadeh’s 
Extension  Principle  [13])  of  type-1  defuzzification  methods 
were  used.  This  extended  defuzzification  gives  a  type-1 
fuzzy  set.  Since  this  operation  takes  us  from  the  type- 
2  output  sets  of  the  FLS  to  a  type-1  set,  this  operation 
was  called  “type-reduction”  and  the  type-reduced  set  so 
obtained  a  “type-reduced  set”. 

To  obtain  a  crisp  output  from  a  type-2  FLS,  we  can 
defuzzify  the  type-reduced  set.  The  most  natural  way  of 
doing  this  seems  to  be  by  finding  the  centroid  of  the  type- 
reduced  set;  however,  there  exist  other  possibilities  like 
choosing  the  highest  membership  point  in  the  type-reduced 
set. 

General  type-2  FLSs  are  computationally  intensive,  be¬ 
cause  type-reduction  is  very  intensive.  Things  simplify  a  lot 
when  secondary  membership  functions  (MFs)  are  interval 
sets  (in  this  case,  the  secondary  memberships  are  either  0 
or  1).  When  the  secondary  MFs  are  interval  sets,  the  type-2 
FLSs  were  called  “interval  type-2  FLSs”.  In  [7],  Liang  and 
Mendel  proposed  the  theory  and  design  of  interval  type- 
2  fuzzy  logic  systems  (FLSs).  They  proposed  an  efficient 


and  simplified  method  to  compute  the  input  and  antecedent 
operations  for  interval  type-2  FLSs,  one  that  is  based  on 
a  general  inference  formula  for  them.  They  introduced  the 
concept  of  upper  and  lower  membership  functions  (MFs) 
and  illustrated  the  efficient  inference  method  for  the  case 
of  Gaussian  primary  MFs.  They  also  proposed  a  method 
for  designing  an  interval  type-2  FLS  in  which  we  tune  its 
parameters. 


C.  Why  Type-2  FLS  is  Necessary? 

FLSs  have  been  extensively  used  in  time-series  forecast¬ 
ing  (e.g.,  [4],  [10],  [7]).  Here  we  designed  a  sensed  signal 
strength  forecasting  scheme  for  wireless  sensors  using 
interval  type-2  FLS.  Why  choose  type-2  FLS?  Acoustic 
amplitude  sensor  node  measures  sound  amplitude  at  the 
microphone.  Assuming  that  the  sound  source  is  a  point 
source  and  sound  propagation  is  lossless  and  isotropic, 
a  root-mean-squared  (RMS)  amplitude  measurement  z  is 
related  to  the  sound  source  position  X  as 


_  a 


+  w, 


(2) 


where  a  is  the  RMS  amplitude  of  the  sound  source,  ^  is 
the  location  of  the  sensor,  and  w  is  RMS  measurement 
noise  [6] .  w  is  modelled  as  a  Gaussian  with  zero  mean  and 
variance  a2.  The  sound  source  amplitude  a  is  also  mod¬ 
elled  as  a  random  quantity,  which  is  uniformly  distributed 
in  the  interval  [a^a/^].  Given  the  location  of  the  sound 
source  X  and  the  sensor  position  c,  p^|f  is  uniformly 
distributed  as  a  is.  Therefore,  z  should  be  modelled  as 
a  Gaussian  primary  MF  having  a  fixed  standard  devia¬ 
tion,  aj.,  and  an  uncertain  mean  that  takes  on  values  in 


lA(xk)  =  exp 
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where:  k  =  1, . . .  ,p;  p  is  the  number  of  antecedents;  l  = 
1 , . . . ,  M;  and,  M  is  the  number  of  rules.  It  is  shown  in 
Fig  6. 


D.  Simulations 

Our  simulations  were  based  on  N  =  1000  samples, 
x(l),  x(2),  2(1000).  The  first  500  data,  2(1),  2(2), 

. . 2(500),  are  for  training,  and  the  remaining  500  data, 
2(501),  2(502), . . .  ,2(1000)  are  for  testing.  In  Fig.  7,  we 
plot  the  sensed  data  that  we  used  for  training  and  testing, 
2(1),  2(2),  ...,  2(1000). 

We  used  four  antecedents  for  forecasting,  i.e.,  2 (k  —  3), 
x(k- 2),  x{k  1),  and  2 (k)  were  used  to  predict  2 (fc  +  1). 
For  type-1  FLS,  the  rules  are  designed  such  as: 

Rl  :  IF  x(k  -  3)  is  F^  and  2 (k  -  2)  is  Fz2  and  2 (k  -  1)  is 
F*i  and  2 (k)  is  F2,  THEN  x(k  +  1)  is  G*. 


1 

0.9 

0.8 

0.7 

0.6 

0.5 

0.4 

0.3 

0.2 

0.1 


Fig.  6.  The  interval  type-2  MFs  with  fixed  std  and  uncertain  mean. 
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Fig.  7.  Sensed  data  for  1024  seconds,  the  sample  rate  is  1024ms/sample, 
x(l),  x{2) . a;(1000). 


For  interval  type-2  FLS,  we  design  the  rules  as: 

Rl  :  IF  2 (k  —  3)  is  F*  and  x(k  —  2)  is  F2  and  x(k  —  1)  is 
Fj  and  x(k)  is  Fg,  THEN  x(k  +  1)  is  &. 

In  this  paper,  we  use  height  defuzzifier  [9]  for  type-1  FLS, 
the  height  of  G*  is  y1.  For  interval  type-2  FLS,  we  use 
Center-of-Sets  type  reduction  [7]. 

As  in  [7],  we  used  only  2  fuzzy  sets  for  each  antecedent, 
so  the  number  of  rules  is  24  =  16.  For  type-1  FLS, 
Gaussian  membership  functions  (MFs)  were  chosen  for 
the  antecedents;  for  interval  type-2  FLS,  Gaussian  primary 
MFs  with  fixed  std  and  uncertain  mean  were  chosen  for 
the  antecedents.  The  initial  locations  of  antecedent  MFs 
were  based  on  the  mean,  mu  and  std,  <ju  of  the  first  500 
points,  2(1),  2(2),  . . .,  2(500). 

We  used  steepest  descent  algorithm  to  train  all  the 
parameters  based  on  the  training  data.  After  training,  the 


rules  were  fixed,  and  we  tested  the  FL  forecaster  based 
on  the  remaining  500  noisy  points,  rr(501),  a(502),  . . 
*(1000). 

We  compared  the  performance  of  the  interval  type-2  FLS 
with  the  type-1  FLS  for  sensed  signal  strength  forecasting. 
For  each  of  the  2  above  methods,  we  ran  100  Monte-Carlo 
realizations  to  eliminate  the  randomness  of  the  parameters, 
and  the  two  FLSs  was  tuned  using  a  simple  steepest- 
descent  algorithm.  We  used  the  testing  data  to  see  how 
each  FLS  performed  by  evaluating  the  root-mean-square- 
error  (RMSE)  between  the  defuzzified  output  of  the  FLS 
and  the  actual  sensor  data  ( x(k  +  1)),  i.e., 


RMSE: 
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where  xfe  =  [x(k  -  3 ),x(k  -  2 ),x(k  -  1),  x(k)]T,  and  T 
denotes  transpose. 

The  RMSE  values  for  each  design,  we  summarize  the 
mean  of  the  RMSEs  and  shown  in  Table  I. 
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SIGNAL  STRENGTH  FORECASTING. 
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System 

RMSE 

type-1  FLS 

9.41 

interval  type-2  FLS 

3.85 

The  results  show  that  the  interval  type-2  FLS  for  sensed 
signal  strength  forecasting  is  much  better  than  using  type-1 
FLS. 


IV.  Conclusions  and  Future  Works 

In  this  paper,  we  studied  the  sensed  signal  strength  using 
the  data  collected  in  MOTE-Kit[3]  testbed  and  show  that 
the  sensor  data  are  self-similar,  which  validate  that  sensed 
signal  strength  is  forecastable  because  self-similar  time- 
series  can  be  forecasted.  Based  on  the  analysis  of  the 
sensor  model,  we  applied  an  interval  type-2  FLS  to  sensed 
signal  strength  forecasting,  and  simulation  results  show 
that  it  performs  much  better  than  does  a  type-1  FLS.  All 
these  studies  are  very  important  for  power  on/off  control 
in  wireless  sensors  to  save  battery  energy,  which  are  the 
future  works  that  we  will  investigate. 
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Abstract—  The  performance  evaluation  is  one  the  most  impor¬ 
tant  research  topics  for  the  Wireless  Sensor  Networks  (WSN). 
Delay  and  energy  efficiency  are  two  important  parameters  to 
evaluate  the  WSN  quality.  In  the  WSN,  the  interference  will 
affect  the  packets  transmission.  When  a  sensor  need  to  send  a 
packet,  we  choose  the  parameter  SIR  as  the  threshold  to  decide 
whether  send  or  not.  If  the  SIR  threshold  is  high,  the  probability 
of  the  sending  a  packet  will  decrease  and  the  packet  need  to  be 
kept  in  the  queue,  and  the  delay  of  the  packet  will  increase.  If  the 
SIR  threshold  is  lower,  the  delay  is  decrease.  However,  in  order 
to  overcome  the  interference,  the  energy  cost  for  each  packet 
will  increase.  Simulation  shows  the  SIR  threshold  can  control 
the  delay/energy  performance  in  WSN. 

I.  Introduction 

Thanks  to  the  rapid  development  in  low  power  wireless 
communication,  microprocessor  hardware  in  conjunction  with 
the  significant  process  in  distributed  signal  process,  wireless 
sensor  networks  (WSNs)  approach  to  a  new  technological 
vision.  While  a  lot  of  research  has  been  concentrated  on  the 
some  important  aspects  of  WSNs,  such  as  energy  efficiency, 
protocol  design  and  network  deployment,  the  performance 
evaluation  in  WSNs  is  rarely  studied  [1].  As  the  WSNs  is 
widely  deployed  in  military  and  commercial  application  [2], 
energy  efficiency  and  delay-aware  become  more  important  for 
WSNs. 

Interference  is  the  major  limiting  factor  in  the  performance 
of  wireless  sensor  networks.  There  are  several  kinds  of  sources 
of  interference.  In  this  paper,  we  only  consider  the  interference 
caused  by  neighboring  sensors.  When  the  neighboring  sensors 
send  out  packets,  they  will  affect  the  packet  received  by 
the  destination  sensor.  In  Figure  1,  the  source  sensor  is 
sending  packets  to  the  destination  sensor,  sensors  B,  C,  E 
are  also  sending  packets  to  their  destination  sensor.  Becase 
every  sensor  has  a  non-directinal  antenna,  it  will  broadcast 
its  packets  to  all  the  sensors  around  it.  Sensors  A,  D  are  not 
sending  packets  at  this  time,  they  will  not  cause  interference  to 
the  destination  sensor.  We  calculate  the  interference  in  every 
small  time  period,  so  we  can  know  the  interference  stays 
constant  during  each  time  period,  but  fluctuates  randomly  in 
different  periods. 
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Fig.  1.  the  interference  caused  by  neighboring  sensors 

The  main  dilemma  that  the  transmitter  faces  is  when  the 
interference  is  very  large,  the  packets  the  destination  received 
will  affected  by  the  interference.  So  we  can  control  the 
packet  transmission  as  the  following:  when  the  transmitter 
observers  high  interference  in  the  channel,  that  is  the  signal-to- 
interference  ratio  (SIR)  less  than  a  threshold,  it  would  be  better 
to  back  off,  buffer  the  traffic  and  wait  for  the  interference  to 
subside  before  it  transmits.  As  it  has  backed  off,  the  buffer  is 
filling  up  with  new  packet  arrivals  and  delay  rises. 

In  order  to  evaluate  the  algorithm,  we  need  to  understand 
the  major  parameters  about  the  wireless  sensor  network. 

A.  Delay 

Because  data  communications  in  the  sensor  networks  has 
trimming  constraints,  it  is  important  to  design  the  network 
algorithm  to  meet  a  kind  of  end-end  deadlines  [3]. 

B.  System  lifetime 

It  is  not  convenient  to  recharge  the  sensor  node  battery, 
so  the  energy  efficiency  is  extremely  important  for  sensor 
network.  And  the  network  should  keep  an  enough  number  of 
“live”  sensor  nodes  to  collect  data,  which  means  the  network 
need  to  keep  the  energy  among  the  sensor  nodes  in  balance. 
We  use  the  remaining  alive  sensor  nodes  as  the  parameter  of 
the  networks  lifetime. 


C.  Network  efficiency 

The  sensor  network  is  used  to  collect  data  and  transfer 
packets.  The  amount  of  packets  transmitted  is  one  of  the 
parameters  to  evaluate  the  networks  efficiency. 

The  rest  of  this  paper  is  organized  as  follows.  In  section  II, 
we  introduce  some  preliminaries.  In  section  III,  we  introduce 
the  packet  transmission  scheduling.  In  section  IV,  we  introduce 
the  performance  analysis.  In  section  V,  we  introduce  the 
simulation  result.  In  section  VI,  we  conclude  the  paper. 

II.  Preliminaries 


C.  Energy 

A  sensor  node  consumes  significant  energy  when  it  trans¬ 
mits  or  receives  a  packet.  But  we  will  not  consider  the  energy 
consumed  when  the  sensor  node  is  idle. 

The  distance  between  two  nodes  are  variable  in  the  WSN 
and  the  power  loss  model  is  used.  To  send  the  packet,  the 
sender  consumes  [6], 

Ptx  =  Pelec  +  £fs  *  d2  (4) 

and  to  receive  the  packet,  the  receiver  consumes, 

Prx  —  P gfec  (5) 


A.  Interference 

Frequency  reuse  implies  that  in  a  given  coverage  area 
there  is  several  wireless  sensors  that  use  the  same  set  of 
frequencies.  These  sensors  are  called  co-channel  sensor,  and 
the  interference  between  signals  from  these  sensors  is  called 
co-channel  interference  [4]. 

Let  N  be  the  number  of  co-channel  interfering  sensors, 
then  the  signal-to-interference  ratio(SIR)  for  a  mobile  receiver 
which  monitor  a  forwarding  channel  can  be  expressed  as  : 
ratio(SIR)  for  a  mobile  receiver  which  monitor  a  forwarding 
channel  can  be  expressed  as: 


S^P^ 
I  Pi 
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Where  Ps  is  the  desired  signal  power  from  the  desired 
sensor  and  P/.  is  the  interference  power  caused  by  the  ith 
interfering  co-channel  sensor. 


B.  Queueing  Model 

We  setup  a  threshold  to  decide  whether  the  source  sensor 
send  out  a  packet.  The  sensor  is  keeping  collecting  data  with  a 
Poisson  distribution,  while  the  service  is  general  distribution. 

The  M/G/l  queue  has  exponentially  distributed  interarrival 
times,  general  process  times,  and  a  single  server.  The  inclusion 
of  general  process  times  makes  this  formula  more  relevant  for 
wafer  fabs,  where  processing  times  are  typically  not  as  variable 
as  the  exponential  distribution  used  with  M/M/1  queues. 

Let  A  be  the  arrival  rate,  ji  be  the  service  rate  and  the 
P  —  ~  be  the  traffic,  if  we  substitute  in  the  variance  plus 
the  squared  mean  for  the  second  moment  of  the  service  time 
distribution  [5],  we  can  get: 

Average  queuing  length: 
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where  Peiec  represents  the  power  that  is  necessary  for 
digital  processing,  modulation,  and  ejs  represents  the  power 
dissipated  in  the  amplifier  for  the  free  space  distance  d 
transmission. 

D.  Delay 

The  packet  transmission  latency  between  the  sensors  in¬ 
cludes  three  parts:  the  wireless  channel  transmission  delay, 
the  Physical/Mac  layer  transmission  delay,  and  the  queuing 
delay  [7]. 

Defined  D  as  the  distance  between  two  sensors  and  C  as 
the  light  speed,  the  wireless  channel  transmission  delay  as: 

Delay  ch  =  ^  (6) 

The  Physical/Mac  layer  transmission  delay  will  be  decided 
by  interaction  of  the  transmitter  and  the  receive  channel,  the 
sensor  density  and  the  sensor  traffic  intensity  etc. 

The  queuing  delay  is  decided  by  the  sensor  I/O  system 
processing  rate,  the  subqueue  length  in  the  sensor. 

In  order  to  make  the  system“stable”,  the  rate  at  which  sensor 
node  transfers  packets  intended  for  its  destination  must  satisfy 
all  sensor  that  the  queuing  lengths  will  not  be  infinite  and  the 
average  delays  will  be  bounded. 

III.  Packet  Transmission  Scheduling 

Time  division  multiple  access(TDMA)system  divided  the 
radio  spectrum  into  time  slots,  and  in  each  slot  only  one 
user  is  allowed  to  either  transmit  or  receive.  The  TDMA 
frame(Figure.2)  has  two  segments:  the  header  segment  consists 
of  request  slots  and  the  trailer  segment  consists  of  information 
slots.  Request  slots  are  minislots  and  an  information  slot  is 
much  longer  in  order  to  transfer  data  packets. 


Average  queuing  time: 

_Lg_  A(Var[50]+p-2) 

9  A  2(1  —  p) 


Fig.  2.  TDMA  frame 

(3)  We  setup  a  threshold  to  decide  whether  the  source  sensor 
send  out  a  packet  or  insert  the  packet  in  its  queue.  The  sensor 


is  keeping  collecting  data  with  a  Poisson  distribution.  When 
one  packet  is  generated,  insert  the  packet  into  the  tail  of  the 
queue  of  the  sensor,  then  compare  the  SIR  of  the  destination 
sensor.  If  the  SIR  large  than  the  threshold,  the  source  sensor 
sends  out  its  packet  which  locate  at  the  head  of  the  queue  to 
the  destination  sensor,  else  keep  idle. 

A.  SIR  calculation 

Every  sensor  calculates  its  SIR  continuously.  It  needs  to 
know  the  signal  power  S  and  interference  power  7*  caused 
by  the  ith.  interfering  sensor.  The  transmission  energy  each 
packets  cost  can  be  decided  by  (4).  It  is  a  function  of  the 
distance  between  the  source  sensor  and  the  destination  sensor. 

D(d>  s )  =  VWd  -zs)2  +  (yd  -  Vs)2  (7) 


Fig.  3.  one-step  Markov  path  model 


Where  x&,  ya ,  xs,  ys  are  X  position,  Y  position  of  the 
destination  sensor  and  the  source  sensor. 


Ph  —  PSi  ~  £fs 


D&k)2 
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D(d,s)2 

Where  PSi  is  the  power  that  the  interference  sensor  used  to 
send  its  own  packets.  Pji  is  the  interference  power  generated 
by  sensor  I,  efs  represents  the  power  dissipated  in  the  amplifier 
for  the  free  space  distance  transmission. 


D2  :  denote  the  single  TDMA  frame  that  the  data  packets 
has  to  wait  in  the  (z+l)th  TDMA  frame. 

£>3  :  denote  the  request  slot  subframe  that  the  data  packet 
has  to  wait  in  the  (i+2)th  TDMA  frame. 

I>4  :  denote  the  time  interval  that  the  data  packets  have  to 
wait  for  other  data  packets  in  the  ( i+2 )  frame. 

The  data  packets  queueing  delay: 

Td  =  Di  +  T>2  +  D3  +  D4  (10) 
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Where  Peiec  represents  the  power  that  is  necessary  for 
digital  processing.  According  to  (5),  we  konw  is  the  power 
the  receiver  comsumes  is  Peiec  . 

According  to  the  SIR  calculation  algorithm,  the  sensors  need 
to  exchange  their  information  when  they  are  sending  packets. 
In  the  head  of  every  packet,  it  will  include  the  X  position  xs, 
Y  position  ys  and  the  signal  power  Ps. 


B.  one-step  Markov  path  model 

The  sensors  are  roaming  independently  with  variable  ground 
speed.  The  mobility  model  is  called  one-step  Markov  path 
model  [8].  The  probability  of  moving  in  the  same  direction 
as  the  previous  move  is  higher  than  other  directions  in  this 
model,  which  means  this  model  has  memory.  Fig.3  shows  the 
probability  of  the  six  directions. 


IV.  Performance  Analysis 

The  delay  of  this  WSN  system  can  be  analyzed  using 
the  M/D/1  queueing  theory.  However,  we  cannot  achieve  the 
variance  of  the  general  distribution.  We  consider  approximate 
analysis  of  ON/OFF  data  traffic  under  light  load  assumptions 
that  we  believe  most  networks  will  operate  in  [9]. 

In  the  following,  we  assume  the  light  load  condition.  To 
calculate  the  mean  data  packet  queuing  delay,  let  us  first  define 
some  variables. 

Di  :  denote  the  time  interval  from  the  time  the  data  packet 
arrives  in  the  ith  TDMA  frame  to  the  beginning  of  the  next 
TDMA  frame. 


Let  us  assume  the  D\  is  uniformly  distribution  between  the 
beginnings  of  two  consecutive  TDMA  frames,  e.g.,  between 
zero  and  Ttdma,  where  Ttdma  is  the  TDMA  frame  interval. 
Furthermore,  let  us  assume  that  there  are  rrid  integer-multiple 
of  one  packet  of  size  equivalent  to  one  information  slot  packets 
generated  at  a  constant  interval  of  Tt dm A/md  during  the  ON 
period  in  a  TDMA  frame  interval,  where  md  is  the  number 
of  information  slots  used  by  a  data  source  in  a  TDMA  frame 
interval  during  its  ON  period.  Its  pdf  is  given  by: 


f  U)  -  /  */  i^-Ttdma  <  t  <  £ Ttdma,  a  =  1 

JDlK  ’  \  0,  othewise 

(11) 

Where  a  is  the  arrival  order  of  the  packets  in  an  TDMA 
frame  time.  Thus  the  mean  of  D\  is: 


E[Di\  -  Ttdma  (12) 

The  pdf  of  D2  is  give  by: 


f  1,*/  t  =  TDMA 
\  0,  otherwide 


while  the  pdf  of  7)3  is  given  by 
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M/  t  =  NRTR 
0,  otherwise 
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Where  NR  is  the  number  of  request  slots  in  the  request  slot 
subframe  and  TR  is  the  time  duration  of  one  request  slot.  That 
is,  the  means  of  D2  and  D3  are  respectively  given  by: 


E[D2]  —  TTdma ,  (15) 


and 


12.5 


E[DS]  =  NrTr  (16) 


Let  Nd  be  the  number  of  data  users  in  the  active  (ON)  state 
that  constitutes  a  self-similar  traffic.  The  pdf  of  Nd  is  given 
by: 


Pr[Nd  =  nd\  «  (  )  Pr:tjl  -  Pon,d)M*~n*  (17) 

Where  Pon,d  ~  *i  is  the  mean  ON  period  for  a 

ON/OFF  data  traffic  and  1 2  is  the  mean  OFF  period  for  the 
ON/OFF  data  traffic.  The  conditional  probability  of  Nd  is 
given  by: 


Fig.  4.  Average  Delay  performance  of  the  three  thresholds 


Pr[Nd  =  nd\Nd  >  0]  =  Pr[Nd  =  nd}/(  1  -  Pr[Nd  =  0]) 

(18) 

Let  S  be  the  total  number  of  data  packets  generated  in  a 
TDMA  frame  interval.  Its  pdf  is  given  by: 

Pr[S  =  s  =  mdNd ]  =  Pr[Nd  =  nd  =  S/md\  (19) 

Each  of  k  x  md  data  packets  is  assumed  to  have  equal 
probability  of  transmitting  in  each  of  the  k  x  md  consecu¬ 
tive  information  slots.  The  pdf  of  Z)4  given  a,  where  a=l,2 
,3 ,...md,  is  given  by: 

/d4 \a(t  =  ((z  -  1  )md  +  a  -  1  )Tj) 

Md 

=  ^  E  hPr\Nd  =  >  0]$[i-  ((z  -  1)7714  -Fa  -  1)T/] 

i  =  1, 2,  ■  •  *  ,  Md 

(20) 

Where  £(•)  is  an  impulse  function.  Its  mean  is  given  by 


md  Md  Md. 

£N=EE^E 


0=1  t=l  k—i 


ipr[Nd  =  k\Nd  >  0]((z  -  1  )md  +  a  -  1)T> 


(21) 


V.  Simulation 


We  implemented  the  simulation  model  using  the  OPNET 
modeler.  The  simulation  region  is  800x800  meters.  In  order 
to  overcome  the  interference,  the  sensor  would  cost  different 
energy  with  different  SIR  threshold.  We  chose  the  SIR  thresh¬ 
olds  as  three  cases:  high(0.5),  medium(0.3)  and  low(O.l). 

There  were  12  sensor  nodes  in  the  simulation  model,  and 
the  sensors  were  roaming  independently  with  variable  ground 
speed  between  1  to  9  meters  per  second.  The  mobility  model 
was  called  one-step  Markov  path  model.  The  movement  would 
change  the  distance  between  sensor  and  it  would  also  change 
the  interfemce  power. 


1)  Average  Delay:  We  used  the  average  delay  parameter 
[see  (22)]  to  evaluate  the  network  performance.  Each  packet 
was  labeled  a  timestamp  when  it  was  generated  by  the  source 
sensor  node.  When  its  destination  sensor  node  received  it,  the 
time  interval  was  the  transmission  delay. 


Average  Delay  = 


(22) 


Observing  Fig.4  the  delay  performance  of  three  thresholds, 
the  lower  threshold,  the  better  delay  performance.  When  we 
chose  the  low  SIR  threshold,  the  sensor  had  a  high  probability 
to  send  out  a  packet  than  high  and  medium  cases.  The  average 
delay  of  the  low  threshod  around  21%  was  shorter  than  that 
of  the  high  case.  However,  the  average  delays  of  high  and 
medium  case  were  almost  the  same.  The  reason  was  that 
in  most  cases,  the  interference  was  between  0.1  and  0.3.  If 
the  service  is  time-sensitive,  such  as  video  or  audio  service, 
we  would  choose  the  low  SIR  threshold  to  meet  their  delay 
performance  requirement. 

2)  Energy  Efficiency:  When  we  chose  the  lower  SIR  thresh¬ 
old,  that  meant  we  need  to  increase  the  transmission  power 
to  overcome  the  interference.  According  to  (4)  and  (5),  for 
high  SIR  case,  the  sensor  node  consumed  Ptx  watts  during 
transmissions  and  Prx  watts  during  reception.  We  assumed 
Pelec  was  equal  to  6.0xl0“4  and  efs  was  equal  to  6.0xl0“4. 
For  medium  SIR  case,  both  the  transmission  and  reception 
energy  comsumed  would  be  1.67  times  of  those  of  the  high 
SIR  case.  And  for  low  SIR  threshold  case,  both  the  trans¬ 
mission  and  reception  energy  comsumed  would  be  5  times  of 
those  of  the  high  SIR  case.  The  lower  SIR  threshold,  more 
energy  comsumed.  In  the  wireless  sensor  network,  we  used 
two  parameters:  the  number  of  sensor  nodes  alive  and  the 
remaining  energy  to  describe  the  energy  efficiency. 

When  the  remaining  energy  of  a  sensor  node  was  lower 
than  a  certain  threshold,  the  sensor  was  considered  as  “dead”. 
In  this  simualtion,  we  chose  1.2xl0-3  as  the  threshold.  A 
sensor  was  “dead”  meant  it  could  not  transmit/receive  packets 
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Fig.  5.  Node  alive  performance  of  the  three  thresholds 

any  longer,  so  it  would  be  ignored  by  the  sensor  network.  The 
sensor  was  used  to  collect  data  and  transmit  the  packets.  The 
number  of  sensor  of  a  wireless  sensor  networks  which  was 
below  a  certain  threshold  means  this  network  does  not  work. 

As  Fig.5  showed,  the  remaining  sensor  nodes  alive  of  the 
lowe  SIR  threshold  was  much  worse  that  of  the  high  case  and 
the  medium  case  was  between  the  low  case  and  high  case.  The 
reason  was,  the  higher  SIR  threshold,  less  energy  comsumed. 
The  time  of  the  first  sensor  node  “dead”  of  the  high  SIR 
threshold  was  over  two  time  longer  of  that  of  the  low  SIR 
threshold  case. 

Fig.6  showed  the  remaining  energy  of  the  three  thresholds. 
We  assumed  that  the  energy  of  each  sensor  was  0.1J  and  the 
packet  size  was  600  bit,  and  the  channel  transmission  rate  was 
1M  bps.  So  when  the  sensor  transmitted  or  received  a  packet, 
it  would  cost  6  x  10-4  second.  If  a  sensor  node  transmitted 
Nums  packets  (each  packet  cost  6  x  10~4  second)  and  receives 
Numr  packets  (each  packets  also  cost  6  x  10-4  second)  and 
it  was  roaming  in  the  network  for  Tm,  we  could  get  the 
remaining  energy  Ei  of  this  sensor  node: 

Ei  -  0.1  —  (Ptx  x 6 x  10-4 x Nums-\-Prx  x 6 x  10”4 x Numr) 

(23) 

The  remaining  energy  Ew  of  the  whole  networks  was 
described  as: 

12 

Ew  =  (24) 

i=  1 

Figure.7  showed  the  remaining  energy  of  the  high  SIR 
threshold  was  not  dropping  as  sharply  as  that  of  the  the  low 
SIR  threshold.  The  medium  was  between  them. 

3)  Network  quality:  The  role  of  the  wireless  sensor  network 
in  the  real  world  was  to  collect  data  and  transmit  packets.  In 
our  simulation,  we  assumed  the  collecting  data  distribution 
of  the  sensor  node  was  Poisson  distribution  and  the  arriving 
interval  was  1  second.  Observing  from  Fig. 7,  the  higher  SIR 
threshold,  more  packets  collected.  It  was  because  the  energy 
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Fig.  6.  Remaining  energy  of  the  three  thresholds 


Fig.  7.  Packets  received  of  the  three  threshold 


cost  by  each  packet  in  the  high  case  is  less  that  of  the  low 
case. 

VI.  Conclusion 

In  the  wireless  sensor  network,  the  interference  will  greatly 
affect  the  packets  transmission.  If  the  SIR  is  high,  the  proba¬ 
bility  of  the  sending  a  packet  will  decrease  and  the  packet  need 
to  be  kept  in  the  queue.  The  delay  of  the  packet  will  increase. 
If  the  SIR  is  lower,  the  delay  is  decrease.  However,  in  order  to 
overcome  the  interference,  the  energy  cost  for  each  packet  will 
increase.  We  realize  the  wireless  commun cation  system  with 
TDMA  and  we  analyze  the  performance  of  delay  and  energy 
according  to  different  SIR  threshold.  Simulation  shows  the  SIR 
threshold  can  control  the  delay/energy  performance  in  WSN. 
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Abstract—  Wireless  Sensor  Networks  (WSNs)  pose  consid¬ 
erable  technical  challenges  in  energy  efficiency  for  resource 
constrained  (energy,  bandwidth,  computation  and  memory).  In 
this  paper,  we  propose  an  asynchronous  contention-based  and 
energy-efficient  MAC  protocol  (ACEMAC).  Our  algorithm  ame¬ 
liorates  the  biggest  limitation  for  existed  contention-based  MAC 
protocols,  time  synchronization  needed.  We  propose  a  designing 
model  to  optimum  power  off  duration,  power  on  duration  and 
system  stage  construct.  We  also  make  theory  analysis  on  average 
latency  (Point-to-Point)  for  our  algorithm.  Simulation  results 
verify  that  our  algorithm  can  successfully  obtain  the  optimum  re¬ 
schedule  duty-cycle,  power  off  duration  and  power  on  duration 
to  achieve  higher  average  successful  transmission  rate,  shorter 
data  packet  average  latency  and  longer  network  lifetime. 

I.  Introduction 

Wireless  sensor  networking  is  an  emerging  technology 
that  has  a  wide  range  of  potential  applications  including: 
military  sensing,  physical  security,  air  traffic  control,  traffic 
surveillance,  video  surveillance,  industrial  and  manufacturing 
automation,  distributed  robotics,  environment  monitoring,  and 
structures  monitoring [15].  However,  individual  sensor  nodes 
are  resource  constrained  (energy,  bandwidth,  computation, 
memory).  So  WSNs  pose  considerable  technical  challenges 
in  energy  efficiency,  spectral  efficiency  and  security.  In  this 
paper,  we  focus  on  the  energy  efficiency  problem  for  WSNs. 

There  are  two  approaches  for  energy  reservation  on  commu¬ 
nication.  One  approach  does  power  control  at  physical  layer. 
The  other  approach  implements  energy  reservation  task  on 
MAC  layer.  Considerable  energy  in  traditional  MAC  protocol 
is  used  by  idle  listening,  i.e.,  listening  to  receive  possible  traffic 
that  is  not  sent.  For  example,  the  digital  2  Mbps  Wireless 
LAN  module  (IEEE  802.1  l/2Mbps[  10])  specification  shows 
the  energy  consuming  ratio  of  idle:  receive:  send  isl  :2:2.5[14]. 
So  about  one  fifth  of  total  energy,  which  should  be  saved,  is 
consumed.  It  is  obvious  that  enforcing  some  sensor  nodes, 
which  are  staying  at  idle  state,  to  turn  off  their  RF  interfaces 
is  an  effective  method  to  implement  energy  reservation 

In  literature,  some  energy-efficient  MAC  protocols  for 
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WSNs  are  proposed.  The  first  energy-efficient  MAC  protocol 
for  ad  hoc  networks:  PAMAS[5]  is  proposed  in  1999.  PAMAS 
conserves  battery  power  by  intelligently  powering  off  users 
that  are  not  actively  transmitting  or  receiving  packets.  But 
the  disadvantage  for  this  protocol  is,  two  different  physical 
channels  are  needed,  i.e.,  control  channel  and  traffic  chan¬ 
nel.  After  PAMAS,  some  solutions  have  been  proposed  for 
WSNs.  They  can  be  classified  into  two  categories:  Schedule- 
based  and  Contention-based.  TRAMA[2],  a  schedule-based 
MAC  protocol,  employs  a  traffic  adaptive  and  distributed 
election  scheme  to  allocate  the  system  time  for  different 
sensor  nodes.  According  to  the  election  scheme,  receivers  are 
selected  basing  on  schedules  announced  by  transmitters.  For 
S-MAC[3],  a  contention-based  MAC  protocol,  borrows  the 
main  idea  from  PAMAS  to  implement  energy  reservation.  But 
S-MAC  has  only  one  common  channel  to  transmit  data  and 
control  information.  Considering  the  fixed  duty-cycle  is  not 
optimal,  S-MAC ’s  improvement  version  T-MAC  is  proposed 
in  [4].  For  S-MAC  and  T-MAC,  they  tightly  depend  on  time 
synchronization.  Accurate  time  synchronization  method  is  the 
premise  to  guarantee  all  nodes  switch  between  sleep  state  and 
active  state  simultaneously. 

In  this  paper,  we  propose  a  novel  contention-based  MAC 
protocol:  asynchronous  contention-based  and  energy-efficient 
MAC  protocol  (ACEMAC).  Our  algorithm  overcomes  the 
disadvantage  for  S-MAC  and  T-MAC,  tightly  depending  on 
time  synchronization^],  [7],  [8],  [9].  ACEMAC  sets  up  system 
time  through  free-run  timing  method.  At  the  same  time,  clock 
drift  aware  method  is  used  to  solve  schedule  un-coincidence  of 
different  sensor  nodes.  And  no  extra  control  channel  is  needed. 
We  use  node  to  stand  for  sensor  node  in  the  following  parts. 

The  remainder  of  this  paper  is  organized  as  follows.  In  sec¬ 
tion  II,  we  provide  some  assumptions  and  the  network  model 
of  our  algorithm.  Our  algorithm  is  described  in  section  III. 
Average  latency  is  analyzed  in  section  IV.  Simulation  results 
are  given  in  Section  V.And  Section  VI  concludes  this  paper. 


II.  Assumptions  and  Network  Model 

We  have  following  assumptions  related  to  our  algorithm: 

•  There  are  only  two  cases  that  can  cause  the  packet 
transmission  fail.  One  is  collision  between  packets.  The 
other  is  that  the  receiver  is  unavailable  for  receiving  when 
transmitter  sending  packets; 

•  The  errors  caused  by  channel  can  be  recovered  by  channel 
coding  at  physical  layer; 

•  We  can  use  certain  clustering  algorithm,  such  as 
LEACH[16],  to  cluster  a  whole  WSNs  into  independent 
groups; 

•  The  traffic  for  each  node  is  the  information  detected  and 
collected  by  this  node; 

•  The  information  transmitted  over  the  network  is  cluster- 
inside; 

•  For  information  transmission  time,  we  ignore  the  process 
time  at  receivers  and  propagation  time. 

The  network  model  that  fits  our  algorithm  is  described  as 
follows.  For  one  cluster,  there  are  some  normal  nodes  and  one 
cluster  head.  In  the  same  cluster,  there  is  only  one  hop  between 
cluster  head  and  normal  nodes.  If  two  normal  nodes  stay 
during  each  other’s  communication  range,  they  are  neighbors 
for  each  other.  All  normal  nodes  stay  within  cluster  head’s 
communication  range.  But  for  different  normal  nodes,  they 
may  not  be  neighbors  with  each  other.  Therefore,  within  one 
cluster,  there  can  be  several  nodes  making  transmission  at  the 
same  time. 

III.  Asynchronous  Contention-Based  and 
Energy-Efficient  MAC  Protocol  (ACEMAC) 

We  have  two  goals  when  designing  our  algorithm.  One 
is  to  implement  energy  reservation  through  scheduling  all 
normal  nodes  powering  off/on  simultaneously.  The  other  one 
is  to  obtain  an  optimal  way  to  design  the  duration  for  power 
off/on.  We  exploit  re-schedule  method,  a  simple  and  easy 
implement  way,  to  archive  the  coincidence  among  normal 
nodes  instead  of  time  synchronization  mechanism.  In  order 
to  resolve  the  problems  caused  by  time  variant  environment, 
network  topology  and  traffic  strength  of  nodes,  we  design  a 
clock  drift  aware  system  to  dynamically  and  adaptively  adjust 
the  re-schedule  duty-cycle.  Our  method  belongs  to  contention- 
based  MAC  protocol  category.  In  the  following  parts,  we  will 
discuss  the  detail  about  our  algorithm. 

A.  System  Time  Scheme 

Fig.  1  presents  the  system  time  scheme  structure.  ACEMAC 


Fig.  1.  System  Time  Scheme  Structure 

* 

divides  system  time  into  different  stages.  During  different 


stages,  node  operates  differently.  The  function  for  each  stage 
is: 

•  ’’Hello  Stage”  is  used  by  normal  nodes  to  exchange 
HELLO  message; 

•  ’’Schedule  Broadcast  Stage”  is  used  by  cluster  head  to 
broadcast  schedules  among  cluster; 

•  ”On  Stage”  is  preserved  for  all  nodes  to  power  on 
communication  interface.  Therefore,  transmitting  and  re¬ 
ceiving  can  be  implemented; 

•  ’’Off  Stage”  is  preserved  for  all  nodes  to  power  off  com¬ 
munication  interface.  Therefore,  there  is  no  transmitting 
and  receiving  happened.  But  data  storage  and  sensing  still 
continue; 

From  above  figure,  we  can  see  that  there  are  several  rotations 
from  power  on  to  power  off  within  the  interval  (T)  between 
two  continuous  schedule  broadcast. 

B.  System  Schedule  Set  Up 

ACEMAC  schedules  all  nodes’  transmissions  into  certain 
stages  of  system  operation  time.  In  other  time,  nodes  reserve 
energy  through  powering  off  their  RF  interfaces.  The  premise 
for  ACEMAC  work  smoothly  is  to  set  up  a  coincident  system 
schedule  correctly  and  effectively  among  nodes  within  on 
cluster. 

Setting  up  a  global  clock  within  one  group,  as  done  for  most 
TDMA  system,  is  not  a  good  choice  for  WSNs.  The  reasons 
are: 

•  Spectral  efficiency  requires  the  overhead  brought  by 
exchanging  information  among  nodes  for  synchronizing 
the  time  of  different  nodes  should  be  as  little  as  possible; 

•  Energy  efficiency  brings  forward  cutting  short  the  energy 
consumed  by  time  synchronization  processing; 

•  Limitation  on  computation  and  memory  resources  make 
it  impossible  to  choose  high  quality,  but  complex  time 
synchronization  methods; 

Our  algorithm  can  overcome  this  defect  of  setting  up  global 
clock.  We  choose  free-run  timing  method.  That  means,  each 
node  simply  times  itself  to  its  own  free-running  local  oscillator. 
In  order  to  guarantee  each  node  in  the  same  cluster  to  switch 
to  the  same  stage  simultaneously  basing  on  local  clock,  special 
schedule  broadcast  is  used. 

Cluster  head  makes  a  broadcast  to  inform  the  system 
schedule  to  every  normal  node  in  this  cluster.  In  Fig.  3, 
schedule  broadcast  message  format  is  shown.  In  this  schedule 

Type  Source  On-Duration  Off-Duration 

Fig.  2.  Schedule  Broadcast  Packet  Format 

broadcast  message,  On- Duration  field  tells  nodes  the  length 
of  time  after  which  they  should  end  transmitting  and  receiving; 
Of f— Duration  field  tells  nodes  the  length  of  time  after 
which  they  should  end  current  power  on-off  rotation. 

From  above  figure,  we  can  see  that  a  schedule  broadcast 
packet  regulates  the  duration  of  different  stages,  instead  of  the 


exact  start  time  and  end  time.  Even  though  each  node  uses  its 
own  local  clock,  basing  on  this  schedule  broadcast,  all  nodes 
could  enter  the  same  working  state  simultaneously  according 
to  our  assumption. 

Overcoming  clock  drift  is  another  function  for  schedule 
broadcast,  besides  informing  nodes  system  schedule.  The 
quality  of  each  local  clock  usually  boils  down  to  its  frequency 
stability-that  is,  the  ability  of  its  frequency  standard  to  emit 
events  at  a  constant  frequency  over  time.  In  general,  as 
frequency  standards’  stability  and  accuracy  increase,  so  do 
their  power  requirements,  size  and  cost.  For  general  WSNs, 
whose  nodes  are  usually  configured  simply  and  inexpensively. 
Therefore,  there  are  always  some  clock  shifts  caused  by  fre¬ 
quency  standard  inaccurate  for  nodes.  Moreover,  the  frequency 
generated  by  a  quartz  oscillator  is  also  affected  by  a  number 
of  environmental  factors:  the  voltage  applied  to  it,  the  ambient 
temperature,  acceleration  in  space,  magnetic  fields,  and  so 
forth.  More  subtle  effects  as  the  oscillator  ages  also  cause 
longer-term  frequency  changes.  Under  these  conditions,  after 
a  period  of  time,  system  schedule  un-coincident  among  nodes 
is  unavoidable.  See  Fig.  3. 
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Fig.  3.  Schedule  Un-coincidence  Due  to  Clock  Drift 


If  the  frequency  standard  for  node  A  is  faster  than  that  of 
node  B.  From  above  figure,  we  can  see  that  the  transmitter, 
node  A,  precedes  the  receiver,  node  B,  into  the  on  stage.  So 
during  Ail,  data  transmissions  cannot  be  done  successfully 
between  node  A  and  node  B.  We  can  solve  this  problem 
through  doing  re-schedule.  See  Fig.  4.  After  each  re-schedule, 
Ail  between  node  A  and  node  B  is  removed. 
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of  coincident,  and  then  re-schedule  only  when  needed.  The  aim 
of  our  algorithm  is  to  reserve  energy  when  ensuring  certain 
successful  transmission  possibility. 

For  designing  re-schedule  duty-cycle  (i.e.,  the  length  of  the 
interval  of  two  continuous  schedule  broadcast),  we  have  two 
options.  One  is  dynamically  adjusting  re-schedule  duty-cycle 
method.  The  other  is  fixed  re-schedule  duty-cycle  scheme. 

It  is  obvious  that  for  fixed  duty-cycle,  there  is  a  large  amount 
of  energy  used  to  send  and  receive  the  schedule  broadcast,  if 
the  re-schedule  duty-cycle  is  short.  On  the  other  hand,  there 
is  a  great  deal  of  energy  wasted  by  fail  information  exchange, 
if  the  re-schedule  duty-cycle  is  long.  However,  when  most 
of  fail  transmissions  are  caused  by  clock  drift,  such  energy 
for  re-schedule  is  reasonable  and  worthy.  So  adaptively  adjust 
method  is  a  better  choice. 

We  build  a  kind  of  clock  drift  aware  scheme  to  dynamically 
adjust  re-schedule  broadcast  duty-cycle.  That  is,  when  the 
cluster-head  detects  the  clock  difference  among  nodes  influ¬ 
ences  the  information  transmission  at  very  small  degree,  an 
acceptable  one,  it  would  prolong  the  duty-cycle.  But  once  the 
cluster-head  finds  that  the  clock  drift  decreases  the  successful 
transmission  rate  sharply,  it  would  shorten  the  duty-cycle.  The 
details  are  discussed  in  the  following  part. 

Each  node  would  send  out  a  re-schedule  request  after  a 
fail  transmission.  In  this  paper,  when  the  retransmission  time 
exceeds  certain  value,  we  call  it  a  fail  transmission.  It  should 
be  noted  that  retransmission  is  only  for  data  packets.  We 
use  this  message  to  inform  the  cluster  head  about  the  fail 
transmission  ratio  for  this  node.  The  format  of  re-schedule 
request  packet  is  shown  in  Fig.5.  We  use  the  same  model  to 
determine  the  value  of  Nret  as  in  [13]. 

Nret  «  1.45/n— ^ —  (1) 

1  -  p 

Where  p  is  the  probability  of  successful  transmission. 


Type 


Source 


Transmission  Fail  Rate 


Fig.  5.  Re-schedule  Request  Packet  Format 


Fig.  4.  Schedule  Un-coincidence  Removing  Through  Re-schedule 

C.  Reschedule  Duty-Cycle  Design 

From  above  discussion,  we  notice  that  keeping  system 
schedule  coincident  among  nodes  can  avoid  transmission  fail¬ 
ures  caused  by  receivers  do  not  power  on  their  RF  interface. 
However,  providing  state  switching  simultaneously  at  all  time 
will  is  not  necessary.  For  example,  there  are  two  nodes.  Their 
frequency  standards  have  great  discrepancy.  So  there  is  quite 
time  difference  between  them.  But  there  is  seldom  informa¬ 
tion  exchange  between  them.  This  kind  of  system  schedule 
un-coincidence  almost  has  no  any  influence  on  information 
transmission.  Therefore,  some  nodes  can  be  allowed  to  go  out 


According  to  the  re-schedule  requests  received,  cluster 
head  adjusts  the  next  re-schedule  duty-cycle  before  doing  re- 


schedule  broadcast.  That  is 

(  2  x  Tj-i , 

when  N  decrease ; 

Ti  =  U  x  7U , 

when  N  increase ;  (2) 

[Ti-i , 

otherwise . 

where  T*  is  the  next  duty-cycle  of  re-schedule  broadcast;  T •_! 

is  the  previous  duty-cycle;  and  N  is  the  number  for  cluster 
head  shortening  re-schedule  duty-cycle.  The  range  of  N  is  [0, 
NmaxY  If  half  number  of  nodes’  fail  transmission  ratio  are 
bigger  than  5%,  cluster-head  during  this  period  of  time,  it 
could  add  1  to  the  last  value  of  N.  Otherwise,  the  cluster-head 
should  reduce  the  last  value  of  N  by  1  until  N=Nmax, 


D.  HELLO  Message  Design 

Hello  stage  is  preserved  for  each  node  to  broadcast  its 
own  HELLO  message.  Through  this  message,  each  node  can 
inform  neighbors  that  it  is  alive.  Moreover,  this  message 
includes  each  node’s  traffic  arrival  rate  (A  packets/second) 
and  average  service  time(ju  packets/second).  This  information 
is  the  considering  factor  for  cluster-head  to  determine  the 
length  of  On-stage  and  Off-stage.  Fig.  2  shows  the  structure 
of  a  HELLO  message.  The  traffic  arrival  rate  and  average 


Type 

Source 

Data  Arrival  Rate 

Service  Rate 

Fig.  6.  HELLO  Message  Format 


service  time  are  statistic  values  basing  on  previous  traffic 
arrival  rate  and  average  service  time.  In  our  algorithm,  we 
define  a  new  terminology:  coherent  time.  Coherent  time  is 
a  special  length  of  period,  within  which  we  assume  that 
the  environment  condition,  the  channel  condition,  network 
topology  and  traffic  strength  are  fixed.  During  every  coherent 
time,  we  carry  on  measurement  on  traffic  arrival  rate  and 
service  time.  In  order  to  reduce  the  dependent  degree  on  the 
cluster  head,  the  estimations  for  traffic  arrival  rate  and  average 
service  time  are  done  by  each  normal  node  locally.  Each  node 
randomly  chooses  a  time  within  the  hello  stage  to  send  out 
HELLO  message.  We  make  this  random  process  complies 
with  uniform  distribution.  Transmission  time’s  randomicity 
could  reduce  the  collision  probability  for  HELLO  message 
and  increase  the  successful  transmission  probability.  In  the 
following  part,  we  will  analyze  the  probability  ( Ps )  for  a 
normal  node  send  HELLO  message  successfully. 


N 

Ps  =  l-'E,P{n  =  i}  (3) 

i= 2 

where  n  is  the  number  of  nodes  whose  transmission  is  fail  for 
collision.  There  are  N  normal  nodes  in  this  cluster. 
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where  r  is  the  transmission  duration  for  a  HELLO  message. 
Th  denotes  the  duration  of  hello  stage.  Then  the  probability  for 
all  N  nodes  send  message  to  cluster  head  successfully  during 
HELLO  state(T/l)  is: 


N 

p*  =  i-E 


i=2  L 


r(2Th  -  r) 


n 


(5) 


We  can  see  that  P$  is  a  function  of  parameter  N  and  TV 


E.  Power  On  Duration  (Tn)  Design 

For  WSNs,  according  to  the  application  field,  there  are  lots 
of  traffic  types.  We  don’t  hope  to  limit  our  parameter  design 
into  one  specific  application.  So  basing  on  the  traffic  arrival 
rate,  we  design  the  sleep  duration.  If  we  know  the  average 


traffic,  arrival  rate  is  A*  for  node  i.  So  the  average  number 
(Ni)  of  packets  arriving  during  the  sleep  part  for  node  i  is: 

Ni  =  A  fat  (6) 

According  to  (1),  the  value  of  Ni  is  proportional  to  TSti.  We 
know  that  the  buffer  size  for  each  sensor  node  is  limited.  When 
the  buffer  is  used  out,  the  following  coming  data  packets  must 
be  discarded.  We  can  obtain  one  criterion  for  determining  the 
value  of  Ts^:  the  number  of  arrival  messages  should  not  be 
bigger  than  the  buffer  size  during  the  Ts>i  period,  i.e., 

Ni  <  K.  (7) 

where  K  is  the  available  buffer  size.  We  assume  that  for  each 
node,  the  buffer  size  is  same.  It  is  reasonable  and  realistic  to 
suppose  each  node  equipped  with  same  resource. 

Since  Ni~X{Ts^.  Then  we  have 

A iT,ti  <  K  —>  TS'i  <  ~(i  =  1,2,  •  •  •  ,M)  (8) 

Aj 

We  assume  that  there  are  M  nodes  in  a  cluster  totally.  The 
data  arrival  rate  for  every  node  of  this  cluster  is  identical  and 
independent.  So,  for  the  same  number  of  data  packets  arriving, 
the  duration  for  different  node  might  be  different.  We  let  the 
Cluster-head  choose  the  shortest  duration  for  the  whole  cluster. 
The  sleeping  duration  for  whole  cluster  is: 


According  to  Little  theorem[15],  the  average  number^)  of 
customers  in  a  queue  system  equals  to  the  product  of  customer 
arrival  rate(A)  and  average  time(u;)  spent  in  the  system: 

Nq  =  Ait;  (10) 

For  node  i,  we  can  get: 

Ni  —  XiWi  (11) 

Compare  (5)  and  (10),  we  get 

At  (TSji  +  Tb)  =  A iWi  — ►  TSij  +  fB=Wi  (12) 

From  this  equation,  we  get  another  criterion  for  choosing  the 
value  of  TSji.  For  all  sensor  nodes,  the  permitted  latency  is 
Wmax.  So  we  can  get: 

Ts,t  4-  Tb  <  Wmax  ^Ts<  Wmax  -  Tb  (13) 

From  (8)  and  (12),  the  optimum  value  of  Ts  should  follow  the 
rule: 

Ts  <  min  |wmoi,mm(^)|  (14) 

In  this  paper,  we  set  Ts  as: 

Ts  =  min  |wmai,min(^)|  (15) 

We  know  that,  the  longer  the  sleeping  period  is,  the  more 
energy  reserved.  However,  the  disadvantage  is  the  average 
latency  of  data  packets  increasing.  Energy  is  one  of  the  most 
important  resources  for  sensor  networks.  Reserving  energy  and 
reduce  the  latency  are  two  contradict  tasks. 


F.  Power  Off  Duration(Tf)  Design 

During  active  part,  nodes  start  to  send  out  data  packets, 
which  are  received  during  last  sleeping  part  and  current 
active  part.  Through  contention,  each  node  obtains  the  right 
to  hold  the  link.  The  contention  progress  is  almost  like 
CSAMA/CA[1].  Using  carrier  sense,  RTS/CTS  shake-hand 
and  ACK  confirmation,  CSMA/CA  ensures  no  collision  for 
data  packet  transmission.  And  also,  it  permits  the  transmitter  to 
know  exactly  that  if  the  data  packet  transmission  is  successful 
or  not.  In  our  sleeping  duration  and  active  duration  designing, 
we  not  only  try  to  extend  the  sleeping  time  to  reserve  energy, 
but  also  ensure  data  packets’  latency  are  not  longer  than 
acceptable  value.  So,  the  criteria  for  active  duration  is  that 
it  (Ta,i)  is  long  enough  for  all  received  data  packets  to 
be  sent  out.  Now,  we  should  estimate  the  average  number 
(Nt,i)  of  data  packets  received  during  previous  sleep  part  and 
current  active  part,  then  we  can  determine  the  value  of  Ta>i. 
Because  the  traffic  arrival  process  is  independent  with  the  data 
transmission  process,  for  Nt,i  we  have 

Tt,i  =  Ns,i  4-  Na>i  =  A  i(Ts  +  Ta)  (16) 

Basing  on  the  Ntj  data  packets  should  be  sent  out  during 
current  active  duration(Ta^),  the  following  equation  existed: 

=  A  i(Ta,i  +  7s,i)  (17) 


model.  In  this  part,  we  analyze  the  average  latency  for  one 
M/G/l  system.  The  average  waiting  time  ( W )  is: 


W  = 


R 

1  -  P 


(20) 


where  R  is  the  average  residual  time,  p  is  the  traffic  intensity. 
There  are 


P  = 


A 

P 


(21) 


and 


_  AX2  (1  -  p)V2 
2  2V 


(22) 


where  A  is  the  data  arrival  rate;  p  is  the  average  service  time. 
X 2  stands  for  second  moment  of  service  time,  X  stands  for 
average  service  time  and  V  stands  for  the  average  sleeping 
time,  p  is  known  if  given  A  and  p.  In  order  to  get  the  value 
of  W ,  we  should  compute  X2.  Since  X=Tb(&)+£s.  Then 


X2  =  E[X2} 

=  E[T2B{k)}  +  t2s  +  2 tsE[TB(k))  (23) 


ts  is  the  average  time  that  the  link  is  sensed  busy  due  to  a 
successful  transmission.  Basing  on  [11],  we  have 


The  active  duration  for  node  i  is: 


pi  —  A  i 


(18) 


If  the  cluster-head  chooses  the  longest  duration  as  the  whole 
cluster’s  active  duration,  for  all  nodes,  each  node  could  have 
enough  time  to  send  all  received  data  packets  out.  Then  we 
get  the  active  duration  for  whole  cluster.  That  is, 


IV.  System  Performance  Theory  Analysis 


Contention-based  MAC  protocol  divides  the  system  time 
into  two  types:  sleep  time  and  active  time.  For  sleep  time, 
each  node  can  be  regarded  as  a  waiting  queue.  Data  packets 
independently  enter  this  queue  at  certain  rate.  But  there  is  no 
any  transmission.  For  active  time,  a  group  of  nodes,  locating 
in  each  other’s  communication  range,  could  be  treated  as 
a  M/M/1  queuing  system.  Within  one  cluster,  all  normal 
nodes  stay  within  cluster-head’s  communication  range.  But 
for  different  normal  nodes,  maybe  they  are  not  neighbors. 
That  depends  on  the  topology  of  the  cluster.  So,  the  link 
could  be  successfully  used  by  more  than  one  transmitting- 
receiving  pairs  at  one  time.  For  this  special  characteristic, 
within  one  cluster,  there  are  several  M/M/1  queuing  systems 
existed  simultaneously.  Looking  through  the  system  time,  we 
cannot  use  M/M/1  queue  model  any  more  for  introducing  sleep 
time.  We  use  M/G/l  Queue  with  vocation  model[12]  as  our 


ft  +  (*-!)«.  <M> 

i=l 

For  E[T2(k)\,  we  have: 

oo 

Efflm  =  Y^fB)\k)p{K=k) 

k=l 

=  fX£f(^-i))2  +  (*-i  )2(tc)2 

fc=i  i= i  z 

k 

+2<c(fc-l)^|(Wi-l)](l-g)fc-\?5) 


i=l 


We  let 


E[TUk)]  =  -(S1  +  S2) 
For  Su  we  have: 


m  k 

k=  1  i=  1 


(26) 


(27) 


For  S2,  we  have: 


oo  m 


S*  =  E  ~  1))2  +  (  E  &  ~  1))2K1  -  9)fc_19  (28) 


fc=77l-fl  1  =  1 


i=m+ 1 


where  m  is  the  ’’maximum  back-off  stage”.  Wi  is  the  minimum 
contention  window  size. 


V.  Simulations  and  Performance  Evaluation 


In  this  section  we  evaluate  the  performance  of  our  method. 
OPNET  is  used  to  run  simulators  on  the  network.  And  the 
size  of  the  network  varies  from  11  nodes  to  22  nodes.  All 
nodes  transmit  to  some  other  node  in  the  network  accord¬ 
ing  to  the  same  source  rate  with  fixed  packet  sizes.  The 
communication  range  for  each  node  is  20  meters.  Nodes  are 
randomly  placed  in  an  area  of  100x100  meters  and  have  no 
mobility.  We  trace  each  node  in  the  network  and  compute  the 
successful  transmission  rate,  data  packet  average  time  delay 
and  the  time  of  first  node  dead.  The  experiment  is  repeated 
for  100  different  seeds.  We  do  that  not  only  for  statistical 
reasons,  but  also  for  the  fairness  problem  inherent  in  the 
DCF  of  MAC  protocol.  Regarding  the  physical  layer,  we  use 
Frequency-Hopping  Spread  Spectrum  (FHSS)  with  a  raw  bit 
rate  of  1Mbps.  Other  physical  layer  parameters  are  same  as 
IEEE802.il  protocol  10]. 

We  separately  use  our  method  and  a  derivate  way,  which 
only  exchanges  the  sleeping  part  and  active  part  of  a  frame. 
Under  different  clock  shift  rate,  we  compare  the  system 
successful  transmission  rate,  average  latency  and  the  time  of 
first  node  dead.  From  Fig.  7,  we  can  see  that  with  the  clock 


Fig.  7.  Data  Packet  Successful  Transmission  Rate 

shift  rate  increasing,  there  are  more  data  packets  transmission 
failed.  When  clock  shift  rate  is  larger  than  0.005(i.e.,  there 
are  0.005  seconds  drift  every  second),  sleep-active  method 
(i.e.,  for  each  frame,  the  active  part  follows  sleeping  part),  the 
performances  decrease  sharply.  The  successful  transmission 
rate  decreases  about  15.39%.  But  for  active-sleep  method, 
there  is  only  5.49%  decrease.  Our  algorithm,  active-sleep 
method,  has  better  tolerance  for  clock  shift.  Fig.  8,  shows 
that  with  the  clock  shift  rate  increases,  the  average  time  delay 
for  data  packet  decreases.  The  reason  is,  we  just  consider  the 
successful  transmission.  And  a  data  packet  would  be  discarded 
after  several  unsuccessful  retransmissions.  The  average  latency 
of  the  active-sleep  method  is  almost  33.3%  shorter  than  sleep- 
active  method.  From  Fig.  9,  we  can  see  that  the  time  of  first 
node  dead  is  brought  forward  with  the  clock  shift  rate  increase. 
These  simulation  results  also  show  that  the  performance  of  our 
system  is  not  sensitive  with  the  network  density  very  much. 


Fig.  8.  Average  Time  Delay  for  Data  Packet 


Fig.  9.  Time  of  First  Node  Dead 


Fig.  10,  Fig.  11,  Fig.  12  show  the  relationship  between 
different  sleeping  duration  and  system  performance.  We  use 
optimum,  longer  and  shorter  sleeping  duration  values  sepa¬ 
rately.  The  sleeping  duration,  basing  on  our  algorithm,  let  the 
successful  transmission  rate  is  almost  same  as  all  nodes  always 
working  without  any  sleep.  Moreover,  it  prolongs  the  time  of 
first  node  dead  about  22.58%.  We  also  can  see  that  with  the 
sleeping  duration  increase,  the  average  latency  increase. 

We  do  the  same  experiments  with  the  active  duration.  From 
Fig.  13  and  Fig.  15,  we  can  get  that,  even  the  successful 
transmission  rate  for  optimum  active  duration  is  a  little  lower 
than  no  sleep  one,  the  time  of  the  first  node  dead  is  prolonged 
about  22.23%.  More  energy  are  reserved. 

In  the  last  simulation,  we  do  experiments  with  fixed  duty- 
cycle  and  adaptively  adjusted  duty-cycle.  There  are  two  cases: 
one  is  the  initial  synchronization  broadcast  duty-cycle  that  is 
2  times  of  frame  time;  the  other  one  is  the  duty-cycle  that  is 
the  4  times  of  frame  time.  Under  different  clock  shift  rate,  we 
compare  the  successful  transmission  rate. 

VI.  Conclusions 

In  this  paper,  we  propose  a  novel  energy-efficient  MAC  pro¬ 
tocol:  ACEMAC.  Our  algorithm  is  an  asynchronous  method, 
which  completely  overcomes  the  biggest  disadvantage,  i.e., 
so  tightly  depending  on  time  synchronization,  of  existing 


Fig.  1 1 .  Average  Time  Delay  for  Data  Pacet 


Fig.  13.  Data  Packet  Successful  Transmission  Rate 


contention-based  and  energy-efficient  MAC  protocols.  We  use 
schedule  broadcast  and  clock  drift  aware  scheme  to  dynam¬ 
ically  solve  the  system  schedule  un-coincidence  problem.  In 
this  paper,  we  also  supply  a  parameter  designing  model  to 
optimize  key  parameters.  Moreover,  we  make  theory  analysis 
on  average  latency  (Point-to-Point)  for  our  algorithm.  The 
simulation  results  prove  that  our  algorithm  can  successfully 
obtain  optimum  re-schedule  duty-cycle,  power  off  duration 
and  power  on  duration  to  achieve  higher  average  successful 
transmission  ratio,  shorter  data  packet  average  time  delay  and 
longer  network  lifetime. 
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Fig.  14.  Average  Time  Delay  for  Data  Pacet 


Fig.  15.  Time  of  First  Node  Dead 


Fig.  16.  Data  Packet  Successful  Transmission  Rate 


